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In previous comparisons of the theoretical and experimental temperature effect for the 
reflection of x-rays, a marked discrepancy was found at high temperatures. It is here shown 


that this discrepancy disappears when proper account is taken of the variation of the charac- 


teristic temperature with temperature. Further, a relationship is found between the temperature 


factors for two allotropic modifications of the same element. 


$1. INTRODUCTION 


HE effect of temperature in decreasing the 

intensity of reflected x-rays has been sur- 
prisingly well interpreted at ordinary tempera- 
tures by the simple Debye model of a solid." 
In this model the observed structure factor F is 
given by fe-”, where f is the structure factor for 
a perfect lattice, and!:* 


6h*T sin? @ sO 
M=-— — -0(=). (1) 
m,k@Q? 2 T 


Here m, is the mean mass of the atoms, @ is the 
glancing angle of incidence of the x-rays on 
the reflecting plane, and ¢ is a factor which 
deviates appreciably from unity only when 


* The senior author was aided in part by a grant from 
the Rockefeller Foundation to Washington University for 
research in science. 

1See Compton and Allison, X-Rays in 
Experiment, pp. 435-445. 

* For a simple derivation, see Zener and Jauncey, Phys. 
Rev. 49, 17 (1936). 
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©/T>1. If © is regarded as independent of 
temperature, as it has always been in the litera- 
ture on x-rays, log F becomes a linear function of 
temperature in the region 7 > ©. The only experi- 
mental values of F at high temperatures are for 
NaCl and for KCI.4 The observed values of log F 
for these two crystals deviates radically from a 
linear function above 400°C. The purpose of the 
first part of this paper is to show that the tem- 
perature factor e~” calculated by the simple 
Debye model is in agreement with the observed 
temperature factor even at high temperatures 
when proper account is taken of the variation 
of © with temperature. 

While in §2 the continuous variation of M for 
a single crystal phase is discussed for the entire 
temperature range of that phase, in §3 the dis- 
continuous variation of M is examined when 
the crystal undergoes an allotropic transfor- 
mation. 


3 James and Firth, Proc. Roy. Soc. A117, 62 (1928). 
‘ James and Brindley, Proc. Roy. Soc. A121, 155 (1928). 
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$2. TEMPERATURE FACTOR AT HIGH 
TEMPERATURES 

The Debye model of a solid can make definite 
predictions about the temperature dependence of 
its properties only when the volume of the solid 
is kept constant. Thus in order to compare the 
observed specific heat with that predicted by 
the Debye model, C, is always first reduced to 
Cy. In order fairly to compare with experiment 
the prediction of the Debye model as to the in- 
fluence of temperature upon the intensity of 
reflected x-rays, proper account must be taken 
of the fact that the observations are made at 
constant pressure and not at constant volume. 

Two methods of approach are available. In the 
first method, Eq. (1) is written in the form! ? 


M=8r'y sin? 6/22, (2) 


where yu is the mean-square displacement of an 
atom. One then constructs a molecular model of 
the solid, including third- and fourth-order terms 
in the expression for the potential energy, which 
expands with temperature. Both C,—Cy and the 
mean-square displacement yu are then directly 
calculated. This approach has already been 
investigated,®> but has not given any quantita- 
tive predictions on the variation from the linear 
dependence of M with T at high temperatures. 

In the second method of approach one starts 
from the phenomenological Debye model, and 
determines that dependence of © upon tempera- 
ture which is necessary to fit the specific heat 
data. This function 0(T) is then substituted 
into Eq. (1). This second method is the one 
adopted in this paper. 

Two methods may be used to determine the 
function ©(7) from specified heat data. The 
simpler one is to:find that 0 which gives the ob- 
served Cy.®° Thus 0(7) is defined by the equation 


Cy(T)=Cp{ O(T)/T}. (3) 


The second method’ is to determine that func- 
tion ©@(7) which makes the observed entropy 
equal to the entropy calculated from the Debye 
model. The observed entropy (per atom) at a 
temperature 7» is 


51, Waller, Ann. d. Physik 83, 153 (1927). 

6 Jauncey and Deming, Phys. Rev. 48, 481 (1935). 

7 This method is in principle identical to that used by 
Fowler, Statistical Mechanics. The derivation here given is, 
however, simpler than that given by Fowler. 
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eC, 
Sovs.(Ta) = f aT. 

J» : i 
In order to calculate the entropy by the Debye 
model, one must pass from the absolute zero to 
the temperature 7) in two steps. First one 
changes the volume of the solid, keeping it at 
the absolute zero, from its value at p=0, T=0, 
to its value at p=0, T=T>. According to the 
Nernst heat theorem,’ the entropy change is 
zero in this change of volume at the absolute 
zero. The solid is now heated, at constant 
‘=7>). Throughout this 
second process the characteristic temperature 


volume, from T=0 to 


remains constant. Hence the entropy change is 
given by 

> Cy{O(To)/T} _ 

So(T) =k | : dT. 
all 1 

If we now equate these two expressions for the 
entropy, and differentiate with respect to 7», 
we obtain 


(d In ©)/dT = —(C,—Cp(9/T))/TCp(0/T). (4) 


If the Debye model gave an exact representa- 
tion of the specific heat, and if the empirical C, 
and Cy were known precisely, Eqs. (3) and (4) 
would give the same function ©0(7). However, 
since Cp(Q/T) becomes very insensitive to © 
when JT becomes comparable with and larger 
than ©, either a small error in the empirical 
values of Cy or a slight deviation of the true 
Cy from Cp produces a relatively large change in 
©(T). This is illustrated in Table I. The function 
©(T) obtained from Eq. (4) is not sensitive to 
these errors. The second method is thus used in 
this paper. 

The integration of Eq. (4) only gives the ratio 
©(T)/0(0). O(0) must be determined by using 
Eq. (3) at low temperatures. The function 0(7) 
thus determined, when substituted into Eq. (1), 
gives the temperature factor at all temperatures 

TABLE I. The second row gives the percent error in © 


due to a 1 percent error in Cy, when ©@(7) is determined 
from Eq. (3). 


0/T 0.25 0.5 1 2 4 


5 ln @/6 ln Cy | 240 35 95 2.7 0.8 
| - 


8 My Plank, Treatise on Thermodynamics, 5th edition, 
74 
p. 274. 
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ras_e Il. Log of ratio @(T)/O(0) for KCl and 


T (% —200 ~ 100 0 
—In {O(T)/O(0)! 
KCI | 0.00346 0.0176 0.0341 


NaCl | 0.00266 0.0160 0.0330 


with constant pressure as predicted by the 
Debye model. 

According to the Debye model Cp and Cy are 
identical. Hence in Eq. (4) Cp may be replaced 
by Cy. Now it is well known that the ratio® 
(C,—Cy)/(TCyC,) is very nearly a constant 
over the entire temperature range of those 
solids with a simple structure, the only solids to 
which the Debye model is applicable. The con- 
stant K may be calculated from the known values 
of C, and Cy at one temperature To, say at room 
temperature. When this relationship is used, 
the integration of Eq. (4) gives 

O(T) =0(0)e-*4 (7), (5) 
aT 
where H(T)= C,dT 


“0 


and K=[C,(To)—Cy(To) ]/TC,(To)Cv(T»). 

The only experiments over a wide range of 
temperatures are for NaCl* and KCI.‘ The ratio 
©(T)/O(0) obtained from Eq. (5) for these two 
crystals is given in Table II. The accepted values 
of ©(0), which are obtained from Eq. (3), are 
281°K for NaCl and 230°K for KCl. The com- 
parison between the calculated and _ experi- 
mental values of log F(T)/F(290°) is given in 
Figs. 1 and 2. A better fit with the experimental 
data is obtained by taking slightly larger values 
for ©(0) than given by Eq. (3). The experimental 
points at the two highest temperatures for KCI 
fall considerably below the best theoretical curve. 
In view of the excellent fit with the other points, 
and in view of the increasing spread of the 
experimental points with temperature, these two 
isolated experimental points are probably too 
low. 

It is to be expected that the (0) which gives 
the best fit to the x-ray reflection data will 
always be somewhat larger than the 0(0) de- 
termined from the specific heat data via Eq. (3). 
This is because on the Debye model the solid is 


° Eucken, Handbuch der Experimentalphysik, Vol. VIII 
1 (1929), p. 204. 
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NaCl at various temperatures as given by Eq. (5 

100 200 300 400 500 600 
0.0513 0.0690 0.0870 0.106 0.123 0.141 
0.0519 0.0717 0.0911 0.111 0.132 0.153 


strictly not specified by a single 0. The 0 used 
in Cp and in Eq. (1) refers to an average 0, 
but not to the same average. In the first case © 
refers to an average of the type (@-*)~!, in the 
second case to an average of the type (@-*)~!. 
Now it may readily be proved, by using Schwarz’ 
inequality,’ that 


The application of this inequality to an elasti- 
cally isotropic solid will now be considered. 
Here © as defined by the first average, say Oc, 
is given by 
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Fic. 1. Experimental points are marked by circles. 
Curve A, theoretical curve with constant @. Curves C 
and B, theoretical curves with © calculated by Eq. (5), 
with 0(0) = 281, O(0) = 298, respectively. 
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Fic. 2. Experimental points are marked by circles. 
Curve A, theoretical curve with constant ©. Curves C 
and B, theoretical curves with © calculated by Eq. (5), 
with @(0) = 230, @(0) = 240, respectively. 


1° Courant and Hilbert, Methoden der Mathematischen 
Phystk Vol. I (1924), pp. 2, 33. 
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TABLE III. The second row gives for various values of 
Poisson's ratio the ratio for an elastically isotropic solid 
of the O which must be used in the Debye factor, Eq. (1), 
for the intensity of x-ray reflection to the @ which must 
be used in the Debye expression for the specific heat. 


o 0.25 0.3 0.35 0.4 0.45 0.5 


0,/9c, 1.025 1.027 1.032 1.040 1.051 1.070 


“ 


where 6@;, 6, refer to the characteristic tempera- 
tures associated with the transverse and longi- 
tudinal vibrations, respectively. The © defined 
by the second average, say ©,, is given by 

@,-? = (26,2+6,-2) /3. 


P 
The ratio 0,/Oc, will in this case be a function 
only of Poisson’s ratio o, and is given in Table 
III. The elastic anisotropy of all actual crystals 
will still further increase the ratio of the two 


types of average. 


$3. RELATIONSHIP BETWEEN THE TEMPERATURE 
FACTORS OF ALLOTROPIC FORMS OF 
THE SAME ELEMENT 


According to the Debye model of a solid, each 
crystal phase is associated with a definite charac- 
teristic temperature (which according to §2 
must, however, be regarded as a function of 
temperature). The value of the characteristic 
temperature ©, of the low temperature modifica- 
tion at the transition temperature 7) may be 
calculated as described in §2. The value of the 
characteristic temperature ©, of the high tem- 
perature modification at ZT) may at once be 
obtained from the entropy equation 

Sp{ Or To} —Sp{O0i1/To} =L/T. (6) 
Here L is the latent heat of the transition (in 
calories per gram atom when S is the entropy 
per gram atom). The solution of this equation 
for ©, may be affected either by use of tables” 
of Sp(QO/T) or, when T>0, by the analytical 
expression" for Sp. When T>1.250, the follow- 
ing equation gives S with an error of less than 1 
percent: 

Sp=R{4+3 In T/O}. 

Substitution of this equation into (6) gives 


QO, >.= el (3RTo) | 


Using Eq. (1), we obtain 


4 An analytical expression for S(x) is given by Gruniesen, 
Handbuch der Physik, Vol. X (1926), 22. 

12 Extensive tables of Sp(x) are given by Simon, Hand- 
buch der Physik, Vol. X (1926), pp. 369-370. 


S. BILINSEY 


raB_e IV. 
HEAT Of 
TRANSITION TRANSFORMATION 
TEMPERATURI alories M 

ELEMENT ; per gram atom) My 
Tin! 298 564 1.54 
Arsenic? 1017 525 1.19 
Selenium® 423 174 1.15 
Zinc* 618 190 1.11 
M — 1108 158 1.05 
Manganese 1317 249 106 
Calcium® 673 100 1.05 
Thallium? 498 51 1.03 

1 Lewis and Randel, Thermodynamics (1923), p. 464 

2 Laschtschenko, J. Chem. Soc. 121, 972 (1922). 

Mellor, A Comprehensive Treatise on Inorganic and Theoretical 


Chemistry, Vol. X, p. 715. 
* Reference 3, Vol. IV, p. 430 
5 Umino Science Reports, Tohoku University 16, 775 (192 
6 International Critical Tables, Vol. VII, p. 293 
? Werner, Zeits. f. anorg. Chemie 83, 275 (1913 


Me/M,=e7!/@kTo (7) 


for the ratio of the two M’s at the transition 
temperature. 

When a crystal phase is anisotropic, the 
Debye function will vary with the orientation 
of the reflecting plane.” In this case Eq. (7) can 
only be valid when 4; and M, refer to some sort 
of average. The correct type of average may be 
uniquely obtained when 7>0. In this case the 
Einstein model of a solid may be used to de- 
scribe the vibrations of the atoms. In this model 
each atom is considered as bound elastically to a 
fixed lattice point. Let pis, iy, wiz ANd per, Mey, Mee 
be the mean squares of the displacement of the 
atoms along the three principal axes, for the low 
and high temperature phase, respectively. The 
change of entropy may then be written in terms 
of the ratio of two phase integrals." 


I SSS eM Fl wizt Flay 2) Dd xd ydz 
—=-—Rln . 
. SSIS COME wast Fl wayt 2°) e2adxd ydz 
The reduction of this integral gives 
MeréM2yKMe: Mizcbiyi2 =e? B2o, 
Using Eq. (2), we see that this equation is equiva- 
§ | | 1 

lent to Eq. (7) where now M stands for the 
geometrical mean (M’ M” M’”’)!. Here M’, M”, 
M’” are the three M’s in which the reflecting 
plane is normal to a principal axis. Table IV 
gives a list of the ratios M:/M, which are pre- 
dicted by Eq. (7). 


13 Zener, Phys. Rev. 49, 122 (1936). 
4 Herzfeld, Kinetische Theorie der Warme (1925), p. 161. 
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The Absorption of X-Rays of Wave-Length 50=\<150 X.U. by Elements of Low 
Atomic Number 


TREVOR R,. CUYKENDALL, Cornell Unive 


rsity, Ithaca, N. Y. 


Received May 11, 1936 


With a 600 kv x-ray source and a two-crystal spec- 
trometer especially designed for short wave-length x-rays, 
the mass absorption coefficients for C(6), Na(11), Al(13), 
S(16), K(19), Ni(28), Cu(29) were determined for wave- 
lengths from 50 to 209 X.U. A small 6A interval was used. 
rhe results are estimated to be in error by no more than 
two percent. The present results are compared with the 
few overlapping measurements of other observers. The 


INTRODUCTION 


HE determination of the absorption coeffi- 

cients of monochromatic x-rays in the 
region of wave-lengths shorter than 150 X.U. is 
of importance (1) in testing the Klein-Nishina 
formula! for the intensity lost by scattering, and 
(2) in providing additional information which 
may be used to check the validity of recent 
formulae for the photoelectric absorption. Fur- 
thermore, absorption measurements for mono- 
chromatic radiation are of value to the medical 
profession in the use of high voltage x-rays for 
therapeutic purposes. 

No systematic study of absorption coefficients 
has been made for the present wave-length 
region. Absorption measurements at wave- 
lengths less than 120 X.U. require a powerful, 
high voltage, x-ray source of constant intensity, 
and an analyzing spectrometer or other device 
which effects a satisfactory compromise between 
intensity and monochromaticity of the emergent 
beam. It is to be noted that there exists no 
practical target material, except possibly ura- 
nium, whose characteristic x-radiation is shorter 
than 120 X.U., and therefore the continuous 
spectrum must be used. 

Early determinations? of absorption coeffi- 
cients include a few measurements in the present 
wave-length region. However the low intensity 


! Klein and Nishina, Zeits. f. Physik 52, 853 (1928). 
* Compton and Allison, X-Rays in Theory and Experi- 
ment, p. 534. 
Richtmyer, Phys. Rev. 18, 13 (1921). 
* Richtmyer and Warburton, Phys. Rev. 22, 539 (1923). 
° Allen, Phys. Rev. 27, 266 (1926). 


values of the electronic scattering coefficient predicted by 
the Klein-Nishina formula are in excellent agreement 
with the experimental values of the scattering by carbon 
in the region 40=A=—100 X.U., which confirms the Klein 
Nishina formula for this wave-length region. The photo 
electric absorption in copper was found to be proportional 
toaA 


of the continuous spectrum produced by the 
x-ray sources then available, combined with 
the difficulties encountered in measuring very 
small ionization currents with the electrometer 
as compared with the FP-54 vacuum tube, 
favored intensity at the expense of resolving 
power. Thus, although a single crystal spec- 
trometer was used as monochromator, the 6A 
interval (as calculated from the spectrometrical 
geometry) was a large fraction of \, which 
means that the uncertainty in the effective wave- 
length was large. A further disadvantage of the 
single crystal spectrometer, resulting from the 
small Bragg angles encountered at short wave- 
lengths, is the problem of adequately shielding 
the ionization chamber from scattered radiation. 

The filtered end-radiation from a_hetero- 
geneous beam has been used recently in the de- 
termination of absorption coefficients.’ The 
results obtained with this method are discussed 
elsewhere.’ 

In consideration of the low resolving power 
used in previous work, and the small number of 
elements studied, the present investigation in 
the wave-length region 50=A<150 X.U. was 
undertaken with a 600 kv, 2000-watt, x-ray 
source, and with a _ two-crystal spectrometer 
especially designed for this wave-length region. 
The results of the absorption measurements for a 
series of elements of low atomic number are given 
herein, and for elements of high atomic number 
are given in the following paper by M. T. 
Jones.’ 


6 Alichanjan and Kosman, Zeits. f. Physik 90, 799 (1934). 
7 Jones, Phys. Rev. 50, 110 (1936). 
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APPARATUS AND METHOD OF OBSERVATION 


The apparatus used in this investigation has 
been previously described.’ Because of the 
greater practical resolving power without undue 
sacrifice of intensity, a two-crystal, rather than 
a single-crystal, spectrometer was used as mono- 
chromator. The particular design of the present 
instrument, which utilizes internal reflection of 
the x-rays as they are transmitted through the 
body of the crystals, permitted adequate shield- 
ing of the ionization chamber from scattered 
radiation. The vertical limiting slits were so 
adjusted that the geometrical resolving power 
was many times greater than that due to the 
crystals. The horizontal divergence was limited 
by the size of the focal spot. The 6d interval, as 
determined from considerations previously dis- 
cussed,*’ was 1.8<6\<1.0 X.U. for 50<\<215 
X.U., respectively. 

The absorbers, placed between the x-ray tube 
and the first crystal, were mounted as described 
in reference 8. Lead shields having holes 
3/4X1/2 inches in the path of the x-ray beam 
protected the absorbers from radiation scattered 
by other objects in the vicinity. It was thought 
possible that an error might be introduced in 
the observations taken at wave-length \,, due to 
radiation of \<),, present in the beam, under- 
going a wave-length shift by multiple Compton 
scattering within the absorber in such a way as 
to be reflected by the crystals as radiation of 
wave-length \;. The error due to this effect was 
calculated and found to be entirely negligible. 

Each measurement of an absorption coefficient 
consisted of 9 to 12 values of (J+J,.) and of 
(Io+J].), where I and IJ, represent, respectively, 
the ionization current with the absorber in the 
path of the beam, and removed, and where J, 
represents the background or residual ion 
current. At intervals during the recording of the 
other ion currents, 4 or 5 values of J, were de- 
termined by stopping the x-ray beam with a lead 
shutter, 3/4 inch thick, placed between the 
target and the first crystal. That this method 
gave the true residual current is shown by the 
fact that for the shutter open or closed no 
difference was observed in the residual ion 
current (J,) when one crystal had been turned a 
few degrees (about 5°) from its position pre- 


§ Cuykendall and Jones, Rev. Sci. Inst. 6, 356 (1935). 
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viously determined by the setting of the wave- 
length screw (see reference 8) for a given value 
of \. This test was made for several wave-lengths, 
including 50 X.U. 

At each wave-length \ the x-ray tube was 
operated at a voltage such that there was not 
present in the beam second-order energy of 
wave-length \/2. Especial care was taken to 
keep the tube voltage constant, but due to the 
type of voltage source employed it is probable 
that occasionally there were random fluctua- 
tions of as much as 3 percent. It is believed 
that this fluctuation was the greatest source of 
error. In order to reduce errors due to this cause, 
a large number of observations were taken. 


THE ABSORBERS INVESTIGATED 

Particular care was taken to obtain absorbers 
of high purity. The effect of such impurities as 
were present in the samples is less than the 
experimental error of the results. The carbon 
absorbers (99.9 percent pure) were regraphitized 
graphite made especially for spectroscopic pur- 
poses. The aluminum used was 99.97 percent 
pure and contained 0.008 percent copper and 
0.008 percent iron as the principal impurities. 
Chile copper, 99.96 percent pure, was used. 
The nickel absorbers were the same high purity 
absorbers used in a previous investigation.‘ The 
sulfur absorbers were carefully cast from the 
highest grade C.P. sulfur obtainable. 

The sodium and potassium absorbers were 
prepared by the R. and H. Chemicals Depart- 
ment of the E. I. duPont de Nemours & Com- 
pany. The metals were cast in blocks of suitable 
dimensions and each block was wrapped in sheet 
of aluminum foil whose edges were sealed with 
adhesive tape and paraffine. The purity of the 
samples was reported to be as follows: The 
sodium was 99.7 percent pure and contained no 
heavy metals and not over 0.04 percent calcium. 
The potassium was 99.0 percent pure with 
sodium as the other main constituent. Samples 
of both metals contained as other impurities 
traces of the oxide, hydroxide, and carbonate. 


RESULTS 


The atomic absorption coefficient, ua, is de- 
fined by the relation 


Ma=Ts—Zoy, (1) 











t ry 





o 
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where 7,4 is the photoelectric absorption coeffi- 
cient, Z the atomic number, and o, the total 
scattering per electron. As discussed presently, 
for elements of low atomic number and for the 
present wave-length region, the scattering, Zo,., 
is the greater part of the observed absorption. 
Several theoretica) formulae have been proposed 
which express the variation of ¢, with wave- 
length. Of these the Klein-Nishina formula! 
predicts values of o, which differ considerably 
from those of the earlier theories, and the Klein- 
Nishina values have been shown to agree with the 
observed scattering at wave-lengths from 20 to 50 
X.U.* The present report is largely concerned 
with the comparison of Ze, as predicted by 
theory and as experimentally observed in the 
wave-length region 40<\’<100 X.U. 

The variation of r4 with Z and \ may be 
written as 


tTa=KZ°x¢4, (2) 


where K is a constant for \<Ax jimit, and p and g 
are constants, approximately 4 and 3, respec- 
tively. [More precisely, as discussed elsewhere,’ 
p and q may be functions of Z and X. | It has been 
shown that for gamma-rays the photoelectric 
absorption in elements Z=1 to Z=16 is less 
than experimental error.'° Hence, because of 
the relation in Eq. (2), it may be expected that 
the photoelectric absorption in carbon, Z=6, 
is less than the experimental error of us at wave- 
lengths considerably greater than those of 
gamma-rays. Therefore, the absorption measure- 
ments for carbon provide a test of the validity 
of the Klein-Nishina formula. 

The experimental results for carbon are given 
in Table I. In Fig. 1 are plotted the experi- 
mental values of u4/Z and also the theoretical 
values of the total scattering per electron, o,, 
as calculated from the Klein-Nishina formula. 
The values of us/Z and o, are in excellent 
agreement from 40 X.U. to the wave-length at 
which the photoelectric absorption becomes 
appreciable, about 100 X.U., which confirms 
within experimental error the Klein-Nishina 
formula for this wave-length region. 

The absorption coefficients determined for Na, 
Al, K, Ni, and Cu are given in Table II, and for 





® Read and Lauritsen, Phys. Rev. 45, 433 (1934). 
10 Rutherford, Chadwick and Ellis, Radioactive Radia- 
tions, p. 473. 
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TABLE I. Measured values of the mass absorption coeffi- 
crent, p/p, and the electronic absorption coefficient, Me, for 
carbon; and comparison with the Klein-Nishina values of the 
Scattering per electron, ce. 





KLE! 
EX! NISHINA 
p/p PROBABLI VAL VALUI ie —« 100 
X.U.) (cm? g™! ERROR pe X10 «, X10 u 
40 0.105 0.001 3.46 3.447 +0.5 
50 0.114 0.001 3.76 3.734 0.7 
60 0.120 0.001 3.96 3.969 —0.2 
70 0.126 0.002 4.16 4.167 —0.2 
80 0.132 0.002 4.36 4.336 +0.5 
85 0.134 0.001 4.42 4.405 0.4 
90 0.137 0.001 4.52 4.483 0.8 
100 0.142 0.001 4.69 4.610 1.6 
110 0.147 0.002 4.85 4.723 2.7 
115 0.149 0.001 4.92 4.778 3.0 
120 0.151 0.002 4.98 4.824 3.2 
130 0.154 0.001 5.08 4.915 3.4 
140 0.158 0.003 5.21 4.995 4.4 
160 0.163 0.002 5.38 5.132 4.8 
184 0.171 0.002 5.64 5.280 6.9 
209 0.177 0.001 5.84 5.395 8 
we =p 4/Z =p/p*(A/ZN 1, atomic weight: Z, atomic number 
N, Avogadro's number 
2In the calculation of the Klein-Nishina values of ¢-, given in the 
ibove table, the following values of the general physical constants 
were employed 
€ 4.767 X10 e.3.u h =6.542 X10? erg: sex 
m =0.9038 X 10°" g¢ 2.998 & 10" cm /se« 
Recently Backlin and Flemberg (Nature 137, 655 (1936)) have reported 
a new value for e (4.800 x 10 which is in good agreement with that 
utilized by Birge (Phys. Rev. 48, 918 (1935 4.8036) in obtaining a 
new value of # (6.6286). Using e =4.800 x 10 und h =6.6286 x 10™ 
an evaluation of the Klein-Nishina formula for ¢, gives 3.821 and 
4.724( X10) at 50 and 100 X.U., respectively. These are each some 


2 percent greater than the values in Table I 
S are given in Table III. The mass absorption 
curves are shown in Fig. 2. The values of 74 were 
calculated from the absorption Eq. (1) in which 
were used the values of o, given in Table I. 
The variation of 74 with \ is shown in Fig. 3, 
from which g=2.83 and 2.6, for copper and 
aluminum, respectively. However the actual 
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Fic. 1. The experimental values of u4/Z for carbon and 
the theoretical values of o, as calculated from the Klein- 
Nishina formula. The agreement in the region 40=\=100 
X.U. is within experimental error, which confirms the 
Klein-Nishina formula for this wave-length region. 
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TABLE II. The observed mass absorption coefficients; and the atomic and photoelectric coefficients 


x (X.U.) Na(11 Al(13 
u/P MA TA u/p MA TA u/Pp 
50 0.113 5.02 0.17) 0.115 
55 0.117 5.20 0.17 
60 (0.119 4.52 0.16) 0.121 5.38 0.22 | 0.136 
65 
70 0.126 4.78 0.20) 0.131 5.79 0.37 0.144 
75 0.135 6.00 0.48 
80 0.134 5.08 0.31) 0.138 6.14 0.50) 0.162 
85 0.142 6.31 0.59 
90 0.140 5.30 0.37 | 0.144 6.40 0.57 | 0.180 
95 0.147 6.54 0.64 
100 0.147 5.58 0.51) 0.150 6.67 0.68 | 0.194 
110 0.154 5.84 0.65 | 0.158 7.03 0.89 | 0.202 
115 0.162 7.20 0.99 
120 0.161 6.11 0.80 | 0.168 7.47 1.20} 0.228 
125 0.175 7.78 1.46 
‘ 


130 0.168 6.38 0.97/ 0.176 83 1.44 | 0.256 
140 0.174 6.60 1.10 | 0.191 8.50 2.01 | 0.297 
160 0.186 7.06 1.41 | 0.213 9.47 2.80) 0.366 
184? | 0.208 7.89 2.09/;0.240 10.67 3.81 | 0.459 
209% | 0.231 8.78 2.84/0.279 12.37 5.36| 0.606 


‘ry =u4—Ze, Values of ¢, are given in Table I 

values of g are of little significance since the 
photoelectric absorption is, in the case of 
aluminum at 100 X.U., only 10 percent of the 
observed absorption, and consequently contains 
an unavoidably large experimental error. 

The present experimental values of u/p for 
carbon and aluminum are within one percent 
of those recently reported by Hahn" for \>140 
X.U., but the present values differ from Allen’s® 
in some cases by as much as 10 percent. It is to 
be noted, however, that Allen’s value for carbon 
at 209 X.U. is but 2 percent greater than the 
scattering predicted by the Klein-Nishina for- 
mula, at 140 X.U. is about 1/2 percent greater, 
and at 81 to 90 X.U. is more than 8 percent 
greater than the predicted scattering. In view 
of the present measurements it is unlikely that 
the photoelectric absorption varies with wave- 
length in this manner. 

In Fig. 4 the author’s results for aluminum 
are plotted with the results of other investi- 
gators.’ "> From the 7,4 curve for aluminum, 


TABLE III. The absorption coefficients determined for sulfur. 


d (X.U.) p/p wa X 10% 74 X10% 
75 0.142 7.51 0.72 
100 0.172 9.09 1.71 
120 0.200 10.57 2.85 
130 0.219 11.58 3.72 


" Hahn, Phys. Rev. 46, 149 (1934). 
2 Chao, Phys. Rev. 36, 1519 (1930). 


Ra & 


2 W &A@B radiation 


UYKENDALL 


x [U4)} 
K(19 Ni(28 Cu(29 
BA TA /Pp MA 4 Mp MA TA 
7.41 0.32 | 0.126 12.2 1.7 | 0.129 13.5 2.7 
0.148 14.3 3.5 
8.76 1.22 | 0.163 15.8 4.7 | 0.160 16.7 5.2 
0.168 16.3 4.9 
9.27 1.35 | 0.187 18.1 6.4 | 0.182 19.1 7.0 
0.195 18.9 7.0 
10.43 2.19 | 0.218 21.1 9.0 | 0.215 22.5 9.9 
11.59 3.07/}0.250 24.2 11.7}0.259 27.2 14.2 


12.49 3.73|0.297 28.7 15.8/0.306 32.1 7 
13.01 4.040.349 33.8 20.6/ 0.371 38.9 25.2 
14.68 5.51/0.420 40.6 27.1)0.434 45.5 31.5 
16.49 7.15/0.490 47.4 33.6)0.513 53.8 39.5 
19.13 9.64/0.594 57.5 43.5,0.606 63.5 49.0 
23.57 13.82 | 0.771 74.6 60.2/0.830 87.0 72 
29.56 19.53) 1.10 106.5 91.7) 1.17 122.6 107 
39.03 28.78 | 1.56 151. 136. | 1.63 170.9 15 


W Aq radiation 


Fig. 3, it appears that at 50 X.U. the value of 
ra is greater than 0.12X10-*. Assuming this 
minimum value, we have ws =4.98 X10-**, which 
is approximately 2 percent greater than the 
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Fic. 2. Mass absorption curves for the elements studied. 
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Fic. 3. The variation of the photoelectric absorption with 
wave-length. 


value reported by Read and Lauritsen.’ At the 
same wave-length their value for carbon is about 
3 percent less than that predicted by the 
Klein-Nishina formula. These discrepancies may 
be due to the fact that Read and Lauritsen 
employed a much smaller resolving power than 
was used by the author. 

Discrepancies exist between the previously 
reported results*: 4: >." for nickel and copper and 
the present measurements in the region 140<) 
<209 X.U. Three independent determinations 
of u/p for copper at 140 X.U. gave values within 
2 percent of each other, but these values are 8 
percent less than the value reported by Hahn." 
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Fic. 4. The present results for aluminum plotted with 
those of Read and Lauritsen,’ Hahn" and Chao” show the 
natural extension of the absorption measurements to the 
gamma-ray region. The form of the curve is explained as 
follows: At wave-lengths less than 75 X.U. the absorption 
is practicall¥ all scattering; at longer wave-lengths the 
photoelectric absorption becomes important, increasing to 
about 60 percent of the total absorption at 200 X.U. 


Read has measured the absorption in copper in 
the region 20 to 47 X.U.: a smooth curve through 
his results gives at 50 X.U. a value within 2 
percent of the present result. 

The probable errors given in Table I were 
calculated in the usual manner from the varia- 
tions of the individual readings obtained during 
the determination of the mass absorption coeffi- 
cient. After an interval of more than two months 
many of the values of u/p were redetermined 
showing good agreement with the earlier results. 
Considering all the sources of error, the writer 
estimates that the values of u/p given in Tables 
I, II, and III are in error by no more than 
2 percent. 
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Employing a specially constructed two-crystal spec- 
trometer and a 600-kilovolt x-ray machine, mass absorption 
coethcients have been measured, with an average probable 
error of 0.8 percent, for Pb(82), Ta(73), Sn(50), Ag(47), 
Mo(42), and Cb(41) in the wave-length range 30 <\<185 


x. Klein-Nishina 


formula for absorption due to scattering, the photo- 


Assuming the correctness of the 


electric absorption rt, was calculated for the above ele- 
ments. In the empirical equation tz=C.Z?X, p increases 


INTRODUCTION 


HEN a beam of monochromatic x-rays 
passes through an absorbing material, the 
law 


I = [pe #al N/A) (02) | (1) 


is obeyed, where Jo is the intensity of the incident 
beam, J the intensity after passing through a 
thickness x of absorber of density p, ua. the total 
atomic absorption coefficient, N Avagadro’s 
number, and A the atomic weight. 

Theoretical and empirical laws have been 
developed expressing the variation of yu, with Z, 
the atomic number, and \ the wave-length of 
the x-rays; the empirical laws take the form 


Ma = Ta(Z, 4) +o(Z, d), with ra(Z, A) = C.Z?X4, (2) 


where 7,(Z, \) is the true or photoelectric atomic 
absorption coefficient, ¢,(Z,\) is the atomic 
absorption coefficient due to scattering, C, is a 
constant over a given wave-length range deter- 
mined by the absorption limits. 

In the wave-length range 0.1<’<0.7A, Richt- 
myer! determined p=4, g=3, 0.1<¢/p<1.0, and 
C,=22.4X10-*7 for \ less than the K limit and 
6<Z<82; ¢/p=o,(N/A). Allen? made a large 
number of measurements in the range 0.081<X 
<2.0A and as a result advocates a value of 
q=2.92 although some of his values came out as 
low as 2.6. Jaeger? measured p/p, the mass 

1F, K. Richtmyer, Phys. Rev. 17, 264 
(1921); 27, 1 (1926); 30, 755 (1927). 

2S. J. M. Allen, Phys. Rev. 24, 1 (1924); 27, 266 (1926); 


28, 907 (1927). 
3R. Jaeger, Zeits. f. Physik 69, 565 (1931). 


1921); 18, 13 


from 3.51 140 to 50 X.U. 
and q increases from 2.60 to 2.80 as Z decreases from 82 


to 3.88 as \ decreases from 


to 41. Agreements with data of other observers for Pb(82 
are given. At 30 and 40 X.U., the present data for Pb(82 
are about 10 percent higher than the theory developed by 
Hulme, McDougall, Buckingham and Fowler. New con 
stants in Gray’s empirical equation for rpp are suggested of : 
a=3.7321, b=1.03 and c=0.44. 


absorption coefficient, for Pb(82) in the region 
30<A<127 X.U., 
method. Later Alichanjan and Kosman‘* made 
similar measurements in the range 25<\<60 
X.U., introducing a correction for the effec- 
tive wave-length when using the end-radiation 
method. Hahn® made very careful absorption 
measurements on Pb(82), Au(79), W(74), Ta(73), 
Ag(47), Cu(29), Al(13), C(6) and paraffin in the 
range 139<\<208.6 X.U. The data at 208.6 
X.U. were taken with a two-crystal spectrometer ; 
the rest of the data with the second crystal 
removed. Recently, Read*® measured absorption 
coefficients for Pb(82) and Sn(50) in the range 
20.2<A<52.5 X.U. rocksalt 
crystal; the width of the band of rays used by 
him about 5 X.U. Read estimates his 
experimental error on the lead data to be 2 


employing the end-radiation 


using a_ single 


was 


percent. 

In the region 4.7<\<47 X.U., Chao’ meas- 
ured absorption coefficients in Pb(82), Sn(50), 
and Al(13) employing gamma-rays. Gray*® sum- 
marized gamma-ray data for Pb(82) and de- 
termined the constants in an empirical equation 
suggested by Ellis, to fit this summary, using 
Allen’s data at 100 X.U.; this equation for r, 
the linear photoelectric absorption coefficient, is 


logio 7=a+logyo A+c(logy A)?, (3) 


* A. Alichanjan and M. Kosman, Zeits. f. Physik 90, 779 
(1934). 

> T. M. Hahn, Phys. Rev. 46, 149 (1934). 

6]. Read, Proc. Roy. Soc. A152, 402 (1935). 

7C. Y. Chao, Proc. Nat. Acad. Sci. 16, 431 (1930); Phys 
Rev. 36, 1519 (1930). 

8 LL. H. Gray, Proc. Camb. Phil. Soc. 27, 103 (1931). 
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where a=3.6505, b=1.0 and c=0.480. This law 
fits gamma-ray data within about 10 percent. 

A large number of theoretical calculations have 
been made for 7(Z, A); the best of them, calcu- 
lated by wave mechanics, gives values for 7 
lower than experiment in the x-ray region. The 
most recent calculation by Hulme, McDougall, 
Buckingham and Fowler® gives values agreeing 
with Gray’s Eq. (3) above. Their formula is 
valid for all wave-lengths and atomic numbers 
but is impractical for all but very short wave- 
lengths. 

Of the theoretical calculations for ¢,(Z, \), the 
absorption due to scattering, the Klein-Nishina 
formula has been checked experimentally by 
Chao,’ by Read and Lauritsen,'® and by Cuyken- 
dall;" it represents the results within experi- 
mental error. These checks were carried out on 
the elements of low atomic number, C(6) and 
Al(13), and at short wave-lengths (A< 100 X.U.), 
because the photoelectric absorption may be 
considered negligible in this range. 

Since for gamma-rays (A<50 X.U.) the ab- 
sorption due to scattering plays the major role 
and for x-rays (A>100 X.U.), the photoelectric 
absorption is the larger part of the total, at 
least for the heavier elements, the intermediate 
region 30<A<150 X.U. for x-rays is an im- 
portant one for the investigation of the variation 
of « and +r with the wave-length and atomic 
number. In order to study the photoelectric 
absorption in this region, it is better to use the 
elements of high atomic number for which ¢ is 
small compared to +. There are almost no 
reliable data in this region for monochromatic 
rays. End-radiation measurements are not very 
accurate and the radiation employed is not 
monochromatic. Any measurements in this region 
with single crystal spectrometers as mono- 
chromators were taken with very low resolving 
power of the spectrometers. Therefore a study 
was undertaken of the photoelectric absorption 
in the range 30<A’A<200 X.U. for elements of 
high atomic number. 


°H. R. Hulme, J. McDougall, R. A. Buckingham, and 
R. H. Fowler, Proc. Roy. Soc. A149, 131 (1935). 

1 J. Read and C. Lauritsen, Phys. Rev. 45, 433 (1934). 
1 T. R. Cuykendall, Phys. Rev. this issue. 
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APPARATUS 


The apparatus used consists of a 600-kilovolt 
x-ray machine and a two-crystal spectrometer 
described by Cuykendall and Jones.” The source 
of the x-rays is a cascade-type Coolidge tube, 
with a tungsten target, which is operated on two 
300-kilovolt induction coils in series. The mono- 
chromator is a two-crystal spectrometer specially 
designed for use in the wave-length region 
30<A<215 X.U. Transmission of the x-rays 
through the body of the crystals and reflec- 
tion from internal planes instead of the usual 
reflection from the surface, is employed. An 
ionization chamber containing argon at 80 
atmospheres pressure, and a F P-54 electrometer- 
tube system is used to measure the intensity of 
the reflected radiation. The width of the wave- 
length band separated from the continuous 
spectrum produced by the x-ray tube, is of the 
order of 1.8<6\<1.0 X.U., for 50<\<215 
X.U." 


RESULTS 


The elements investigated were Cb(41), 
Mo(42), Ag(47), Sn(50), Ta(73) and Pb(82) in 
the range 30<\<185 X.U. All were of high 
purity and known to be free of any impurities 
which would affect the results beyond experi- 
mental error. The silver samples were 99.98 
percent pure containing no lead; most of the 
impurity was oxygen. Table I gives a summary 
of the observed mass absorption coefficients y/p, 
atomic absorption coefficients y., and photo- 
electric atomic absorption coefficients 7.. Each 
value is obtained from 9 to 12 pairs of observa- 
tions with an average probable error of 0.8 
percent. 

The observed values of the atomic absorption 
coefficients for Pb(82) are plotted in Fig. 1 
along with the data of other observers. Richt- 
myer’s values at wave-lengths less than the K 
limit agree with the observations. Richtmyer'’s 
and Hahn’s data at 208.6 X.U. agree with the 
observed value within experimental error; this 
was taken as a check. Allen’s values, in the 
range 50<\<110 X.U., are from 40 to 60 


2 T. R. Cuykendall and M. T. Jones, Rev. Sci. Inst. 6, 
356 (1935). 
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TABLE I. Observed absorption coefficients. 
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Pb(82 Ta(73 Sn(5 
p/p tg X10 Tq X10 y/o la X 10% x 102 Pa 
0.22 0.752 0.499 
0.394 1.346 1.064 
0.600 2.052 1.745 0.458 1.372 1.100 0.482 
0.900 3.077 2.751 0.689 2.064 1.774 0.615 
1.455 4.975 4.634 1.018 3.049 2.745 0.842 
1.829 6.254 5.896 1.357 4.065 3.747 1.189 
2.490 8.513 8.145 1.819 5.448 5.120 1.567 
3.248 11.11 10.73 2.435 7.293 6.956 1.998 
4.079 13.94 13.56 3.014 9.027 8.681 2.394 
5.150 17.61 17.21 3.770 11.29 10.94 3.108 
6.412 21.93 21.52 4.553 13.64 13.28 3.745 
5.476 16.40 16.03 4.524 
5.108 
\g(47 Mo(42) Cb(41 
ale lig X 102 x 10 n/p a X 102 Tq X10 Hy x 
0.229 0.407 0.232 0.177 0.280 0.124 0.278 
0.283 0.504 0.317 0.227 0.360 0.193 0.356 
0.400 0.712 0.516 0.319 0.505 0.330 0.469 
0.544 0.969 0.764 0.412 0.653 0.470 0.617 
0.715 1.273 1.062 0.535 0.847 0.658 0.814 
0.900 1.602 1.386 0.650 1.029 0.836 1.025 
1.081 1.925 1.702 0.821 1.300 1.102 1.237 
1.362 2.425 2.198 1.036 1.643 1.440 1.636 
1.656 2.948 2.717 1.243 1.969 1.762 1.982 
2.038 3.628 3.392 1.513 2.397 2.185 2.389 
2.108 3.339 3.123 3.420 
3.193 5.057 4.837 4.750 
| T 
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Fic. 1. Comparison of present values of ua for Pb(82) with those of other 


observers. 
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percent higher than the present data. This 
disc:epancy could be accounted for by assuming 
that the radiation he employed was not mono- 
chromatic due to low resolving power of his 
apparatus in this region. Jaeger’s data agree with 
Allen’s but, as pointed out by Alichanjan and 
Kosman,‘* a correction should be made to the 
wave-length of the short-wave limit of their 
operating voltage, when employing the end- 
radiation method, to take into account the fact 
that a homogeneous beam of x-rays is never 
really reached by this method; the effective 
wave-length is longer than the short-wave limit. 
This would move Jaeger’s data more in agree- 
ment with the present observations. Alichanjan 
and Kosman’s data, for which this correction 
has been made, agree with the data herein 
reported. Read’s data agree with the observed 
data within experimental error on the wave- 
length region covered by him. The gamma-ray 
measurements of Chao agree within 15 percent 
where they overlap with the present data. 

Hahn’s data for Ta(73) are lower than an 
extrapolation of the present data in the region 
140<\<160 X.U. beyond the experimental 
error. His data for Ag(47) are higher in the same 
range. The data of Allen are consistently 10 to 
50 percent higher than the present data. Read's 
data for Sn(50) fit on to an extrapolation of the 
observed data for \<50 X.U. 

The Klein-Nishina formula for the absorption 
due to scattering has been checked within 
experimental error;’: '°." therefore, if one as- 
sumes that the scattering per electron o, is 
independent of the atomic number, the photo- 
electric absorption r, per atom may be calculated 
from ta=ua—Zeo,.* Plotting these values of 1, 
against wave-length on a double logarithmic 
scale should give a straight line if g in Eq. (2) is 
constant; the slope of such a straight line gives 
the value of g. Fig. 2 is such a graph for the 


TABLE II. Experimental values of q in Eq. (2). 


42 47 50 73 82 


41 
2.79 80 2.63 2.63 2.60 2.60 


nN 


<a 


* As an extreme case, if, due to coherent scattering, a, 
were 5 percent greater for Pb(82) than calculated by the 
Klein-Nishina formula, at \=30 X.U., ra would be less 
than 2 percent higher which is of the order of the experi- 
mental error and would not seriously affect the conclusions 
drawn herein. 
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F iG. 2. tr. plotted against wave-length for various elements 


elements studied. The data for each element fall 
nearly on straight lines. Table II gives the values 
of q found from these curves. 

In a similar manner, the exponent p, of Z in 
Eq. (2), may be determined from a double 
logarithmic plot of +r, against Z. Fig. 3 is the 
graph obtained. Again the data fall nearly on 
straight lines. Table III gives the values of p. 
The variation of p agrees qualitatively with the 
theory" of photoelectric absorption which predict 
that p=5 as \=0 and should decrease as the 


? 


TABLE III. Experimental values of p in Eq. (2). 


4 (X.U.) p A (X.1 p 
140 3.51 90 3.69 
130 3.64 80 3.73 
120 3.62 70 
110 3.68 60 3.88 
100 3.69 50 3.84 


8F. Sauter, Ann. d. Physik 9, 217 (1931); 11, 454 
(1931). 
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FG. 3. tr. plotted against atomic number for various wave- 
lengths. 


wave-length increases. In the wave-length region 
here studied these data indicate that if there is 
a variation of g with X, it is much less than that 


predicted by theory. 

The agreement with Gray’s empirical Eq. (3) 
is not good for the values of the constants given 
by him. Gray used Allen’s data at 100 X.U. in 


r. JONES 
TABLE IV. Values of r,X 10” for Pb(82) 
CALCULATED FROM Eo. (3) 
using using CALC 
a =3.6505 a =3.7321, FROM 
\ (X.1 OBSERVED b=1.0, 0.48 b =1.03 0.44 THEORY 
30 0.50 0.45 0.49 0.47 
40 1.06 0.92 0.98 0.91 
50 1.74 1.64 1.70 
60 2.75 2.62 2.69 
70 4.63 4.06 4.07 
SO 5.90 5.90 5.82 
90 8.14 8.26 8.02 
100 10.73 11.20 10.74 
110 13.56 15.51 14.63 
120 17.21 19.26 17.98 
130 21.52 24.45 22.63 


the determination of his constants which partly 
explains the lack of agreement. Suggested values 
of the constants of Eq. (3) to fit the present 
data are: a=3.7321, b=1.03 and c=0.44. Table 
IV gives the observed values, values calculated 
from Gray’s equation, and values calculated 


using these new constants. Of course it is 
obvious that an equation in which g varies with 
wave-length can fit only approximately the 


present data for which g has been shown to be 
constant; therefore such an equation has no 
theoretical significance, but may be of practical 
use in irradiation work with hard x-rays. 

Calculations taken from the theoretical photo- 
electric absorption curve given by Hulme, 
McDougall, Buckingham and Fowler? for Pb(82) 
are given in the last column of Table IV. The 
calculated values are about 10 percent lower than 
the observations. 

The author wishes to express his appreciation 
to Professor F. K. Richtmyer for his assistance 
and encouragement throughout the work and to 
the Raritan Copper Works, Perth Amboy, N. J., 


for donating the silver samples. 























PHYSICATI 





VOLUMI 50 


REVIEW 


Mechanical Detection and Measurement of the Angular Momentum of Light 
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The electromagnetic theory of the torque exerted by a 
beam of polarized light on a doubly refracting plate which 
alters its state of polarization is summarized. The same 
quantitative result is obtained by assigning an angular 
momentum of fh h) to each quantum of left (right) 
circularly polarized light in a vacuum, and assuming the 
conservation of angular momentum holds at the face of the 


plate. The apparatus used to detect and measure this effect 


ELECTROMAGNETIC FIELD THEORY 


HE moment of force or torque exerted on a 

doubly refracting medium by a light wave 
passing through it arises from the fact that the 
dielectric constant K is a tensor. Consequently 
the electric intensity E is, in general, not parallel 
to the electric polarization P or to the electric 
displacement 


D=KE=E+4rP 


in the medium. The torque per unit volume 
produced by the action of the electric field on 
the polarization of the medium is 


1=PxXE=(DxXE)/4r. (1) 


Following Sadowsky and Epstein' we may 
calculate this torque for a simple case as follows: 
Assume a doubly refracting medium of perme- 
ability unity and take as x, y, and s axes the 
principal axes of the tensor K for the light fre- 
quency in question. Denote the principal values 
of K by 2, n,2 and n/ so that n,, n, and n, will 
be the principal indices of refraction. Let the 


electric components of a plane light wave 

*Member of the Physics Department, Worcester 
Polytechnic Institute. Preliminary work and final calcu- 
lations were carried on at Worcester. Experimental work 
was done on leave-of-absence from Worcester as Research 
Associate at Palmer Physical Laboratory, Princeton 
University. Additional measurements were made by Mr. 
W. Harris at Princeton, as will be explained later. 

1A, Sadowsky, Acta et Commentationes Imp. Uni- 
versitatis Jurievensis 7, No. 1-3 (1899); 8, No. 1-2 (1900). 
P. S. Epstein, Ann. d. Physik 44, 593 (1914). The deriva- 
tion is repeated here because the first articles are in 
Russian and relatively inaccessible, while in the last there 
seems to be a misprint in the result. I wish to thank 
Dr. Boris Podolsky for translating parts of Sadowsky’s 
articles for me from the Russian. I am also deeply indebted 
to Professor A. Einstein not only for his advice in checking 
Professor Epstein’s calculation but also for several inter- 
esting discussions about the experimental part of the 
work here reported. 


was designed to enhance the moment of force to be 


measured by an appropriate arrangement of quartz wave 


plates, and to reduce interferences. The results of about 120 
determinations by two observers working independently 
show the magnitude and sign of the effect to be correct, 
and show that it varies as predicted by the theory with 
each of three experimental variables which could be 


independently adjusted. 


propagated in the +s direction be given by 


E=A cos @ cos (Z,+A), 
A sin @ cos (Z, —A), QO, 


D=KE=n.,7A cos @ cos (Zj+A), 
n,2A sin @cos(Z,;—A), 0, (2) 


where Z,;=w(t—nz/c); n=(n,+n,)/2, 


A=n2(m,—mn,)/A; w=2c/X. 


For any value of s (in particular s=0) for which 
the wave plate thickness A/z is a whole number 
the 
angle +@ with the x axis. 

The torque per unit volume at any point in 


light is plane polarized and E makes an 


the crystal is calculated from (2) according to (1) 
1=0, 0, 


(A?/8x)(n/—n,?) sin 26 cos (Z;+A) cos (Z;—A) 
and the time average value of the s component is 
(A /42)?(n,2—1n,7) sin 20 cos 2A. 


Integrating from 2; to s2 we get the torque per 
unit area on the crystal between these values: 


9 


= —(A/47)*nX sin 20(sin 2As—sin 2A;). (3)? 


This torque, which has hitherto been generally 
considered too small for experimental detection, 
was detected and measured in the present ex- 
periment.® 


* J. H. Poynting, Proc. Roy. Soc. A82, 560 (1909), infers 
from an ingenious mechanical analogy that circularly 
polarized light should exert a torque equal to the light 
energy per unit volume times A/2m on unit area of a 
quarter-wave plate which makes the light plane polarized. 
This result, which neglects surface reflections, is contained 
in (3). 

>R. A. Beth, Phys. Rev. 48, 471 (1935); also Annual 
Meeting of the Am. Phys. Soc., St. Louis, Missouri, 
January 1, 1936. A. H.S. Holbourn, Nature 137, 31 (1936), 
also reports successful measurement of the effect. 
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QUANTUM THEORY 


It is worth observing that (3) may be derived 
from the quantum theory by assigning an angular 
momentum or spin of i(—h)* to each photon of 
left (right) circularly polarized light, the spin 
axis being in the direction of propagation of the 
light.° 

Consider an elliptically polarized light wave, 
propagated in the +2 direction in a vacuum, 
whose components in the directions of the prin- 
cipal axes of the ellipse (phase difference = 7/2) 
have amplitudes Xo and Yo. Let the components 
of the same when resolved along an 
arbitrary pair of perpendicular x and y axes in 
the plane of the ellipse, be given by 


Ycos(Z—A), 0, 


wave, 


E=VNX cos (Z+A), (4) 


where Z=w(t—3/c) ad 2A is the phase angle by 
which the y component lags behind the x com- 
ponent of the wave. It may be shown that® 


XYsin2A=Xo¥o, X*°+Y°=X7+ Ye". (5) 


If we now regard the wave (4) as the super- 
position of left and right circularly polarized 
components of amplitudes L and R, respectively, 
calculation shows that 


L?-R?=XYsin2A, L?+R?=(X?+ Y?)/2. 


The number of left circularly polarized photons 
transmitted per unit area per second is the 
Poynting energy flow cL?/42 divided by the 
energy per photon 2rhv=2rhc/d, or AL?/87°h. 
Multiplying by # gives \L?/8x*, the angular 
momentum transmitted per unit area per second 
by the left circularly polarized component. 
Taking account of the similar expression for the 
right circularly polarized component, the angular 
momentum transmitted per unit area per second 
by the wave (4) is 


M =)X(L?2— R?) /8r?=XAXY sin 24/87? 
=AXnVo/8r*. (6) 


It is well known that the linear momentum 


{h=h/2r. 

5A. E. Ruark and H. C. Urey, Proc. Nat. Acad. Sci. 
13, 763 (1927); Harnwell and Livingood, Experimental 
Atomic Physics (1933), p. 81; Brillouin, Les Statistiques 
Quantiques, Chap. III (1930). 

6 See, e.g., Schuster and Nicholson, Theory of Optics, 


third edition, pp. 14, 15 (1924), especially Eqs. (15) and 
(17). 
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transmitted per unit area per second is the 
energy flow divided by c, in this case: 


(L?+ R?) /4a = (X?°+ VY?) /8a = (Xe?+ Yo?) /8x. 


Thus the angular and linear momenta are pro- 
portional to the natural invariants (5) of the 
wave (4). For X»9=Y, their ratio is X=X/2z, 
which is another form of Poynting’s result.* 

The value (6) of the angular momentum, which 
is independent of the quantum constant /, can 
be derived from the linear momentum (Poynting 
vector divided by c) in a finite beam, but the 
elementary method given is more descriptive in 
bringing out the fact that the quantum theory 
leads to the same value (3) as the wave theory 
for the torque on a section of a crystal. 

To show this equivalence we have now to 
obtain the angular momentum for a wave (2) 
in the crystal from the expression (6) for a wave 
in a vacuum, by placing the face of the crystal 
at z=2, (vacuum for z<z;) and using the con- 
servation of angular momentum and Fresnel’s 
expressions for the amplitudes of the reflected 
and transmitted waves. In the reflected wave 
the y component of amplitude Y(m,—1)/(m,+1) 
lags behind the x component of amplitude 
X(nz—1)/(nz+1) by the phase angle 24. We 
get the angular momentum transmitted through 
the face z=2, by calculating (6) for the direct 
and reflected waves in the vacuum (taking 
account of direction of transmission of each) and 
assuming the conservation of angular momentum 
at the face 

M,=X Y(n,+n,)d sin 24/42°(n,+1)(m,+1). 
Using Fresnel’s expressions for the amplitudes, 
we identify the transmitted wave in the crystal 
with (2) at by choosing 


A cos 6=2X/(n.,+1); Asin 0=2Y/(n,+1); 


=e 
2=2) 


A=A)= 72;(",—m,)/X. 


Hence the angular momentum transmitted by 
(2) at s=2; is 


M,=(A/47)?nX sin 26 sin 2A}. (7) 


It can now be seen directly that the torque (3) 
is equal to the excess of the angular momentum 
(7) per unit area per second flowing into a 
section of the crystal at s=2,; over that flowing 
Hence both the field and the 


=? 
2=2Z. 


out at 
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Fic. 1. Diagram of apparatus. 


quantum theories predict the same value for 
the effect which has been measured in the present 
experiment. 


APPARATUS 


The basic idea used in detecting and measuring 
the effect’ is to observe the deflection of a quartz 
wave plate hung from a fine quartz fiber when 
suitably polarized light is sent through the plate. 
In order to increase the effect to be measured and 
to avoid as far as possible interfering effects due 
to mechanical and electrostatic disturbances, 


7See Poynting, reference 2, A. Kastler, Société des 
Sciences physiques et naturelles de Bordeaux, Jan. 28, 
1932. R. A. Beth, abstract, Boston Meeting, American 
Physical Society, Phys. Rev. 45, 296 (1934). 
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Fic. 2. Photograph of apparatus. 


stray heating effects, radiometer and gas effects, 
and variations in light pressure, the apparatus 
was constructed as shown in Figs. 1, 2, 3, and 4. 
Some of the precautions taken to avoid these 
interferences may have been unnecessary, but 
they were sufficient. 

The whole apparatus is supported on the 
heavy cast iron brackets I in Fig. 1, each 
12”"X12”, which are bolted to a brick pier to 
avoid mechanical disturbances, especially tor- 
sional vibrations of a period comparable to that 





PLANE OF ELECTRIC VECTOR 
OF LIGHT FROM NICOL 


Fic. 3. Wave plate arrangement. 
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Fic. 4. Detail showing method by which light beam could 
be given lateral adjustment. 


of the torsional pendulum formed by the one- 
inch circular wave plate ./ hanging at the bottom 
of the quartz fiber (about ten minutes). The fiber 
is about 25 centimeters long and hangs in a 
quarter inch hole bored lengthwise through a 
copper cylinder two inches in diameter. The 
upper end of the fiber is attached with a bit of 
flake shellac to a conical copper plug which was 
carefully ground in place and sealed with 
Apiezon grease. Thus the zero position as well 
as the amplitude of the torsional pendulum can 
be changed while the chamber is evacuated. 

The hanging plate and the fixed plate above 
it, Fig. 3, are located in the cylindrical vacuum 
chamber C which is of copper, three inches in 
diameter with walls one-half inch thick. The 
heavy copper chamber is intended to eliminate 
unequal heating effects from the room as well as 
to shield off stray light and electrostatic dis- 
turbances.§ 

The chamber was evacuated through wide 
tubing by a high speed diffusion pump and three 
liquid-air traps. A two-stage diffusion pump and 
forepump provided the forevacuum for the high 
speed pump. According to the rated air calibra- 
tion of the Western Electric ionization gauge, 
pressures below 10~-® mm of mercury existed in 
the chamber while the effect was being measured; 
no noticeable fluctuations appeared in the actual 


*] wish to thank Professor G. P. Harnwell for this and 
many other helpful suggestions concerning the mechanical 
and optical design of the apparatus. 
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measurements which might be ascribed to radi- 
ometer or gas effects. 

Light from a three-millimeter tungsten ribbon 
filament F is focused by the fused quartz lens L, 
10.7 centimeters in diameter, through a large 
Nicol prism N, the lower plate B (Fig. 3), the 
fused quartz window W, and the hanging plate 
M, on the reflecting layer of aluminum on the 
top of the upper plate 7 in the chamber. The 
maximum deviation of the light from the vertical 
was about 10°. W has a circular aperture one inch 
in diameter. A round one-quarter inch copper 
disk d fastened in the middle of W prevents light 
from passing through the hole // in the top plate 
T. With this arrangement no light energy 
reaches the fiber and most of the energy is 
reflected out of the vacuum chamber altogether, 
thus minimizing fluctuations in the position of 
the pendulum due to unequal heating of the 
fiber and light pressure on the small mirror m. 
These undesirable effects are found to be present 
if the shield d is removed. 

The lower plate B is placed outside the vacuum 
chamber in a brass ring bushing held so it is free 
to turn in a brass frame attached to the upper 
cast iron bracket. The position of the lower wave 
plate may be read on a circular scale by means 
of a pointer attached to the bushing. A circular 
brass plate P attached to the under side of the 
bushing carries adjustable stops S, and the plate 
may be rotated from one stop to the other while 
observing the swing of the pendulum (with a 
telescope placed in front of the apparatus) by 
means of cords running in grooves on the edge 
of the plate P. Frame, plate P, cords, stops, 
pointer and circular scale are visible at the top 
of Fig. 4. 

The ideas underlying the wave plate arrange- 
ment, Fig. 3, may be qualitatively understood 
as follows. The light coming upward through the 
plates is reflected by the aluminum player on the 
top side of 7, passing downward again through 
the plates. For a certain wave-length A», B and 
T are quarter wave plates with their axes as 
shown at 90° to those of the hanging half-wave 
plate M. If light of wave-length >» from the 
Nicol, plane polarized at @=45°, enters the 
bottom plate B, both direct and reflected beams 
will be circularly polarized in each space between 
the plates in the directions indicated by the 














ANGULAR 


curved arrows. Thus the angular momentum 
delivered per second to the half-wave plate WM 
is almost four times what could be obtained by 
simply absorbing the same amount of circularly 
polarized light. 

For wave-lengths \ different from Xo it will 
be seen that, because of the alternate advancing 
and retarding of each component with respect 
to the other, light in both direct and reflected 
beams is plane polarized at the initial angle @ at 
the bottom of B, the middle of M, and the top 
of T. If the plates vary somewhat from the exact 
specifications these levels of plane polarization 
will be correspondingly shifted. In general the 
light will be elliptically polarized in each space 
between the plates. Careful consideration shows 
that the angular momentum delivered to is in 
the same sense as that for Xo for all wave-lengths 
\ for which 


(My—Mz),/A<2(My—Nz)y,/do 
or, approximately, \>Apo/2. 


If the lower plate B is rotated 90° from the 
position shown, the torque due to light of wave- 
length Xo will be just reversed. For this position 
the contributions to the torque of other wave- 
lengths is not easily considered qualitatively, but 
detailed calculation shows, as might be expected, 
that the integrated torque for all values of \ 
used is smaller in magnitude and opposite in 
direction compared to the torque similarly cal- 
culated for the position shown in Fig. 3. The 
torque is a continuous function of the position 
angle of the bottom plate, and its value can, with 
some labor, be calculated for each position. 

The design wave-length \» is determined on 
the basis of (3) by noting that A? is practically 
proportional to the energy in the light at each 
wave-length, and that \ enters as a factor. The 
energy distribution has a rather sharp peak in 
the region just above 1.0 yw for the tungsten 
filament temperatures to be considered. By a 
series of trial calculations \y>=1.2 » was chosen 
to make the light torque a maximum. 

The actual plates used, of course, vary some- 
what from the ideal specifications. The most 
accurate plates obtained were made to within a 
few percent of the specified thickness by the 
Bryden Company of Waltham, Massachusetts, 
whom I wish to thank for this excellent coopera- 
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OF LIGHT 


tion. For the purposes of the exact calculations 
the plates were measured by means of a Babinet 
type compensator with sodium light. From the 
known values of the refractive indices of quartz® 
the retardation calculated for each 
plate for all wave-lengths contributing to the 
effect. 

Besides almost quadrupling the torque com- 
pared to what would be obtained if the light 


was then 


were simply absorbed, and furthermore reducing 


heating, radiometer, and effects in the 


vacuum chamber, this plate arrangement greatly 


gas 
reduced difficulties due to radiation 
The disk M cannot be mounted exactly at right 
angles to the fiber axis, nor can the light be 


pressure. 


projected exactly in a vertical direction. Thus 
the resultant light pressure integrated over the 
disk will not be quite vertical and will produce a 
“light pressure torque” unless its line of action 
lies in a vertical plane containing the fiber axis. 
Without special precautions this torque might 
very easily mask the effect to be detected and 
measured.!° 

A tungsten filament light source was chosen, 
at a sacrifice of light intensity compared to an 
arc, in order to keep the resultant light pressure 
on the disk M as steady as possible, both in 
position and magnitude, and to give a steady 
and reproducible spectral energy distribution for 
purposes of calculation. The undesirable result- 
ant light pressure on the disk M is decreased in 
the plate arrangement described, first, because 
most of the energy is allowed to pass through the 
disk and, secondly, because the pressure due to 
surface reflections and absorption for the upward 
beam is largely cancelled by the similar pressure 
for the beam reflected downward. 

Finally, the moment arm for the residual light 
pressure can be varied by shifting the entire 
light beam parallel to itself. This is done by 
mounting the entire optical system up to A —A 
in Fig. 1 on the double slide D, which in turn is 
supported in the lower of the two cast iron 
brackets J. The detail view, Fig. 4, shows the 
two micrometer screws by which the lateral 


9 International Critical Tables, V1, 341-342. 

10 Besides the smallness of the effect itself, the light 
pressure torque is perhaps the greatest difficulty in this 
experiment. I wish to thank Professor A. Kastler of the 
University of Bordeaux for a personal communication 
emphasizing this fact. 
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motion in two coordinates is produced. These 
are calibrated in thousandths of an inch so that 
any particular setting is easily repeated at a 
later time. By using this arrangement and leaving 
the lower wave plate B in a fixed position, the 
change in the sum of light torque, and light 
pressure, radiometer, and effect torque 
produced by a ten-percent change in filament 


gas 


current was measured for various settings of 
one of the double slide screws, the other being 
kept fixed. The results when plotted give a 
smooth curve crossing the zero torque line when 
the filament image lies across the middle of the 
window W as nearly as the eye can judge. This 
indicates that the spurious torques are probably 
of the nature anticipated and, in any event, can 
be held steady and made very small by appro- 
priate setting of the double slide. In making the 
measurements of the light torque effect itself, 
the filament current was held constant to better 
than one percent, and the filament image was 
placed across the middle of the window. Care 
was also taken to use the same lengthwise 
portion of the filament in all cases. 

The entire optical system can be rotated in its 
support on the double slide. Changes in the 
angle @ so produced are read on the circular scale 
visible in Fig. 4 on top of the double slide and 
below the Nicol housing NV. 

As described below, changes in the light torque 
effect corresponding to definite rotations of the 
bottom plate were measured and found to be 
reproducible for the same settings to within a 
few percent. When the bottom plate was re- 
moved, or when it was replaced by a fused quartz 
plate, rotation in the same manner produced no 
measurable effect. Hence it may be assumed that 
rotating the lower plate leaves the light energy, 
light heating, and gas 
effects sensibly unchanged. In other words, the 
apparatus used successfully eliminates the inter- 


pressure, radiometer 


fering effects, and enhances the light torque 


effect to be measured. 


MEASUREMENTS 


The apparatus constructed according to the 
plan described showed the effect in the predicted 
direction and of about the predicted magnitude 
at the first trial on July 10, 1935. With this, the 





A. BETH 
first aim, that of detecting the effect, was accom- 
plished. The second aim, that of measuring the 
effect, was carried out by a resonance method 
(see reference 7). One of the first series of meas- 
urements taken is shown in Fig. 5, which illus- 
trates both the definiteness of the effect observed 
and the kind of minor fluctuations encountered. 
At the top of Fig. 5 are given the readings 
a,(v=0, 1, 2, -) in centimeters on the scale 
(338 centimeters from the fiber axis) correspond- 
ing to the end points of the swing of the torsional 
pendulum. Tenths of millimeters on the scale 
could be estimated by means of the cross hairs 
in the telescope. The bottom plate B was in the 
position shown in Fig. 3 (@=45°, right stop R) 
while the readings ado to a, were being taken, but 
was turned 90° to the left (to @= —45°, to left 
stop L) at the instant at which the end point a; 
was read. The change in light torque produced 
by this rotation of B was such as to oppose the 
swing of the pendulum from a; to as. At ag the 
plate B was turned back to R, the torque change 
again opposing the ensuing swing from d¢ to a. 
The plate was turned similarly in ‘‘antireson- 
ance” at ad; and ag as indicated, but was then 
left in the position R when dy, di) and a, were 
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Effect of light torque in resonance and out of 


Fic. 5. 
the torsion 


resonance on the vibration amplitude of 
pendulum. 
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read. At dy. the plate was turned R to L, the 
change in torque now aiding the swing of the 
pendulum between dy. and a3. At a;3 the plate 
was turned L to R. This turning 
was continued at a,, and a@,;, and then the plate 
was again left in the position R for the drift 


“resonance” 


readings dig to doo. For most of the readings 
taken another ‘‘antiresonance”’ series was added 
to the schedule described, and, for the readings 
taken by Mr. Harris, the “drift” periods were 
shortened while the ‘‘resonance”’ 
middle of the schedule was lengthened to six 
half-cycles. 

From the end-point readings taken according 
—d,1+2a,—a,4, and 
calculated. The ab- 
A, are plotted 


series in the 


to such a schedule 4A,= 
4Z,=a,-1+2a,+4,41 were 
solute values of the ‘‘amplitudes” 
in the upper part of Fig. 5 and the “zeros” Z, 
in the lower part. For small values of 
slopes) the arithmetic 


(circles) 
damping (see “‘drift’’ 
means used differ by much less than the experi- 
mental error from the ideally correct geometric 
means for the quantities calculated. 

From each such set of readings we wish to 
determine the value of the change in light torque 
on the hanging plate 1 when the bottom plate 
B is turned R to L or L to R. Since the torsional 
amplitude of the pendulum was at most a few 
degrees, we may assume that the light torque 
on M was a constant for a given position of the 
bottom plate B. We may consider a change in 
the position of the bottom plate to be equivalent 
to a change in the equilibrium position about 
which the torsional pendulum executes its oscil- 
lations. Let x be the corresponding change in 
scale reading. Careful consideration shows that, 
leaving damping out of account, the change in 
A, per cycle must then be 2x during “‘resonance”’ 
or “‘antiresonance’”’ sequences. With a small 
amount of damping, x is therefore very nearly 
equal to the average of the rates of change of A, 
per half-cycle during 
" For example, in Fig. 


Ragan and “‘anti- 
resonance.’ , the magni- 
tudes of the slopes A and B, by he squares, 
are 0.341 and 0.289 centimeter per half-cycle, 
respectively. Hence x =0.315 centimeter. 

This method of taking readings and calculating 
affords several checks on the self-consistency of 
the sequence of readings a, used to determine 


one value of x. For example, the damping slopes 
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NTUM OF 


a, b, and 
(A—B)/2.in magnitude. 
during the “resonance” and 


should be approximately equal to 
Furthermore Z,+x/2 
“antiresonance”’ 
intervals (plotted as crosses) should form a 
smooth sequence with the values of Z, during the 
drift intervals. Slight disturbances, such as 
between dy and @;9, may be taken into account in 
evaluating these diagrams. 

The quartz fibers were made in a gas-oxygen 
blast flame. Of those which seemed straight and 
uniform over a length of a foot or so, the finest 
(small K)" was taken. K for the first pendulum, 
with which the measurements of Fig. 5 were 
taken, was 8.52 10-* dyne-cm/radian, or 1.26 
Hence the 


4 


x10-* dyne-cm/scale centimeter. 
change in light torque in Fig. 5 was 3.9710 
dyne-cm. This first fiber was accidentally broken 
and most of the measurements were taken with 
the second fiber for which K was 10.110~-° 


dyne-cm/radian or 1.50X10-* dyne-cm/scale 


centimeter. The period in the latter case was 
about nine minutes, so that most of the deter- 
minations of x required consecutive observation 
of end-points for practically two hours. I made 
about 80 determinations during July, August 
and September, 1935, and Mr. Wilbur Harris 
made about 40 determinations during October, 
November and December, 1935, after I had left 
Princeton. Each of us worked in the absence 
of the other, and each made determinations of 
the three types to be described. Mr. Harris 
used a different bottom plate from the one 


1! With regard to the choice of the torsion constant K 
(in dyne-cemtimeters per radian) for the fiber, the follow- 
ing is of interest. During a ‘‘resonance”’ interval the ampli- 
tude (neglecting damping) increases by an amount 2 Nx in 
N cycles. For a given torque change, x is inversely propor- 
tional to K. The number of cycles in a given time is directly 
proportional to ¥ K, assuming the moment of inertia of 
the pendulum isa constant. Hence the change in amplitude 
2Nx in a given time is inversely proportional to y K. In 
other words, the time required to achieve a given increase 
in amplitude by the resonance method varies directly as 
v K, other factors being constant. Increasing K in this 
experiment would therefore involve an increase in the 
time required for each measurement, a disadvantage which 
was not present in the corresponding case of the Einstein- 
de Haas experiment, in which the impulse per half cycle 
was independent of the length of the cycle. Nevertheless, 
Professor Einstein has pointed out that it would also be 
worth while in this experiment to increase K considerably, 
in spite of the increase in time required to achieve a given 
change in amplitude, because then the ratio of the energy 
in the effect desired to the energy of random disturbances 
would become much more favorable. Furthermore, it is 
an advantage to have the end-point readings a, come at 
shorter intervals insofar as the possible segregation of 
the larger random disturbances is concerned. 
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Fic. 6. Comparison of theory and experiment. Type I 
measurements: light torque as a function of filament 
temperature. 


which I used, which seems to eliminate the 
possibility that the results, which were through- 
out consistent, could have been influenced by 
irregular patches of ‘‘infrared dirt’’ on the clear 
crystal quartz plates. Mr. Harris also investi- 
gated very thoroughly whether any effects were 
obtainable in the absence of the bottom plate, 
and found absolutely nothing. These checks 
were invaluable in substantiating the results of 
the present experiment, and Mr. Harris deserves 
much credit for the thoroughness and enthusiasm 
with which he carried out the tedious observa- 
tions. 

Three types of measurements, presented in 
Figs. 6, 8 and 9, were taken to show that the 
effect varies as required by the theory. The 
crosses represent my determinations and the 
circles those of Mr. Harris. Each point is derived 
from observations over a period of one and one- 
half to two hours as described in connection with 
Fig. 5. The smooth curves give the corresponding 
theoretical variation in each case. 

For type I (Fig. 6) the torque change was 
measured for various filament temperatures as 
determined by the current through the lamp. 


A. BETH 


The bottom plate was turned between stops at 
6=45° and @= —45°, and the plane of polariza- 
tion of the Nicol was at 45° to the axes of the 
plates as in Fig. 3. Mr. Harris’ values (circles), 
taken consecutively, without shifting the posi- 
tion of the lamp, form the smoothest sequence 
obtained. The scattering of my values (crosses) 
taken with a different bottom plate is partly 
explained by the fact that I did not make these 
determinations consecutively, but shifted the 
lamp position a number of times for other 
readings. The total amount of light entering the 
plate system changes rather sharply with the 
position of the lamp in its socket and the setting 
of the double slide because of the presence of the 
small light shield d in the middle of the window 
W. The measurements previously reported" of 
this type, were made with a different lamp and 
fiber, and were made with a somewhat larger 
light shield on the window W. The values were 
therefore somewhat smaller than those shown 
in Fig. 6, but the magnitudes increased in a 
similar way with temperature. 

The torque to be expected at various filament 
temperatures was calculated from (3) and the 
spectral energy distributions as follows. From 
the Poynting vector, the energy flow in the plate 
M is A?(n, cos? 6+, sin? @)c/87 or, practically, 
A’nc/8x=2nc(A/4r)*n ergs/cm?-sec. Hence, if 
J,'dd is the total energy flow (integrated over 
the area of the plate) in W/ in the wave-iength 
range \ to A+dxX, the net torque on the plate 
from the energy flow in this wave-length range 
will be L,’d\ dyne-centimeters, where 

L)’ =sin 26,'(sin 2A,;—sin 2As)AJ}//2 2. 


Here tan @,’ is, in general, the ratio of the 
amplitudes of the components of the light in the 
directions of the axes of the plates JT and M. 
6,’= 8 only when the position of the bottom plate 
B is such that the directions of its axes coincide 
with those for M and T. If the total energy flow 
entering the plate system in this wave-length 
range is J,dd ergs/sec., we set J)’=J,F), for 
each of the light beams contributing to the effect 
(enumerated by the index 7), and compute F,, 
using Edison Pettit’s values for the reflectivity 
of aluminum (top surface of 7)" and the re- 


2 R. A. Beth, Phys. Rev. 48, 471 (1935). 
13Edison Pettit, Pub. Astronom. Soc. Pacific 46, 27 
(1934). 
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fractive indices of quartz, which enter into 
Fresnel’s expressions for the transmitted and 
reflected intensities at the plate surfaces. The 
absorption in the quartz plates may be neglected 
beczuse the light has been filtered through the 
fused quartz lens L and the calcite and Canada 
balsam of the Nicol prism. From the theory for 
elliptically polarized light, the angles 6,’ and 
the phases 24 may be calculated for each beam 
and each wave-length and for any position of 
the bottom plate, because the incoming light 
from the Nicol is plane-polarized, and the plate 
thicknesses, as measured by the compensator in 
wave plates for sodium light, were: 7=0.649, 
MJ=1.051, B; =0.535 (my torque measurements), 
B,.=0.454 (Harris’ torque measurements). The 
whole torque on .V is then 


= J L)J)dd dyne-centimeters, (8) 
where the torque factors L) are given by 


L,=>-\F) {sin 26,'(sin 2A;—sin 2As)}y 2nc 
a ” ” 
od seconds. (9) 


The torque factors L, were calculated for the 
sum 7 of the effects of the main beam going 
upward through VW, the main beam reflected 
downward from the aluminized top surface of 7, 
and for the reflections of these beams at the 
surfaces of B, W, M, and 7, repeated surface 
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Fic. 7. Torque factor as a function of wave-length. 
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reflections being neglected. These values of (9) 
are shown in Fig. 7, the angles denoting the 
azimuthal position of the bottom plate; at 45° 
the direction of the ‘‘fast’’ axis of B coincides 
with the direction of the “slow”’ axis of M and 
the “‘fast’’ axis of T (see Fig. 3). The full curve 
refers to the calculations for Bo, the dotted curve 
to those for By. 

Using a rocksalt spectrometer and thermo- 
couple arrangement Mr. Harris compared the 
amount of radiation from a black body at a 
known temperature (radiation from the interior 
of a furnace made out of an alundum cylinder 
23 inches in diameter and a foot long at 766°C 
and 784°C was used) with the amount of radi- 
ation coming through the Nicol in the same 
solid angle at intervals of 0.1 uw for each of six 
different lamp currents. The energy factors J) 
were computed in absolute terms from these 
measurements using calculated values of black 
body radiation, Lambert’s law, and constants of 
the optical system to determine the solid angle 
of the cone of light entering the plate system 
from the Nicol. 

Replacing the integral in (8) by the corre- 
sponding sum for 0.1 uw intervals from 1.0 yu 
to 2.0 uw inclusive, the change in total torque 
from @6=45° to 6= —45° for each of the six fila- 
ment currents and for the bottom plate Bz was 
obtained. These theoretical torque values are 
plotted against the temperatures, 1948°, 2069°, 
2297°, 2506°, 2701° and 2879° absolute, corre- 
sponding to the filament currents used, in Fig. 6 
and connected by straight lines A. The irregu- 
larity of this broken line A indicates the dif- 
ficulties encountered in making the energy 
determinations. The principal sources of error 
are: The spectrometer was not set up especially 
for this purpose, but was originally designed for 
somewhat longer wave-lengths; the solid angle of 
the light cone is difficult to determine with 
accuracy; the neglected effect below 1.0 yu is 
probably not exactly compensated by assuming 
that the Nicol actually polarizes all the light out 
to 2.0 uw. There is reason to believe, from the 
energy measurements themselves, that the ener- 
gies computed for the upper two temperatures 
are somewhat too low. The root-mean-square 
deviation of the measured torques (circles) from 
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the computed values A is 0.29X10~-° dyne- 
centimeter. 

A second attempt was made to get the values 
of the energy factors J, in absolute terms by 
extrapolating the tungsten emissivities given by 
Forsythe and Worthing'* (assuming the emis- 
sivity to be a linear function of temperature at 
each wave-length), multiplying by the black 
body radiation at the temperature of the fila- 
ment, taking account of the solid angle of radi- 
ation collected from each point on the filament 
by the lens L, and making an allowance for the 
reflections from the surfaces of the bulb of the 
lamp, the lens, and the Nicol prism. The change 
in torque at each temperature was calculated 
from (8) as before, the integral in this case being 
replaced by the corresponding sum for 0.2 yu 
intervals from 0.8 un to 2.0 yu inclusive. The 
results must be reduced by about 15 percent to 
give the curve B in Fig. 6. This 15 percent may 
be caused by the absorption of the Canada 
balsam in the Nicol, by nonpolarization of the 
light near 2.0 u, and by inaccuracies in the 
emissivities used. The root-mean-square devi- 
ation of the circles from the curve B as drawn 
is 0.24X10-° dyne-centimeter. It should be 
emphasized that the uncertainties in B leave the 
absolute value in considerably greater doubt than 
in the case of A. 

An important conclusion from Fig. 6 is, 
however, that the slope and shape of the theo- 
retically determined curves A and B agrees 
within the limits of error with the trend of the 
measured torque values. The trend of B relative 
to the torque values (measured for certain lamp 
currents) depends on the current-temperature 
calibration of the lamp which was extrapolated 
from one value at 2941° absolute by comparison 
with the calibration for a similar lamp, while 
the trend of A relative to the torque values 
measured does not depend on this calibration, 
the energy distributions having been measured 
for given lamp currents also. The dotted smooth 
curve C, empirically fitted to the upper nine of 
Harris’ torque measurements, appears as a 
natural compromise between the trends of A 
and B. Half of the circles lie within one percent 
of this curve. Their root-mean-square deviation 
from the curve is 0.054 10~° dyne-centimeter. 


44 Forsythe and Worthing, Astrophys. J. 61, 151 (1925). 
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Fic. 8. Type II measurements: variation of light torque 
with angle between plane of polarization of light and the 
axes of the plates. 


If we choose to think of C as an empirical cali- 
bration curve, this indicates that we may deter- 
mine the temperature of the lamp in this range 
with a probable error of about 7° by measuring 
the torque with the accuracy attained by Mr. 
Harris, the lamp position being left unchanged. 

The type I measurements may be regarded as 
a test of the factors \J,/27c in the calculation 
of the torque from (8) and (9). 

For the type II measurements (Fig. 8) the 
lamp current, and therefore the light intensity, 
was maintained constant and the plate Bz was 
rotated between the same stops as in type I, but 
the plane of polarization of the light was rotated 
to various angles @ with the axes of the plates. 
These measurements are intended to test the 
factor sin 26 in the formula for the torque. 
The root-mean-square deviation of the meas- 
ured values from the computed sine curve 
shown is 0.34X10-*° dyne-centimeter. A small 
systematic deviation from the sine curve 
may be due to the difficulty encountered in 
lining up the plate axes accurately in setting up 
the apparatus. In setting the angle @ the whole 
lamp housing and optical system, including the 
Nicol, was rotated and the double slide reset as 
described. Because of the presence of the light 
shield d this causes a greater scattering of the 
measurements than would be obtained if the 
lamp could be left fixed in position. 

The measurements shown were made by Mr. 
Harris. In my measurements of this type only a 
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“1G. 9. Type III measurements: light torque as a function 
of position of the right stop. 


quarter cycle of a sine curve was covered, but at 
15° intervals. These points also form a smooth 
segment of a sine curve; though the root-mean- 
square deviation is considerably less, the test is 
not as critical and the result not as interesting as 
the full cycle curve shown. 

For the type III measurements (Fig. 9) the 
light intensity was again maintained constant 
and the plane of polarization of the optical 
system was again put at 45° to the directions of 
the axes of plates / and T as in type I. The left 
stop was fixed at —94}° and measurements made 
for different positions of the right stop. Since the 
torque for B, at the left stop was constant, the 
curve shows how the light torque in this plate 
system varies with the position of the bottom 
plate. The elliptical polarization produced by the 
bottom plate depends on its position, and these 
measurements are, therefore, in the main a test 
of the phase factors (sin 24;—sin 242) in the 
formula. From the torque factors of Fig. 7 the 
smooth curve shown in Fig. 9 was calculated 
from the energy distribution. The root-mean- 
square deviation of the measured values of the 
torque from is 0.073X10-* dyne- 
centimeter. This very good agreement is again 
to be expected because the lamp position was 
left fixed. Mr. Harris also made measurements 


this curve 
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of this type, but did not cover a whole cycle. 
They are in similar agreement with the predicted 


values. 
(CONCLUSION 


The change in the angular momentum of 
polarized light in passing through a crystal plate 
has been detected and measured with the appa- 
ratus described. This conclusion rests on the fact 
that the effect 
within the limits of error, in magnitude with the 


measured agrees in sign and, 
effect as predicted both by the wave and by the 
quantum theories of light. In particular it has 
been shown that the effect varies as required by 
the theory with each of three experimental 


variables studied: filament temperature (light 
intensity), azimuth of the plane of polarization, 
and position of the bottom plate. 

I am deeply indebted to Professor A. W. Duff 
of the Worcester Polytechnic Institute who first 
suggested this problem to me and who personally 
helped me with the purchase of apparatus in the 
preliminary stages of the investigation. I wish 
to thank the National Research Council for a 
Grant-in-Aid for the purchase of apparatus and 
the Palmer Physical Laboratory of Princeton 
for providing all the specially made apparatus 
and technical assistance as well as for the hos- 
pitality of the Laboratory. In particular I wish 
also to thank 
Condon, Robertson, Shenstone, Bleakney, and 
Dr. Barnes of Princeton for their willing helpful- 


Professors Harnwell, Einstein, 


ness in discussing the work with me on numerous 
occasions. Dr. S. N. Van Voorhis and Messrs. 
C. W. Curtis and C. W. Lampson of Princeton 
helped me both with suggestions about experi- 
mental details and in actually taking readings. 
Mr. Curtis shared the work of taking readings 
with me for many weeks. The very valuable part 
of the work done by Mr. W. Harris has been 
noted in detail above. Finally, I want to thank 
Mr. A. H. S. Holbourn of the Clarendon Labora- 
tory, Oxford, for his interesting correspondence 


with me about the problem. 
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It has been possible, by a judicious choice of field 
strength to separate the magnetic components due to the 
199 isotope of mercury. Five lines have been studied, 
AA3125(°P; D2), 5461, 4358, 4047(°P 2S). A practically 
complete Paschen-Back effect is observed, each ordinary 
Zeeman component being accompanied by two satellites 


HE hyperfine structure of the mercury lines 

has been studied by numerous investi- 
gators, particularly Schiiler and his collabo- 
rators.' They have been able to show that the 
very strong principal component (which is 
complex under very high resolution) can be 
attributed to the even isotopes of mercury, 
while the weaker satellites are due to the two 
odd isotopes, 199 and 201. A careful analysis 
seems to prove that the complexity of satellites 
can be correctly accounted for if we assign a 
nuclear spin of } to the 199 isotope and a nuclear 
spin of 3/2 to the 201 isotope. 

The effect of the application of an external 
magnetic field to a line showing hyperfine 
structure is to increase greatly the complexity 
of the pattern, each line giving its Zeeman 
pattern to the total; but as the field increases, 
and finally becomes very strong, certain simplifi- 
cations occur which may be described briefly as 
follows.” 

For very strong magnetic fields the energies of 
the various levels may be represented by 
(A <X<w) 


W=mjgwt+Amjm;, (1) 


the first term on the right side of the equation 
being the ordinary Zeeman level without hyper- 
fine structure and the second term due to the 
hyperfine structure. This means that instead of 
each single Zeeman level there will be 27+1 
levels (corresponding to the 27+1 values of m,) 


! See Schiiler and Jones, Zeits. f. Physik 77, 801 (1932) 
for a complete list of references. 

2 For a complete account, see Green and Wulff, Zeits. f. 
Physik 71, 593 (1931). 


due to the nuclear spin i=} of Hgisg. There is excellent 
agreement between the theory and experiment even to the 
slight asymmetries observed. Some comments on the 
recent work at the French Academy express doubt of 


the reality of their findings. 


separated by amounts proportional to Am;, and 
the lines due to transitions between these levels 


will be given by 


Li) 


W, — W2= w(m jg: — mj, 22) ( 


+(A\m; —A2m;,)m 


with selection rules m;,=m;,+1 or m;,=m,, and 
Am;=0. The first term gives the position of the 
ordinary Zeeman component, and the second 
term shows that instead of this one component, 
there will be 27+1 components separated by 
amounts 
A\m;,—A,m 

Thus, in the case of an element containing only 
one isotope, the Paschen-Back effect of hyperfine 
structure gives an immediate means of de- 
termining 7, the nuclear spin, provided A is 
large enough to allow the resolution of the 
separate components. This was the case for TI 
and Bi where one isotope predominates. But in 
the case of Hg, where we have a large number of 
isotopes each with its own nuclear spin, it is 
not so easy to separate the components belonging 
to each of them. Whereas in Tl, where 7=3, 
each of the ordinary Zeeman components yielded 
2i+1 or 2 components; in Hg, we should find 
the ordinary Zeeman components due to the 
even isotopes, 2 components due to isotope 199, 
and 4 components due to isotope 201. 

But the hyperfine structure constants A of 
Hgoo, are much smaller than those of Hgj99 and 
the four components of each line of Hgo are 
therefore crowded together and practically un- 
resolvable by the grating, and should merge with 
the lines due to the even isotopes, whose isotopic 
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displacements are also extremely small, and 
should appear as one line. The pattern to be 
expected is therefore the ordinary Zeeman 
pattern, but each line accompanied by a satellite 
on each side, due to isotope 199. 

An examination of Eqs. (1) and (2) shows that 
while the first term is proportional to the field 
strength, the second term depends only on the 
hyperfine structure constant. To be sure, this 
is true only when Aw, but experimentally, the 
condition is reached for comparatively low field 
strengths. The picture will then appear, as we 
increase the field strength, of a group of lines 
whose separations increase proportionally with 
the field strength (the ordinary Zeeman com- 
ponents) each accompanied by two satellites 
whose separation remains almost constant. Thus, 
unless the field strength is properly chosen, the 
satellites of one particular line may actually be 
beyond one of the other lines of the pattern. 
Furthermore, the satellites are quite faint in 
comparison with the principal lines, so that the 
latter must be very strongly overexposed in 
order to bring out the weaker satellites. This 
always involves the difficulty that in developing, 
the overexposed images tend to spread and 
unless the faint satellites are well-separated 


from the principal lines, they will be lost in the 
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background. It is therefore necessary to choose a 
field-strength and lines with such values of the 
hyperfine structure constant A, that the satellites 
fall about midway between the main Zeeman 
components. This is not possible, of course, for 
all of the components and all of the satellites 
simultaneously, but the conditions have been 
fulfilled in a sufficiently large number of cases to 
make the results worth reporting. 

A field strength of about 36,000 gauss was 
used, (w=1.70 cm™') and the *S*P triplet 
(AA5461, 4358, 4047) was studied, as well as 
§P, *Dz (43125) which is particularly interesting 
because of the fact that the *D, has a negative 
hyperfine structure constant.® 

The results are shown pictorially in Figs. 1 to 4. 
The upper half of each figure is a microphotogram 
of the line, and below is the pattern as calcu- 
lated from the exact formulas of Goudsmit and 
Bacher.‘ The exact calculations show asym- 
metries not mentioned in the discussion above, 
and these asymmetries are apparent on both 
the observed and calculated pictures shown in 
the figures. That the principal lines are all over- 
exposed is shown by the fact that all the principal 

>The apparatus and method used have been described 


by Green and Loring, Phys. Rev. 46, 888 (1934). 
*Goudsmit and Bacher, Zeits. f. Physik 66, 13 (1930). 
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peaks are of the same height except for the 
weakest of the The satellites 
appear on the microphotograms as subsidiary 


lines patterns. 
peaks where they are well-resolved and as points 
of inflection where the resolution is not very 
good. Of course, all the satellites do not appear. 
Some of them are too close for the resolving 
power available. 

The Zeeman effect of the Hg lines is being 
studied in detail at the French Academy of 
Sciences’ using their large magnet which pro- 
vides fields up to about 60,000 gauss. Certain 
irregularities seem to show up, but there is a 
doubt in the author’s mind of the reality of 
these, principally because of the effect of the 
unresolved h.f.s. components. If, for example, we 
label the principal components of \4358 with 


the numbers 1 to 6, the components 2 and 5 


I 
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x is. 


separate proportionally with the field strength, 
while 1 and 6 seem to separate more rapidly. 
This would seem to indicate that something was 
affecting *S, but not *P;. The same sort of thing 
happens with 5461, indicating perhaps a similar 
cause for the behavior of *P». The meagerness of 
the report makes it very difficult to comment 
upon, but in view of all of the work done on 
h.f.s., it would seem that any discussion of the 
peculiar behavior of the Zeeman components of 
a line cannot entirely neglect the effects of the 
Paschen-Back effect of the h.f.s. components, 
which tend to spoil any effect that concerns 
itself only with the linearity of position and 
field-strength. 

In conclusion, the author wishes to express 
his thanks to Professor R. A. Loring, of the 
University of Louisville, for his assistance in 
part of the experimental work connected with 


this paper. 
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The Diffusion and Absorption of Neutrons in Paraffin Spheres* 


> Dupouy and Jacquinot, Comptes rendus 201, 543 
(1935). Jacquinot and Belling, ibid. 201, 778 (1935). 
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Measurements are described of the activities imparted 
to silver and to indium detectors placed in paraffin spheres 
of various sizes with a radon-beryllium source of neutrons 
located at the center. The results are discussed with the aid 
of diffusion theory. Reasonably good quantitative agree- 
ment is obtained at distances from the source greater than 


OLLOWING the discovery of Fermi! and 
his co-workers that hydrogen containing 
substances are effective in slowing down the 
neutrons from a radon-beryllium source, several 
investigators have made experiments to deter- 
mine the velocities of the slow neutrons pro- 
duced. The results of various experimenters!: ?: * 

* Contribution No. 545. 

! Fermi, d’Agostino, Rasetti and Segré, Proc. Roy. Soc. 
A146, 483 (1934); Amaldi, d’Agostino, Fermi, Pontecorvo, 
Rasetti and Segré, Proc. Roy. Soc. Al49, 522 (1935). 

? Westcott and Bjerge, Proc. Camb. Phil. Soc. 31, 145 


(1935); Westcott and Niewodniczanski, Proc. Camb. Phil. 
Soc. 31, 617 (1935); Bjerge and Westcott, Proc. Roy. Soc. 


A150, 709 (1935). 
3 Moon and Tillman, Nature 135, 904 


1935). 


about 8 cm. Various independent estimates are made of 
the mean free path of the neutrons in paraffin, their 
specific absorption rate constant, and the ratio of the latter 
to the diffusion constant. The mean free path is found as 
0.53 cm; the ratio of the specific absorption rate constant 


to the diffusion constant as 0.026 cm 


have shown that the velocities are comparable 
with thermal velocities. Dunning‘ and others, 
by means of a velocity filter, found that the 
neutrons slowed down by paraffin had a velocity 
distribution at room temperature substantially 
that predicted by Maxwell’s law. These results 
suggested that, under definite and simple geo- 
metrical conditions, the behavior of slow neu- 
trons in paraffin might be more closely studied 
from the point of view of diffusion. The method 
of experimentation adopted consisted in deter- 


‘Dunning, Pegram, Fink, D. P. Mitchell and Segré, 
Phys. Rev. 48, 704 (1935); Fink, Dunning, Pegram, and 
Mitchell, Phys. Rev. 49, 103 (1936). 
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mining the activity produced in indium and 
silver detectors placed inside paraffin spheres of 
different diameters at measured distances from 
a source located at the center. 


EXPERIMENTAL PROCEDURE 


Four paraffin spheres were used; their diam- 
eters were 20, 25, 35, and 51 cm, respectively. 
The smaller ones were spherical to within two or 
three millimeters and the largest to within five 
millimeters. A cylindrical hole, 7 cm in diameter, 
was bored radially to the center of each sphere, 
and fitting easily into this hole was a stack of 
paraffin cakes. The two larger surfaces of each 
cake was so curved that a thin flexible detector 
placed between any two cakes would have all 
points on its surface equidistant from the source 
located at center of the sphere. 

The detectors were thin, circular disks of silver 
and indium with several small sectors cut out in 
order to make them conformable to the curved 
surfaces of the paraffin cakes. The thickness of 
the indium detector was 0.005 cm and that of 
the silver 0.02 cm. The spheres and detectors 
are shown in Fig. 1. 

A 3-mm X5-mm cylindrical brass capsule con- 
taining powdered beryllium and radon (300 
millicuries or less) served as the source of 
neutrons. 

In an experiment the detector used was left in 
some one position for a time sufficient for the 
activity to become practically saturated (7 to 
12 half-lives). The activity was measured with a 
quartz fiber electroscope. Readings were begun 





a b 


Fic. 1. (a) Cross section of the 10-cm paraffin sphere 
showing arrangement of the cakes between which the de- 
tector was placed. (b) The detector drawn to the same 
scale. 


at a fixed time (Ag, 30 sec.; In, 60 sec.) after 
removal of the detector from the sphere and 
were continued at suitable intervals for a definite 
period of time (Ag, 6 min.; In, 15 min). The 
total reading for the period, after correcting for 
nonlinearity of scale, for background, and for 
radon decay, was taken as the activity, 8, for 
that experiment. 

The half-lives of the two activities of silver 
were determined with the electroscope used and 
were found to be 23.5 sec. and 2.15 min., re- 
spectively. 


RESULTS OF THE EXPERIMENTS 


The results of the experiments are shown 
graphically in Figs. 2 and 3 where @r is plotted 
against r the distance of the detector from the 
center of the sphere in cm. Because of the strong 
activity and relatively long half-life of indium 
the best accuracy was obtained with it as the 


detector. 


DISCUSSION OF THE RESULTS 


If it is assumed that the neutrons diffuse 
through the paraffin according to the ordinary 
diffusion laws, and that the absorption by the 
thin detector does not cause appreciable dis- 
turbances, and further that the neutrons are 
absorbed in the paraffin at a rate, ac, proportional 
to the concentration, then 


DV*c —ac=dc/dt 


where D is the constant diffusion coefficient and 
a is the specific absorption rate constant. Under 
the steady state conditions of the experiments 
dc/dt=0. The equation becomes then 


DV*c—ac=0. (1) 


The general solution, for the present case of 
spherical symmetry, may be written in either of 
the forms 


cr=(A/a) sinh ar+ Be-2" 
=B sinh a(r+k)/sinh aki; a&’=a/D. (2) 


A, B and hk, are constants of integration and 
are related to each other by the equation 
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Fic. 2. Plots of 8r vs. r for the indium detector in spheres 
of radii 10, 12.5, 17.55, and 25.55 cm, respectively, from 
left to right. The curves for the spheres of various radii 
were drawn using Eq. (5). The top curve is for a sphere of 
infinite radius. 


A=aB/sinh ak. B is directly proportional to 
the strength of the source when the source is 
placed at the center of the sphere. In case a=0, 
cr=Ar+B. 

There are two particular solutions contained 
in (2) which will be required later, namely, 


cr=Be-*" (3) 


when the source is central and the sphere is of 
infinite radius, and 


cr=(A’/a) sinh ar (4) 


when the source is uniformly distributed over 
the surface of a finite sphere. 

In order to relate the concentrations in one 
sphere with a central source to those in a sphere 
of different radius a boundary condition may 
be introduced as follows. The net outward flux 
of neutrons across a spherical surface just inside 
the paraffin is equal to the rate of escape of 
neutrons from the surface. The first quantity is 
equal to —4rr*Ddc/dr; if the rate of “‘evapora- 
tion” of the neutrons per unit area be assumed 
proportional to the concentration prevailing a 
short distance under the surface, then 


— D(dc/0r) par, = RoC rary, 


where fo is the radius of the sphere. The use of 
this boundary condition leads to the particular 
solution 


I 
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Fic. 3. Plots of Br vs. r for the silver detector. The solid 
curves are identical with those of Fig. 2. 


cr=B sinh a(r+l—ro)/sinh a(l—ro), (5) 
where /=(1/a) sinh~'[a/{(1/ro—k2/D)*—a?}!). 
The equations derived involve the concentra- 
tion of neutrons, while the measured quantity 
is B, the activity of a detector. The number of 
neutrons per unit time entering the detector 
will depend on both the concentration in its 
neighborhood and the velocity distribution; of 
the neutrons entering, the fraction captured by 
the detector may again depend on the velocity 
distribution.’ Proportionality between 6 and <¢ 
may safely be expected, then, only in regions of 
constant velocity distribution; this doubtless 
excludes regions close to the central source. 
Moreover, in regions where the velocity dis- 
tribution is not constant it is not likely that D, 
a, and ky remain constant as assumed in Eqs. 
(1) to (5). However, in order to see to what 
extent the observations are in harmony with the 
hypothesis, we have assumed Eqs. (1) to (5) and 
have placed 8=k’c, with k’ constant for a given 
detector. Similarly, with the detector placed on 
the surface of a sphere the activity 8) is assumed 
proportional to fo, the flux of neutrons from the 
surface : By = ko’ fo. If one assumes that the capture 
cross section of indium or silver is a constant, 
and further that the velocity distribution of the 
neutrons is Maxwellian, then it may easily be 


5 For recent experiments and theory on absorption see 
Rasetti, Segré, Fink, Dunning and Pegram, Phys. Rev. 49, 
104 (1936). Breit and Wigner, Phys. Rev. 49, 519 (1936) 
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TABLE I. Total fluxes measured by indium. 
SPHERE RADIUS, ro Activity, Bo (Mean) FLUX, ro%8 
10.0 18.7 1870 
12.5 9.83 1540 
17.55 2.61 800 
25.55 0.53 340 
shown that k’=bno(2kT/xm)* and ko’ =bno 


where 1 is the number of nuclei per cm? in the 
detector and 3d is an instrument constant. If, on 
the other hand, the capture cross section ¢ = a0/2, 
then k’=bnoo/2 and ko’ =(bn/2)(rm/2kT)'. In 
either case k’/ko’ =(2kT/xm)}. 

There appears to be some justification for the 
assumption of constancy of k’ in the fact that 
the curves for silver and indium are superposable 
throughout. Silver and indium are known to be 
activated by neutrons of different qualities.® 
Also the ratio of the 23.5-sec. activity of silver 
to the 2.15-min. activity remains substantially 
constant. It does not follow, but it does appear 
possible that the neutrons which cause activa- 
tion are all at the same ‘‘temperature,” and, 
from Dunning’s velocity filter experiments, the 
temperature corresponds to that of the paraffin. 
It is recognized, of course, that in the region 
near the source the velocity distribution is 
presumably not Maxwellian. 

On the assumption that k’, D and a are 
constants, independent of r, it may be immedi- 
ately concluded that, since the values for Br do 
not fall on a straight line, some absorption of 
the neutrons in the paraffin occurs. Moreover, 
since the total outward flux, 7928), as measured 
by detectors placed on the surfaces of the spheres, 
decreases with increasing radius of the spheres it 
follows that absorption takes place in the paraffin 
if ko’ remains constant after the neutrons have 
traversed 10 cm or more of paraffin. The total 
fluxes are presented in Table I. 

In order to compare the predictions of the 
theory with the observations in the region r=>8 
cm, the following procedure was used. Eq. (2) 
was fitted to a smoothed curve drawn through 
the points for the 51-cm sphere. A satisfactory 
fit was obtained with 0.16 cm as the value 
for a. Eq. (5) was then used to construct the 
predicted curves for the other spheres. The 


®C. H. Collie, Nature 137, 614 (1936). 
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agreement with the observations may be seen in 
Figs. 2 and 3. It is satisfactory for all but the 
25-cm sphere where the deviations are somewhat 
greater than the experimental error. Evidently 
then, had independent values of a been con- 
versely calculated from the relative activities 
near the surfaces of the various spheres, sub- 
stantially the same value of a as before would 
have been obtained. 

A value for a may be obtained in still another 
manner as follows. At a given distance r from 
the center, the activity produced in a detector 
increases with the radius 7 of the sphere. The 
increase must be due to neutrons which have 
been at larger distances and then diffused back 
to r. The difference in the activity (times 7) at r 
for two different spheres of radii 7) and 7’ will 
evidently be 

A’'(ro’) —A’’(r9) 


A= $'’r— r= 


sinh ar 
a 


and (1/A)(dA/dr) =a coth ar. (6) 


Application of relation (6) to the observations 
yielded an average value of 0.23 for a. This is 
appreciably higher than the value given above, 
but the errors involved are also considerably 
greater. For the present it can only be said that 
the two results are quite closely of the same 
order of magnitude. 

An estimate of the value of a was made in the 
following manner. The total rate of absorption 
of neutrons in a sphere is equal to /o"42r°acdr; 
and the total flux of neutrons from the sphere 
equal to 4rro?fy. The sum of these is equal to 
the rate of issue of neutrons from the source, so 
that for two experiments with equal sources but 
different spheres of radii 7» and ro’ the sums are 
equal. Assuming a constant 


ref = ro*fo 
a=— 


Sort dr— fo’ ’r’c'dr 
ry’ Bo’ —To*Bo k’ 


Sor Bdr — fo r’B'dr he’ 


The last member also assumes constancy of ky’ 
and k’. At first sight these assumptions seem 
especially objectionable since the integrations 
extend from the origin; however, at the smaller 
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TABLE II. Mean free path of neutrons in paraffin from indium 
measurements. 


ro(cm 8 0B dr \ 

10.0 18.7 48.9 0.58 

12.5 9.83 26.9 0.55 

17.55 2.61 8.3 0.47 

25.55 0.53 1.5 0.53 
Mean 0.53 cm 


r’s, B and 8 do not greatly differ, and the 
contribution to 
under 6 or 7 cm is at most only a few percent. 
From graphical evaluation of the integrals, 
values of a of from 0.99 to 1.22 10-°k’/ko’ with 
a mean of 1.1310-*k’/ky’ were obtained. With 
k’/ko’ equal to that for a Maxwellian distribution 
of 295°K, a becomes 1.41 X10°+0.08 sec.-. 
Finally, from a consideration of the behavior 


the denominator of distances 


of the neutrons at and near the surface of a 
sphere, a value of their mean free path \ and the 
diffusion coefficient, D=4(8k7T/xm,,)', may be 
calculated. By kinetic theory considerations it 
may be shown that the net flux, in the direction 
of increasing r, of neutrons through unit area of 
a detector placed normal to the concentration 
gradient dc/dr is —4dvdc/dr, where v is the mean 
velocity ; similarly the gross flux is }cv. The net 
flux is measured by a detector placed on the 
surface of the sphere and the gross flux by a 
detector just inside the surface. From this it 
can readily be shown that 


A= —(3/2)Bo/(08/dr), 


where the derivative is taken a short distance 
inside the surface. The values of \ obtained for 
the several spheres are given in Table II. The 
fact that \ is about the same for all spheres does 
not necessarily imply a constancy of the temper- 
ature of the slow neutrons. On the assumption 
that T=295°, D=4.4X 10+. 

From the independently determined values of 





DICKINSON 


ms See 


a and D, a value of a=(a/D)'=0.18 cm™! was 
calculated which is in very good agreement with 
the value 0.16 cm found above. 

The measurements also permit calculations 
for scattering, and cross 
sections for capture.’ That for 
may be obtained from the mean free path 
through the relation ¢,=3Nd where N is the 
number of scattering nuclei per unit volume. If 
hydrogen alone be assumed to scatter, then 
o,=18X10-** cm*. If, 
cross section of 4 10-** be allowed for carbon," 
that for hydrogen becomes 16 X 10-4. The cross 
section for absorption may be obtained from the 


of cross. sections 
scattering, 


Cs, 
however, a_ scattering 


specific absorption rate constant through the 
relation ¢,=(a/N)(rm,/8kT)! where m,, is the 
mass of the neutron. The results of Fleischmann® 
and of Kikuchi, Husimi and Aoki!® indicate that 
protons are mainly responsible for capture in 
paraffin. On capture by 
hydrogen alone, our value of a leads to ¢,=7.3 
X10-** cm’; Kikuchi obtained 3.0 to 8.3 10-*' 
then oy=a/N, inde- 


the assumption of 


cm*. If we take ¢,.=00/v 
pendently of any assumed temperature. The 
value calculated oo=1.8X10-°° cm? 
1, and this would be the cross section for 


becomes 
sec. 
neutrons of unit velocity. These values for o, 
and oo are presumably lower limits. 

The writers are indebted to Professor C. C. 
Lauritsen for helpful suggestions and the use of 
an electroscope, and to Mr. J. E. Walstrom for 
assistance in the preparation of the neutron 


sources. 


7 It has been emphasized by A. C. G. Mitchell that, in 
most cases, absorption and scattering cross sections have 
not been separated. Mitchell, Murphy and Langer, Phys. 
Rev. 49, 400 (1936). 

8’ Dunning, Pegram, Fink and D. P. Mitchell, Phys. Rev. 
48, 265 (1935). Dana P. Mitchell, Phys. Rev. 49, 453 
(1936). A. C. G. Mitchell, Murphy and Langer, Phys. 
Rev. 49, 400 (1936). 

* Fleischmann, Zeits. f. Physik 97, 242 (1935). 

1 Kikuchi, Husimi and Aoki, Proc. Phys. Math. Soc. 
Japan 18, 35 (1936); Nature 137, 30 (1936). 
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The scattering of slow neutrons from various materials 
rhe 


scattering cross section for a number of elements has been 


has been measured with several different detectors 


determined and the results given in a table. With suitable 
filters, the scattering curves for several of the Fermi groups 
have been investigated. The differences between the groups 
show up most markedly when the scatterer is itself a good 


1. INTRODUCTION 


N several recent papers': *:* we have reported 

on some investigations of the scattering of 
slow neutrons by various metals. The experi- 
ments were originally designed to see whether the 
elements showing large ‘“‘absorption cross sec- 
tions”’ also showed large scattering cross sec tions, 
since the original theories of interaction‘ pre- 
dicted such an effect. Experiments on the metals 
Ag, Cd and Hg as scatterers failed to show any 
large scattering and indeed the shape of the 
that 


sorption cross section was due to capture, the 


scattering curves indicated the large ab- 


scattering cross section being normal. These 


facts are in agreement with the theory of Breit 
Wigner” that 
capture are complementary. 


and who show scattering and 


Further experiments® with Rh and In detectors 
marked the 
curves for the same metal with various detectors. 


showed differences in scattering 
A large part of these differences can be accounted 
for by absorption in the detector, since many 
substances have been shown to exhibit marked 
absorption bands by Fermi and his collaborators. ' 
In other cases, discussed below, we have been 
able to show marked differences in the form of 
the scattering curves due to the various groups 
discussed by Fermi. 


' A. C. G. Mitchell and E. J. Murphy, Phys. Rev. 48, 
653 (1935). 
4. C. G. Mitchell, E. J. Murphy and L. M. Langer, 


Phys. Rev. 49, 400 (1936). 
A. C. G. Mitchell, E. J. Murphy and M. D. Whitaker, 
Phys. Rev. 49, 401 (1936). 


‘H. A. Bethe, Phys. Rev. 47, 747 (1935); E. Amaldi, 
©. d’Agostino, B. Pontecorvo, F. Rasetti and E. Segre, 
Proc. Roy. Soc. A149, 522 (1935); Perrin and Elsasser, 
Comptes rendus 200, 450 (1935). 

°G. Breit and E. Wigner, Phys. Rev. 49, 519 (1936). 

®E,. Amaldi and E. Fermi, Ricerca Scientifica; (VI) 2 


(1935). 
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absorber of one of the groups. The effect of temperature on 


the activation of In has been investigated, and the 
activity at liquid air temperature 
temperature ratio Joo 
activity at room temperature 
found to be 1.08 for the D neutrons and 1.21 for the ¢ 
neutrons. 
Finally, using the method of Moon and 
Tillman? and Rasetti and Fink,’ we have 


investigated with an indium detector the temper- 
ature coefficient of the C and D neutrons 
2. APPARATUS AND METHOD 

The apparatus used was the same as that 
A cylindrical block of 
paraffin, 15 cm in diameter and 17 cm high, 
contained a Ra-Be bulb located on the axis 6 cm 
from the top surface. The detectors were thin 
sheets of the metals Ag, Rh and In of dimensions 


previously employed.! 


6 by 10 cm. These were placed on top of the 
paraffin cylinder and the scatterer, usually blocks 
of metal the the detector, 
placed on top of the detector. In some cases a 


same size were 


as 


sheet of Cd metal, containing 0.218 ¢/cm~*, was 
interposed between the paraffin and the detector. 
This served as a filter to cut out neutrons of the 
C group. After irradiation for a given length of 
time the detectors were placed around a Geiger- 
Miiller tube counter 
with the help of a 


and the activity obtained 
of 


some experiments a “‘scale 


thyratron ‘‘scale two” 
recording system. In 
of four’’ counter was also used and gave results in 
substantial agreement with those obtained with 
the of 


communications the percent scattering, defined 


“scale two” counter. In the previous 
as the percentage ratio of the activity of the 
detector with a thickness x of scatterer minus the 
activity without scatterer divided by the activity 
without scatterer, was plotted as a function of 
the thickness of the scatterer. 

In view of the rather large variations in 


’P. B. Moon and J. Roc. Soc. A153, 
476 (1936). 
’F. Rasetti and G. 


Rev. 49, 642A (1936). 


R. Tillman, Proc. 


New York meeting, Phys. 


A. Fink, 
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absorption coefficient for the various groups of 
slow neutrons exhibited by the several detectors 
which we used, we measured their absorption. 
This was accomplished by using the actual 
detectors as absorbing sheets and measuring the 
fraction transmitted with the help of a detector of 
the same material. In certain cases a Cd filter 
(0.218 g/cm) was used to absorb the C neutrons. 
The time of irradiation of the materials was the 
same as that used in the scattering experiments 
(see below). The results of the absorption 
experiments are given in Table I. 

To get the true scattering fraction from the 
curve of percentage scattering, as defined above, 
it is necessary to correct the results for ab- 
sorption in the detector. If y is the fraction of 
neutrons absorbed, then 1—y is the fraction 
getting through to the scatterer. If a(x) is the 
fraction of these scattered from a thickness x of 
scatterer, then the number of these recorded by 
the detector is a(x)y(1—y). The percentage 
scattering as it is defined above, and has been 
used in our previous papers is a(x)(1—y)-100. 
Hence to get the true fraction of neutrons scat- 
tered the results must be divided by (1—y), the 
fraction transmitted by the detector. All curves 
published in this paper have been corrected and 
a corrected table of relative cross section is also 
given. 

The time of irradiation depended on the 
period of the detector used and was also governed 
to some extent by practical considerations. The 
In (54-minute period) detectors were irradiated 
for 15 minutes and readings were taken at 1, 5 
and 9 minutes after the end of irradiation. The 
Ag (22 sec., 2.2 min.) was irradiated 6 minutes 
and readings were taken at 1 and 5 minutes after 
irradiation; while the Rh (44 sec., 3.9 min.) was 
irradiated 5 minutes and readings taken at 1 and 
5 minutes also at 2 and 6 minutes after irradiation. 

In order to test whether there was a difference 
in the scattering as detected by the different 
periods of Ag, the following experiments were 
performed. The detector was irradiated for 10 
seconds and counting began 30 seconds after 
removal of the foil. In this manner scattering 
curves were run for Ni and C and gave sub- 
stantially the same results as were obtained for 
the combined periods. Furthermore, one run 


was made to measure the scattering of a certain 
thickness of Ni, in which the Ag detector was 


MURPHY 


AND WHITAKER 


raB_e I. Absorption measurements 
FRA 
D IR THICKNESS (g l TRANSM 
\g 0.105 0.93 
\g (filtered 0.105 0.73 
Rh 0.124 0.80 
Rh (filtered 0.124 0.55 
In 0.109 0.82 


irradiated 6 minutes and the counting was not 
started until 6 minutes after irradiation, so that 
the activity of the thirty second period had 
practically disappeared. This point lay on the 
curve determined by the other groups of measure- 
ments. The scattering from the metals Ni and C, 
with an Ag detector, is, therefore, independent 
of the period used in the detection. This is 
probably also true for the other scatterers with 
the possible exception of Ag. 


2 


3. RESULTS 


Figs. 1, 2, 3 and 4 show the scattering curves 
obtained for the elements Fe, Ni, Cu and C using 
the several detectors as shown. The form of the 
curves is such that they rise rapidly out of the 
origin and eventually reach a saturation value. 
This form of curve can be explained on greatly 
simplified assumptions as follows. Assuming that 
a plane parallel beam of neutrons falls on a 
scatterer of thickness x, then the fraction of 
neutrons scattered baek to the detector will be, 
if one assumes a combination of single scattering 
and capture, 

n Nao, 
= [1 —e72Nbostocdz), 
no 2Nb(o,+<¢,-) 


where N is the number of atoms/cm* of the 
scatterer, o, and o, the scattering and capture 
cross sections, respectively, and a and } quantities 
which come in because of the solid angles 
involved. The quantity N(o,+0¢,) is the usual 
absorption coefficient, and the factor 2 arises 
since the neutrons traverse the scatterer twice. 
Such a simple picture does not, of course, 
correspond to the facts since the source does not 
give a plane parallel beam of neutrons, and 
integration over the various angles involved will 
eventually have to be carried out. The calcu- 
lation does show, however, that the slope of the 
curve for small thicknesses is proportional to 
Ne, and the saturation value is governed by 
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a,/(o;+o.). Thus the relative scattering cross 
section can be obtained and some idea as to the 
ratio of scattering cross section to total cross 
The 


sections of various scatterers, obtained with a Ag 


section. values of the scattering cross 
detector, together with the saturation value of 
the curve is given in Table II. 

The curves of Fig. 1 show the scattering by Fe 
with Ag, In and Rh detectors with and without 
filtering through Cd. The curves for neutrons 
filtered through 0.218 g/cm? of Cd_ should 
represent the scattering of the A and B neutrons 
(Ag detector) and the D neutrons (Rh detector). 
It will be seen that these groups are scattered 
equally well from the iron. The curves for the 
unfiltered neutrons lie below those obtained 
when Cd filtering was employed and in the 
order of the number of C neutrons detected— 
Ag (50 percent C, 25 percent A, 25 percent B), 
In (60 percent C, 40 percent D), Rh (72 percent 
C, 25 percent D). This seems to indicate that the 
C neutrons have a somewhat smaller ratio of 
scattering cross section to total cross section in 
other The above 


iron than do the groups. 


explanation holds, of course, only in the case 
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that the velocity distribution of the neutrons is 
not changed on scattering from the iron. Such an 
effect in Fe does not appear to occur to any 
appreciable extent.! 

The curves for the scattering from Ni (Fig. 2 
show that the results with Rh and In detectors 
are the same, whereas the curve for the silver 
detector shows cereater scattering. In the case of 
scattering by Cu (Fig. 3) the curves for Ag, Rh 


) 


rasie II. Relative illtering cross section 


of various elemen 


Max 

ELEMENT At. No a, X10-™% cm SCATTERING 

Cc 6 3.4 61 

Meg 12 3.0 37 (at 6cm 

d 13 1.0 21 (at 6 cm 

S 16 0.9 18 (at 6 cm 

Cr 24 1.4 13 

Mn 25 2 17 

Fe 26 10.6 95 

Ni 28 18 100 

Cu 29 8.3 59 

Zn 30 3.7 40 

Ag 47 6.3 Q 

Cd 48 1.2 2.4 at 4 mm 

Sn 50 41 34 

lig sO 4.7 21 

Pb 82 7.7 56 

Bi 83 10.2 43 
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Fic. 5. Scattering of neutrons from 
silver. (Ag detector.) 


and In detectors give approximately the same 


results. For the experiments on the scattering by 


C, the curves for the Ag and In detectors are the 


same whereas the scattering curve using a Rh 
detector lies lower. No definite explanation of 
the differences in these scattering curves can at 
present be given on account of the complicated 
nature of the phenomena. Collisions with carbon 
atoms are known to cause a change in the 
velocity distribution of the neutrons,’ and this 
will probably be the governing factor in the 
explanation. 

The scattering from Ag and Hg present 
interesting cases. Both of these metals are good 
absorbers of ‘“‘slow’’ neutrons. The cross sections 
for absorption in the whole range of the “slow 
region” have been found to be!” for Ag 55 10-"4 
cm~?, and for Hg 38010-*4 cm-~?. Furthermore, 
it is known that Ag shows marked selective 
absorption. 

We have measured the scattering from Ag 
with an Ag detector with and without filtering 
through 0.218 g/cm-? Cd. The curves are shown 
in Fig. 5. It will be seen that when no filter is 
used the scattering is small, about 9 percent, 
and the curve reaches saturation quickly. This 
indicates that the ratio of scattering to total 
cross section is small, most of the neutrons being 
absorbed in the scatterer before getting back to 
the detector. When the Cd filter is employed, on 
the other hand, the maximum observed scattering 
is 25 percent and the slope of the curve indicates 
a larger ratio of scattering to total cross section. 
The neutrons responsible for this scattering are 
probably the B group since the C group is 


?E. U. Condon and G. Breit, Phys. Rev. 49, 229 (1936). 
G. A. Fink, J. R. Dunning and G. B. Pegram, Phys. Rev. 
49, 340 (1936). 

1 J. R. Dunning, G. B. Pegram, G. A. Fink and D. P. 
Mitchell, Phys. Rev. 48, 265 (1935). 


Fic. 6. Scattering of neutrons 
irom mercury. 


absorbed in the Cd and the A group in the Ag 
scatterer. 

The scattering from He with an Ag detector 
with and without filtering is shown in Fig. 6. 
When no filtering is employed the maximum 
scattering is 20 percent showing that a large 
fraction of the neutrons are absorbed in the 
scatterer. For those neutrons which have been 
filtered through Cd, however, a marked differ- 
ence in the form of the curve exists. The maxi- 
mum observed scattering is 55 percent and the 
curve is still rising at this point. The curve shows 
that the ratio of scattering to absorption cross 
section is larger in the case of the filtered neutrons 
(presumably groups A and B) than for the 
unfiltered. It is apparent, therefore, that the 
absorption coefficient is large in the region of 
thermal velocities (C region) and is not so large 
at those points in the higher energy region which 
correspond to the A and B group. This does not 
mean, of course, that the absorption may not be 
large in some other region of the spectrum not 


registered by a silver detector. 


4. TEMPERATURE EFFECT FOR THE NEUTRONS 
Wuicuh Errect INDIUM 

We have measured the temperature effect on 
the neutrons which effect In by the method of 
Moon and Tillman’ and Rasetti and Fink.’ A 
block of paraffin 22 cm in diameter and 15 cm 
high was hollowed out so that a Dewar flask 
could be placed along the axis of the block. 
Inside of the Dewar there was located a cylinder 
of paraffin of 1 cm wall thickness, on the inside 
of which a cylindrical In detector could be 
placed. Arrangements were also made to inter- 
pose a cylinder of Cd (0.410 g/cm) between 
the detector and the small paraffin cylinder of 
paraffin. The neutron source was embedded in 
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ase III. Effect of liquid air temperatures. 
\RRANGEM ( Ss 
Liquid air, no Cd 52.816 
Liquid air, Cd filter 8,027 
No liquid air, no Cd 44.309 
No liquid air, Cd filter 7,473 


the large block 9 cm off the axis so that the 
thickness of paraffin in a direct line between the 
source and the Dewar was 4+ cm. The detector 
was activated when the inner paraffin cylinder 
was at liquid air temperature and when it was 
at room temperature, with and without the Cd 
filter. The In was irradiated 45 min. and counting 
was carried out from the fifth to the eightieth 
minute after irradiation. The results are shown in 
Table III. 

Dividing the results obtained at liquid ait 
temperature by those obtained at room tempera- 
ture, when the Cd filter was present, gives the 
temperature ratio for the D group. This value 
turns out to be 1.08, and is the same as that 
obtained by Rasetti and Fink using Rh as a 
detector. The differences between the counts 
obtained with and without the Cd filter give 
that part of the activity due to the C group. The 
ratio of the activities due to the C group is 1.21. 

Recently Rasetti, Fink, Goldsmith and 
Mitchell'' have determined the energy of the 
various groups of neutrons with boron absorp- 
tion. They showed that the energy of the D 
group activating Rh is approximately 0.6 ev, 
while the energy of those neutrons accounting 
for that part of the activity of In which is not 
effected by Cd filtering is about 0.35 ev. This 
small difference in energy between the two 
classes of neutrons would presumably now show 
up in the measurement of the temperature effect. 

The authors are indebted to Mr. L. M. Langer 
for help with the readings and to Mr. C. B. 
Braestrup of the Physical Laboratory of the 
Department of Hospitals of the City of New 


York for many favors. 


Vote Added with the Proof. R. Fleischmann" has recently 
published experiments in which he obtains much lower 
results for the scattering of slow neutrons from Fe and Cu 


than we do, and criticizes our results on the following 


1 F. Rasetti, G. A. Fink, H. H. Goldsmith and D. P. 
Mitchell, Washington meeting, Phys. Rev. 49, 869A 
1936). 

2 R, Fleischmann, Zeits. f. Physik 100, 307 (1936). 
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grounds. He supposes that in our neutron beam there are a 
considerable number of fast neutrons which do not affect 
a silver detector on their way from source to scatteret 
rhese fast neutrons are then supposed to be scattered in 
the Fe or Cu, return through the detector to the paraffin, 
and there suffer reflection with a loss of velocity which 
makes them capable of activating the silver detector on 


their second trip toward the scatterer hat this effect 


is playing no predominant role in our experiments can be 
shown by several experiments which we have already 
published but which we shall discuss again, in more detail, 
to allay this criticism 


munication, 


1 a group of experiments, described in our first com 


we showed that, using an Ag detector, 4.3 
cm of Cu as scatterer caused an increase of 52 percent in 
the activity of the detector. In a second experiment, a 
sheet of Cd 1 mm thick (0.86 g/cm? 


source from the detector. With this arrangement about 25 


shielded the neutron 


percent as many counts were obtained with Cd as without 
\ thickness of 4.3 cm of Cu caused an increase of about 50 
percent in the count 

In another series of experiments the Ag foil was not 
shielded by Cd from the source, but a sheet of Cd 1 mm 
thick was placed on top of the foil. With this arrangement 
no change of activity of the detector was noted. With 4.3 
cm of Cu placed above the Ag and Cd there was only an 
11 percent increase instead of the 52 percent obtained with- 


out Cd. To be explicit we can use the following figures, 


taken directly from experiment, calculating the activity of 


the detector to be 1000 (the actual count was 1043 


Percet 
\rrangement Activit Scatterir 
Source, Ag 1000 
Source, Ag, 4.3 cm Cu 1540 54 
Source, Cd, Ag 253 
Source, Cd, Ag, 4.3 cm Cu 392 55 
Source Ag, Cd 1000 
Source, Ag, Cd, 4.3 em Cu 1100 11 


rom the results of the experiments with and without the 
Cd between source and detector, one can calculate the 
result to be expected if the Cd is placed between detector 
and scatterer, assuming that the second reflection in the 
paraffin is negligible. It has been shown that, for those 
neutrons which activate Ag, the fraction transmitted by 
Cd is constant beyond 0.3 g¢/cm?,*:' so that the only ab 
sorption to be considered takes place when the neutrons 
pass through the Cd on their way to the scatterer. The 
absorption on the way back to the detector may be 
neglected. Using the figures given above one would expect 
an activity of 1000+139=1139 for the arrangement 
source, Ag, Cd, Cu), or an increase of 12.5 percent in 
substantial agreement with experiment." If, on the other 
hand, a large number of fast neutrons were scattered from 
the Cu, returned to the paraffin, and were reflected back 
to the detector at a lower velocity, one would expect a 
much greater activity tor the above arrangement than was 
actually found. 

'8 A.C. G. Mitchell and E. J. Murphy, Phys. Rev. 47, 
S881 (1935). 

'* That no appreciable change of velocity takes place on 
scattering by Cu has been shown by us. (See ref. 1, p. 655.) 
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Inelastic Scattering of Electrons from Solids 
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The energy distribution of electrons scattered inelasti- 
cally from solid targets has been studied for primary 
energies from 50 to 400 volts, using a method of magnetic 
deflection. Target surfaces of Cu, Ag and Au were deposited 
by evaporation in vacuum immediately before the measure- 
ments were made. For certain values of the energy loss 
suffered by the scattered electron the curves exhibit 
maxima, which are characteristic of the material of the 
target and independent of the primary voltage in the 
range studied. All three metals have two maxima below 
10 volts, the most pronounced one being that for Ag at 
3.9 volts. These maxima appear to correspond roughly to 
the regions of high optical absorption for the same sub- 
stance. An estimate of the depth of penetration of the 


jl is a well-known fact that when the molecules 
of a gas are bombarded by electrons of suffi- 
cient energy, there is a re-emission of electrons 
from the gas which can be regarded as composed 
of three different parts: Elastically reflected 
electrons, inelastically scattered electrons and 
low energy secondary electrons, resulting from 
ionizing collisions. Since the pioneer work of 
Franck and Hertz the energy distribution of the 
inelastically scattered electrons has been studied 
with the object of securing information about 
the energy states of the molecules in question, 
and the transitions between such states. In the 
same way the secondary emission from a bom- 
barded solid may be described as consisting of 
elastically reflected electrons, inelastically scat- 
tered electrons and secondary electrons proper. 
Some years ago it was shown! that the energy 
distribution of the inelastically scattered elec- 
trons exhibited a structure of maxima and 
minima, characteristic of the solid surface under 
bombardment. These occurred at a constant 
distance from the peak caused by reflected elec- 
trons on the energy scale, independent of the 
primary voltage. If the zero of the energy scale 
is determined by the position of this reflection 
peak, the curves evidently represent the dis- 


1 Brown and Whiddington, Leeds Phil. and Lit. Soc. 
Proc. 1, 162 (1927). 

2 Whiddington, Leeds Phil. and Lit. Soc. Proc. 1, 242 
(1928). 

3 Rudberg, K. Svenska Vet. Akad. Handl. 7, No. 1 
(1929). 


scattered electron was obtained from energy distribution 
measurements on films of Ca, CaO, Ba and BaO of known 
thickness, deposited on Ag. The results indicate, that 
the inelastic scattering is determined by the first few 
atomic layers near the surface. The fresh deposits of Ca 
and Ba show a very high rate of oxidation, even when the 
total pressure is <10~? mm; the distribution curves for 
the oxide are entirely different from those for the parent 
metal and exhibit a rather complicated structure. Measure 
ments of the energy distribution of scattered electrons 


should in some cases prove a sensitive tool for the study 
of surfaces. A deposit of 5 X10" atoms/cm? of Ca on CaO 
less than a monatomic layer in the Ca-crystal—is easily 


detected by this method. 


tribution of the energy losses suffered by the 
bombarding electrons. Such curves were ob- 
tained‘ from velocity distribution measurements 
by three different methods, employing focusing 
magnetic fields, a combination of magnetic and 
electric fields as well as focusing electrostatic 
fields. The targets studied were Cu, Ag, Au, Pt 
and the oxides MgO, CaO, SrO and BaO; some 
unpublished results have also been obtained for 
Ni. Later Haworth has reported similar results 
for well outgassed Mo. As in the case of a gas, 
the maxima in these curves were regarded as 
indicating that definite excitation processes were 
taking place in the bombarded substance, the 
position of a maximum representing the energy 
involved in that particular excitation. 

A possible way of accounting for these struc- 
tures theoretically was indicated in some work 
by Kronig and Penney’ and in the general theory 
of metals by Bloch ;’ certain general results which 
can be derived from Bloch’s work were used by 
Fréhlich® in connection with a special model in 
a theory of secondary emission. The theory of 
inelastic collisions in solids is discussed in a 
following paper? by Professor Slater and the 
present writer. 

4 Rudberg, Proc. Roy. Soc. A127, 111 (1930). 

5 Haworth, Phys. Rev. 48, 88 (1935); see also 37, 93 
(1931); 42, 906 (1932). 

® Kronig and Penney, Proc. Roy. Soc. A130, 499 (1931). 

7 Bloch, Zeits. f. Physik 52, 555 (1928). 


8 Frohlich, Ann. d. Physik 13, 229 (1932). 
® Slater and Rudberg, Phys. Rev., following article. 
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The early experimental work* was unsatis- 
factory because of the presence of vapor from a 
waxed joint in the tube; in order to keep the 
bombarded surface clean it was held at incan- 
descence during the measurements. It now ap- 
pears that this method, while successful in the 
case of the oxides, did not always give a well- 
defined surface in the case of metal targets. 
The present paper is a report on some work 
during the last few years, in which the distribu- 
tion curves for Cu, Ag and Au have been re- 
measured under more satisfactory vacuum con- 
ditions.'® In connection with the problem of the 
depth of penetration for inelastically scattered 
electrons, the scattering from films of Ca, Ba, 
CaO and BaO on Ag has also been studied. 


APPARATUS AND METHODS 


The essential features of the apparatus are 
indicated in Fig. 1. Most metal parts are made 
of sheet copper. A beam of electrons of nearly 
uniform speed, obtained from an electron gun G, 
is allowed to fall on the target T in question 
under an angle of incidence of about 45°. The 
target is surrounded by a cylindrical shield S to 
retain the bulk of the secondary emission. The 
total current to 7 and S thus measures the bom- 
barding current and is kept constant during a 
run. A narrow beam of the re-emitted electrons 
from JT escapes through a small opening in S, 
in a direction at right angles to the primary beam. 
Through a narrow slit it enters the deflection 
box B, where the energy distribution in the beam 


!© The material in this part of the paper was presented 
before the Washington Meeting of the Am. Phys. Soc., 
April 28, 1934. Phys. Rev. 45, 764 (1934). 
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is analyzed. The entrance slit and two other 
slits in the box define a circular path, which 
electrons of the proper speed are made to follow 
by means of a magnetic field, perpendicular to 
the path. Two extra partitions with slits serve 
to trap electrons scattered from the walls of the 
box. The electrons emerging from the third slit 
are caught by the collector C. With this arrange 
ment, the energy distribution can be measured 
in two different ways. In the first, the box is kept 
at the same potential as 7 and S, and the 
current 7 received by C is measured as a function 
of the magnetizing current J in the Helmholtz 
coils. The energy distribution is then given by 
i/1*? vs. I*, except for constant factors. In the 
second method, the magnetic field is kept con- 
stant at a suitable value, and 7 is measured as a 
function of a variable voltage v applied between 
the target and shield unit and the box. The dis- 
advantage with the second method, which yields 
results on an energy scale directly, is that the 
effective volume within which the electron tra- 
jectories must fall inside S in order that the 
electrons may reach C, varies both in shape and 
size with the applied voltage v, in a way which 
is difficult to correct for. This is most serious for 
low bombarding voltages V, when the range over 
which v is varied is not negligible compared with 
lV’. Hence the first method has been used except 
for a few measurements on Ag, where small 
parts of the energy distribution were studied with 
a variable accelerating voltage v. 

The primary voltage V, the primary current 
drawn by JT and S, as well as the magnetizing 
current in the second method, were kept sensibly 
constant during each run, the constancy being 


checked by compensating devices with galva- 


























Fic. 1. Apparatus. 
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first method, 7 was measured 


with a The 
was determined by compensation of the drop 
produced in a high resistance R, with an FP-54 
electrometer null instrument. 
The electrical connections are shown in Fig. 2. 


In the 
potentiometer. 


nometers. 


collector current 7 


tube serving as 


PREPARATION OF TARGET SURFACES 


All the targets studied were prepared by 
evaporation in vacuum directly in the tube. A 
small disk 7 of sheet molybdenum is attached to 
a magnetically operated slider, Fig. 1, one end of 
Thus 7 
can be withdrawn from the position for measure- 
S into the bulb to the left, where it 
can be outgassed by bombardment from the 
filament F or covered with a fresh evaporated 
layer of the metal to be studied. In the case of Cu, 
Ag and Au the metal is contained in a molyb- 
denum crucible 7, which can be heated by bom- 
Two such crucibles are 
this way 


which is resting on a track of glass rod R. 


ment inside 


bardment from F. pro- 
vided." The films obtained in 
sufficiently thick, so that no effect of the backing 
As a matter of fact 


were 


substance could be detected. 
the backing material soon became quite complex, 


Skinner and 
rard to 


1'[ am much indebted to Dr. H. W. B. 
Dr. H. M. O'Bryan for valuable suggestions in reg 


crucible design. 
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Fic. 3. Energy distribution of scattered electrons 


since a new surface was frequently formed by 
evaporation and surfaces of different metals were 
In one such a 


often studied alternatingly. case 


target with alternate layers of Cu and Ag was 
heated to incandescence, so that a reaction took 
place. The resulting surface upon cooling gave 
distribution curves quite different from those for 
Cu and for Ag. 

Ca and Ba were evaporated at accurately con- 
trolled temperatures molybdenum 
furnaces of a type described in connection with 


from small 
the measurement of the vapor pressures for these 
metals. > The rate at which the material was 
deposited could therefore be computed from the 
geometrical arrangement and the known vapor 
pressure. 

The apparatus was subjected to the customary 
outgassing by baking at 400°C high 
quency induction heating for the metal parts. It 
The pressure with the 


and fre- 


was pumped continuously. 


tube operating did not exceed the order of 
110-7 mm, as measured with an ionization 
gauge, calibrated for air. 
2? Rudberg, Phys. Rev. 46, 763 (1934). 
Rudberg and Lempert, J. Chem. Phys. 3, 627° (1935). 
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REsuLTs FOR Cu, Ag AND Au 


Typical distribution curves for Cu, Ag and 
Au are shown in Fig. 3,'for'a primary energy of 
about 150 volts. The shiatp peak to the left in 
each curve represents the elastic scattering. This 
serves to locate the zero of energy loss in the 
subsequent curves. The primaries themselves 
have a distribution, due to temperature and 
filament drop, about 0.5 volt wide (Halbwerts- 
breite). The width of the reflection peak is satis- 
factorily accounted for by this original distribu- 
tion and the finite resolving power of the analyzer 
which corresponds to a width AV = 0.008V on 
an energy scale. It is evident that the bulk of the 
emission is made up of low energy secondaries 
which exhibit the broad maximum to the right 


for low energies. It is possible that a closer 
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examination of this low energy part might reveal 
a structure such as has recently been reported 
by Haworth in the case of molybdenum after 
prolonged heat treatment.® Inelastic collisions 
are indicated by the part of the curve imme- 
diately to the right of the reflection peak. When 
the primary voltage is changed this part shifts 
bodily, whereas the low energy side of the curve 
remains essentially the same. In the following 
Figs. 4, 


elast ically scattered electrons. As may be inferred 


5 and 6 attention is confined to the in- 


from Fig. 3 the top of the reflection peak would 
fall well outside the figures on this scale; this 
peak, although not indicated, was always meas- 
ured in order to locate the zero position on the 
energy loss axis. 

For reasons which are not quite clear it appears 


to be more difficult to obtain well-reproducible 
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Fic. 5. 


curves with copper than with the other metals, 
but Fig. 4 gives a fair representation of the 
results. The curves are labeled by the serial 
number and the primary voltage, the latter in- 
creasing downward. The base metal is silver, 
molybdenum or gold. The first hump is always 
present, the second is present but varies in im- 
portance and in some curves, like the group at 
the bottom of the figure, there is a third hump, 
which in the earlier stages of evaporation actually 
rose above the second maximum. This may pos- 
sibly indicate that it is characteristic of some 
impurity distilling over with the metal. 

Curves for silver are shown in Fig. 5. They 
exhibit the sharpest maximum found so far at 


3.90 volts, and a broader maximum at about 7.8 


Silver. 


volts. The range of primary voltage is from 50 
to 370 volts. The broadening of the peak for 
higher voltages is of course due to the decreasing 
resolution of the velocity analyzer for higher 
energies. 

In the case of gold, Fig. 6, there are again two 
main maxima, the first one for an energy loss of 
3.05 volts. For this metal there is also a very 
broad hump around 24 volts. The relative pro- 
portion of higher energy losses is greater than 
for copper and silver. For high primary voltage, 
however, the gold curves bend down in this 
region like the others. The yield per primary, in 
the direction studied and for the voltage range 
employed, is less than for Cu and Ag by a factor 
of 2 to 4. 
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Typical curves for all three metals are collected 
in Fig. 7, which also illustrates the spacing of 
experimental points. Those for Cu and Ag have 
been shifted upwards by amounts indicated to 
the left. It will be seen that the curves for these 
three metals are similar in that there is a deep 
valley next to the reflection peak, which indi- 
cates the absence of inelastic collisions with smal! 
energy losses; this is followed by two main 
maxima. The positions of these in the case of Ag 
and Au are such that the second maximum could 
result from two of the processes giving rise to the 
first; however, the wide difference in shape, at 
least in the case of Ag, makes it very unlikely 
that this is the sole cause of the second maximum. 

The fact that the position of these maxima is 
independent (in the range investigated) of pri- 
mary voltage, of the potentials of other parts of 
the tube such as the electron gun, and of whether 
they are measured by a varied accelerating 
voltage or by a varied magnetic field—should 
prove conclusively, that these maxima are con- 
nected with characteristic transitions in the tar- 
get. A separate experiment showed that the sharp 
silver maximum could be reproduced within 40 


seconds after a fresh surface had been formed by 
distillation. At a pressure of 1X10-’? mm the 


. Gold. 


formation of a monomolecular layer by the 
residual gas requires at least that time if every 
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Fic. 7. Distribution of energy losses for Cu, Ag and Au. 
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Values of the energy loss for the first maximum in 
stlver and tn gold. 
Target Ag larget Au 


TABLE I. 


apenas ENERGY 
Ser. No. VOLTAGE volts) Ser. No. VOLTAGE Pe It 
93 50 3.95 130 127 2.95 
92 77 3.85 125 139 3.00 
63 109 3.65 148 140 2.95 
89 123 3.85 126 154 3.10 
62 141 3.75 133 154 3.15 
90 160 3.95 132 164 3.15 
212 175 3.85 127 165 3.05 
252 175 3.90 131 169 2.95 
87 188 3.80 122 183 2.90 
56* 194 3.80 144 185 3.10 
57 194 3.75 123 186 2.95 
64* 196 3.90 128 214 3.05 
65* 196 3.90 124 217 2.95 
66* 196 3.90 129 248 3.15 
67* 196 3.85 
68* 196 3.90 
61 222 3.85 
232 250 4.10 
173 250 4.00 
174 250 4.05 
60 267 3.80 
88 296 4.10 
59 320 3.80 
58 371 3.90 
mean of 24 values 3.88 mean of 14 values 3.03 
or 3.90 or 3.05 


probable error in mean +0.05 volt 


* Measured by method of variable accelerating voltage. 


hitting molecule sticks. Even granting that the 
pressure of residual gases during evaporation 
might have been somewhat higher than 10-7 
mm this result makes it difficult to ascribe the 
maximum to adsorbed gases. 

The values obtained for the position of the 
first maximum for silver and for gold are sum- 
marized in Table I. 


COMPARISON WITH OPTICAL DATA 


In a given system showing linear absorption, 
the number of transitions per second in unit 
frequency range at a frequency v may be written 
as T,E,? and the rate of absorption of energy as 
hvT,E,’, if E, is the electric vector of the radiation 
field. If the complex refractive index of the 
substance for this particular frequency is n—tk, 
the current density for the component in phase 
with the field is. ymkE,, so that the rate of ab- 
sorption of energy per unit volume is vnkE,?’. 
Hence 7, is proportional to the product nk. 
Kronig'* and others have shown, in the case of 


M4 Kronig, Proc. Roy. Soc. Al24, 409 (1929). 
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Fic. 8. Distribution of optical transitions (upper) and of 


electron energy losses (lower). 


the Bloch-Sommerfeld model of a metal, that 
there exists a selection rule for the electronic 
transitions which can be produced by absorption 
of radiation. According to this rule, which results 
from the space periodicity of the lattice and 
therefore holds for crystalline solids in general, 
transitions can only take place between states 
having the same (3-fold) momentum quantum 
number: ko’=Kp. The probability of such a 
transition is proportional to L*E,?, where L is 
the moment /Yyry'*dr for the two states y and 
y’ in question. Thus 7,=Const->°L*. It can be 
shown that, when the energy loss is very small 
compared with the energy of the bombarding 
electron, a similar selection rule should be 
approached for the case of excitation by electron 
impact, and that the probability of excitation is 
determined by L? also in this case.? Thus one 
would expect the curve for nk as a function of 
the energy hv to be similar to the distribution 
the same 


curve for small losses for 


substance, at least when the bombarding voltage 


energy 


is high. Fig. 8 shows a comparison of mk and the 
energy loss curves for Cu, Ag and Au, the 
frequency being measured in equivalent electron 
volts. The constants used are those given by 
Minor™® (Cu, Ag) and Meier'® (Au), in all cases 
measured in air. It is known that optical con- 
stants vary considerably with the surface condi- 


1903). 
1910). 


'§ Minor, Ann. d. Physik 10, 581 
6 Meier, Ann. d. Physik 31, 1017 
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tions, which probably accounts for the con- 
flicting results often found by different observers. 
It is possible, therefore, that the constants used 
could refer to a material quite different from 
that for which the distribution of energy losses 
has been measured. This is more likely to be the 
case with copper than with the other two; the 
nk-curve given for this metal in Fig. 8 may 
actually be the curve for cuprous oxide. The 
other two metals show a fair agreement as to the 
position of the main peak, although the optical 
maximum appears shifted to higher energies for 
Ag. The optical measurements by Pfestorf' show 
for Au a maximum at a somewhat higher energy 
than that given by Meier’s data, but his points 
are considerably more scattered. Pfestorf’s data 
for Cu give indication of the minimum in the 
nk-curve of Fig. 8, but do not show the maxi- 
mum. 


WITH FILMS CALCIUM 
AND THEIR OXIDES 


EXPERIMENTS 6) | AND 


BARIUM 


The theory of these collision processes, which 
is given in a following paper, is based on the 
model of the inside of a metal infinite in extent. 
In the experiments, however, electrons are shot 
into the metal (and escape from it to be meas- 
ured) through the surface. It is therefore con- 
ceivable that the collisions take place in a 
region at the surface, where the wave functions 
are quite different from those in the interior of 
the metal. No satisfactory theory is available 
for this region at present. 

In order to obtain some information about the 
depth of penetration, it was planned to cover a 
substance of well-defined energy loss character- 
istics, such as silver, with some other material, 
and study the changes in the energy loss curve 
as a function of the thickness of the latter. The 
cover substances chosen were calcium and 
barium, since the useful part of their vapor 
pressure curves appeared to fall in a range 
available for controlled evaporation, yet high 
enough so that one might expect the reflection 
from a clean, cold surface to be negligible. 
Special long exposure tests with the target 
withdrawn to such a position, that it could be 
reached by atoms reflected from the walls but 


'? Pfestorf, Ann. d. Physik 81, 906 (1926). 
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Fic. 9. Showing ‘“‘aging” of Ca-deposits. White circles: 
fresh deposit, black: 30 minutes old. Lower curves with 
continuous Ca-evaporation. 


not by the direct beam, actually indicated that 
such reflection was negligible. In the case of 
calcium, the target rested against the flat bottom 
of a reentrant tube (not shown in Fig. 1) during 
to the 
Some of the films were deposited while this tube 


exposure the atomic beam from gun. 


was cooled by COs-ice and acetone, others with 


the tube at room No difference 


could be detected. 


temperature. 


For convenience the amount of metal per unit 
area has been expressed as the number of single 
layers to which it would correspond, if arranged 
as in the planes of closest packing in that partic- 
ular metal, the (111)-planes for the face-centered 
Ca, the (110)-planes for Ba. It is not supposed 
that this the of the 
deposits, the structure of which is as unknown 


is actual 


arrangement 


as that of the underlying substance. 





146 ERIK R 














Jo 40 volts 









BaO 











Fic. 10. Distribution curves for ‘‘aged’’ Ca- and Ba 
deposits (white circles) compared with old data for CaO 


and BaO (black). 


It was found both in the case of calcium and 
barium, that a fresh deposit gave a distribution 
curve quite different from that obtained from 
the same target after a few hours standing in a 
high vacuum. The change was quite noticeable 
even in the time (about 30 min.) required for a 
run; the curve obtained from a deposit a few 
hours old could be reproduced a week later, if 
no further evaporation took place. This behavior 
was suspected to be due to a reaction between 
the freshly distilled Ca- or Ba-metal and some 
constituent of the residual gas in the tube, at a 
pressure =10~’ mm, the distribution curve after 
“aging’’ being due to a compound. The following 
facts substantiate this hypothesis and are be- 
lieved to prove that the compound formed is 


the oxide. 

Fig. 9 shows two pairs of curves for deposits 
of Ca. The white circles represent measurements 
immediately after the film was formed, starting 
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Fic. 11. Calcium on silver. 


with low energy losses. The black points were 
measured 30 min. later, also from left to right. 
Conditions in the case of the upper pair were the 
same as when the lower pair was measured, 
except that in the latter case the Ca-gun was 
left on in the evaporation chamber during the 
measurement. In this way the walls of the 
chamber were continuously being covered with 
a fresh coating of Ca, which acted as a getter 
removing the active gas from the tube. This 
evidently greatly increased the “‘lifetime’’ of the 
first high peak, which is the most characteristic 
feature of the pure calcium curve. Several years 
ago energy distribution curves for CaO and BaO 
were measured with a different apparatus.‘ The 
nitrates were deposited on platinum, by evapora- 
tion of the water from aqueous solutions, and 
these converted into the oxides by further heat- 
ing in air. Fig. 10 shows a comparison of the 
present measurements (white circles) on “aged” 
deposits of Ca and Ba and the old curves (black 
circles) for the corresponding oxides. The latter 
ones were measured according to the second 
method described (with an accelerating voltage 
v), which accounts for the gradual rise with 
increasing energy loss. Evidently every impor- 
tant feature of the old curves reappears in the 
present measurements on the ‘“‘aged’’ deposits. 
These will therefore be called oxides. 
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Fic. 12. Calcium oxide on silver. 


In the measurements on Ca- and Ba-films the 
gun was usually kept at evaporation temperature 
to slow down the rate of reaction, and the curves 
were taken as fast as possible, starting, as a 
rule, from the end of low energy losses, where the 
peaks are sharper. In this way the distribution 
for higher energy loss will exhibit, to some degree, 
features characteristic of the oxide. These can 
be distinguished from those of the metal itself 
by measuring each particular part of the curve 
separately, immediately after evaporation. In a 
new apparatus, now under construction, it will 
be possible to have the target sprayed con- 
tinuously with the metal during the measure- 
ments. 

Fig. 11 shows results for Ca-films at two 
voltages. The curves are labeled by their serial 
number and the thickness of the film in the unit 
defined above (thickness 1=7.4610" atoms 
cm*). The backing is freshly distilled Ag, except 
in the top curve in a, where the film is backed 
by an old oxidized Ca-deposit. The dotted curve 
No. 169 shows the result of “‘aging’’ on the film 
which yielded No. 168 immediately after evap- 
oration. Evidently the strongest oxide peaks are 
already beginning to appear in the curve for the 
fresh deposit. It will be seen that Ca has a strong 
peak very near the position of the sharp Ag- 
maximum and a broader hump in the region of 
the second maximum for silver. This makes it 
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Fic. 13. Barium on silver. 


impossible to state at what thickness the Ag- 
characteristics are obscured by the Ca-deposit, 
but it is evident that the film affects the distribu- 
tion curve already for a thickness less than 1. 

Energy loss curves for CaO on Ag are given 
in Fig. 12. The thickness values refer to the Ca- 
film before oxidation. In each case the deposit 
was allowed to “‘age’’ for a sufficiently long time 
before measurement, so that no further change 
in the curve would take place. The oxide maxima 
begin to appear for a thickness less than 0.5, but 
the sharp silver peak is not completely obscured 
by a thickness 10 times greater than this. 

The results for Ba-films on Ag are shown in 
Fig. 13 (thickness 1=5.6310" atoms/cm’). 
These films appear to oxidize even more rapidly 
than the Ca-deposits. Thus the BaQO-features 
appear in all the curves at higher energy losses. 
The curves were measured from left to right, so 
that the part for small energy losses, in which the 
characteristic features of Ba fall, is representative 
of Ba-metal. A thick deposit of Ba, or a film of Ba 
on BaO, shows no maximum near the silver peak, 
although thin films of Ba on Ag seem to enhance 
that peak. Ba has a low energy maximum at 
1.7 volts; this is so close to the reflection peak, 
that the maximum is only resolved at low pri- 
mary voltage. At 250 volts it is suppressed by 
overlapping from the strong reflection peak, 


except in the top curve, although its presence is 


indicated in the other curves by the trend in this 
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Fic. 14. Barium oxide on silver. 
TABLI ~ Measured energy losses corresponding to the top 
"Ag peak for Ca and Ba films on silver. 
ENERGY ENERGY 
THICK- SER. Loss rHick SER Loss 
NESS No. VOLTAGE (volts NESS No. VOLTAG v 
Ca film Ba film 
0.04 182 250 3.95 0.6 253 175 3.85 
0.04 183 250 4.10 0.8 213 175 3.60 
0.15 184 250 3.90 1.1 254 175 3.90 
0.36 186 250 3.85 2.2 255 175 3.85 
0.36 190 90 3.75 4 14 175 345 
9? 75 75 # 2 75 3.45 
—— a es Oe ie 
0.73 194 175 3.60 33 217 175 3.50 
0.73 198 90 3.60 3.5 238 250 3.75 
0.83 163 250 3.65 4.1 219 175 3.10 
1.3 164 250 3.50 4.4 256 175 3.60 
1.8 199 175 3.60 4.9 215 175 3.05 
2.2 166 250 3.65 4.9 220 175 3.20 
3.9 168 250 3.50 4.9 233 250 3.60 
4.4 170 250 3.55 4.9 236 250 3.40 
4.4 179 90 3.35 6.5 239 250 3.45 
4.4 180 175 3.55 6.5 259 175 3.40 
— : — 6.5 261 175 3.35 
mean of 10 values 3.56 6.8 217. 175 3.50 
or 3.55 | 8.7 257 175 3.55 
CaO film AS 260 175 3.40 
0.18 185 250 3.95 mean of 16 values 3.43 
0.36 187 250 3.80 | or 3.45 
0.36 188 175 3.75 cinabiaail 
0.36 189 90 3.60 BaO film 
— 4.9 216 175 3.05 
0.73 195 175 3.50 | 6.5 240 250 3.45 
0.73 196 250 3.60 8.7 258 175 3.35 
0.73 = 197 90 3.45 | 13 241 250 3.30 
1.8 = =200 0 «1753.50 | 13 265 1753.45 
18 = 201-175 3.50 | 23 242-250 «= 3.05 
2.2 165 250 3.40 | 46 243 250 = 3.40 
4.4 175 250 3.60 244 250 3.15 


mean of 7 values 
or 


3.51 
3.50 


mean of 8 8 values 3.28 


or 


R | 
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TABLE III. Energy losses for Ba targets. 
EN! y Loss 
Tat THICKNESS Ser. No \ AGI v 
Ba 4.1 219 175 1.80 
4.4 56 175 1.45 
4.9 215 175 1.75 
4.9 220 175 1.65 
on BaO 6.5 249 175 1.70 
6.5 261 175 1.65 
8.7 257 175 1.70 
on BaO 9.8 269 175 1.60 
13 260 175 1.75 
20 221 175 1.80 
49 263 175 1.45 
202 266 175 1.85 
202 267 175 1.50 
mean of 14 values 1.68 
or 1.70 
region. The dotted curve shows the “aging” 


effect on the deposit of No. 239 and the extent to 
which BaO-maxima appear in the Ba-curves for 
higher energy The Ag-peak shows up more 
strongly in No. 240 than in No. 239, indicating 
that the deposit becomes less opaque to the 
Definite Ba-char- 


lc SS. 


electrons when it is oxidized. 
acteristics appear at a thickness about 4, al- 
though of the Ag-peak can 
probably be noticed at a thickness of 2 

14 shows energy loss curves for BaO on 


the enhancement 


Fig. 
Ag. The oxide maxima appear for thicknesses 
less than the smallest one shown here. 

The measured energy losses corresponding to 
the top of the Ag-peak and the peaks for Ca 
and Ba are collected in Table II and III. 


DISCUSSION 


In the cases where sharp maxima do not 
overlap, CaO, Ba and BaO on Ag, 
indicate that the excitation process responsible 
for the Ag-peak still takes place when the 
exciting (and scattered) electron has to penetrate 
a layer several atoms deep. If the penetration 
is of the same order in other materials as in these 


the results 


three substances, the distribution of small energy 
losses can not be determined entirely by the 
surface atoms, but must depend on the structure 
of the solid below. In that case a theory of these 
processes, based on the electronic model for the 
interior of a large crystal, would be expected to 
give a reasonable approximation to the actual 
situation. 

It takes several layers of Ba before the 1.7-volt 
The maxima for CaO are present 


peak appears. 
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in their proper places for a thickness <1 and 
show no appreciable change (other than a gain 
in strength) as the thickness is increased. These 
maxima, however, occur for higher values of 
energy loss and they are rather broad, though 
distinct. Thus they are possibly chiefly deter- 
mined by the nature of the atom or molecule in 
question, and only to a lesser degree by the forces 
from neighboring elements in the structure; the 
shifts in energy due to changes in the latter are 
probably unimportant compared to the width of 
these maxima. The 
different in the case of Ca. For a thickness <1 


situation appears to be 
of this metal on CaO the sharp peak, which in 
Fig. 11 is superposed on the prominent Ag- 
maximum, shows up in the same position as for 
thick films. This may indicate, that a backing of 
CaO happens to act in just the same way on a 
single layer of Ca, as a backing of metallic Ca. 
If the distribution of energy losses results from 
the cooperation of several atom layers below the 
surface, one would in general expect a difference 
between the curve for a single layer and that for 
the thick metal. However, if we imagine a solid 
of N* atoms gradually built up from free atoms, 
so that first N atoms are brought together to 
form a string, then N of these strings lined up to 
make a layer and finally N layers piled up in 
the proper way—then the greatest change in the 
electronic energy levels would be expected to 
take place in the first step, and a relatively small 
change in the last. 

That the surface layer in general influences the 
position of the energy loss peak for the underlying 
substance is evident from Table II. In each case, 
all the values for films above a certain thickness 
(as indicated by a horizontal line) were included 
in the mean. The spread of the individual values 
is fairly wide, but the shift from the pure Ag- 
value of 3.90 volts is evidently much larger than 
the probable error in the mean, which is of the 
order +0.05 voit. Since at least part of the peak 
with Ca-films represents a characteristic energy 
loss for this substance, there is definite evidence 
of the effect in question only for films of CaO, 
Ba and BaO, where the shift is of the same 
magnitude, but increasing in this order. It is 
believed that this represents a general shift in 


the energy levels of the underlying silver target, 
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due to the forces exerted by the adsorbed film. 
This may be in the nature of a general change 
of the potential field in the target, or may 
result from a change in the interatomic distances, 
produced by the of the film. For 
“nearly free’’ electrons in a crystal the energy 


presence 


associated with a given set of quantum numbers, 
and therefore the energy difference for two states 
of given quantum numbers, is inversely pro- 
portional to the square of the lattice spacing. 
If the energy shift of 10-15 percent for the Ag- 
peak is to be accounted for in this latter way, 
it would indicate a linear expansion of 5-8 
percent. A calculation by Lennard-Jones and 
Dent'® for ionic crystals shows that the lattice 
spacing at the free surface should be smaller 
than that in the interior by an amount of the 
order of 5 percent. It is possible that a contrac- 
tion of the same order exists at the surface of a 
metal. When this is covered by some other 
substance, so that it is no longer free, the lattice 
at the surface may effectively approach an 
interior lattice, with a corresponding increase in 
spacing. 

In the preceding discussion it was assumed, 
that the evaporated deposits are spread uni- 
with the 
apparent surface area of the target. This must 


formly over an area comparable 
certainly be true at the moment the deposits are 
formed, since there is no reflection. It is con- 
ceivable that patches of high concentration, 
separated by large bare regions, are subsequently 
formed by the process of surface migration. 
While this can 
appear likely in view of some of the experimental 


not be disproved, it does not 


results. Thus it is hard to understand the high 
yield of CaQ-collisions with deposits <1, if the 
greater part of the surface is bare. The rate of 
migration over a bare surface can not be high at 
the temperature used. In some experiments the 
back of the target was exposed to the metal gun, 
so that a deposit several units thick should have 
The 


tested always gave the curve characteristic of 


been formed. front side, however, when 


Ag, showing the absence of migration round the 
edge. On the other hand the apparent electron 


18 Lennard-Jones and Dent, Proc. Roy. Soc. A121, 247 
(1928). 
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transmission through thick BaO-deposits (Table 
Il) may indicate the presence of thin regions, 
formed in the process of oxidation of thick Ba- 
films. That these regions are not bare is evidenced 
by the shift of the Ag-peak. 

The distribution of energy losses appears to 
offer a new, sensitive tool for the study of sur- 
faces and surface layers. If the substance of the 
film has a distribution curve reasonably different 
from that of the backing material, layers only a 
few atoms thick This test 
appears to be somewhat more sensitive, and less 
likely to alter the studied, than 
present electron diffraction methods employing 


can be detected. 


conditions 


15, 1936 
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high voltages.’ In the case of Ca, the metal can 
be detected in a thickness decidedly <1 on a 
backing of the oxide. This is now being used to 
study the surfaces of thermionic emitters made 
from CaO and BaO. For the detection of small 
amounts of foreign material on a surface this 
method is probably never as sensitive as thermi- 
onic measurements, but in favorable cases it has 
one great advantage compared with these: it not 
only detects, it identifies. 

’ The diffraction beams obtained with single crystals 


at low bombarding voltages appear to be even more sensi- 
tive than the energy loss distribution to a small coverage 


of foreign material on the surface: Farnsworth, Phys. Rev. 
49, 605 (1936). 
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The general properties of the single electron wave func- 
tions for a crystal are discussed, and an expression derived 
for the probability of excitation by electron bombardment 
as a function of the two electronic levels involved in the 
transition. When the excitation energy is small compared 
with the energy of the bombarding electron, the relative 
transitions probabilities approach those for excitation by 
radiation. A detailed calculation has been carried out for 
the case of solid copper, for which the energy levels are 
now known from other work. Approximate wave functions 


of the atomic type for the lower state, free-electron-like 


HE problem of inelastic scattering of elec- 

trons in matter, early discussed for relatively 
fast electrons by Thomson! and others as a 
problem of atomic scattering, was treated very 
completely by Bohr? on the basis of classical 
electron theory. More recently, quantum theory 
has been applied particularly to the scattering by 
free atoms. Fréhlich* has given a wave-mechan- 
ical treatment of the inelastic scattering from 
metals, using a special Bloch-Sommerfeld model 
for the metal. He was able to deduce from this 
the characteristic features exhibited by the dis- 


1 J. J. Thomson, Conduction of Electricity through Gases, 
343 (1903). 
2N. Bohr, Phil. Mag. 25, 10 (1913); 30, 581 (1915). 
H. Frohlich, Ann. d. Physik 13, 229 


1932). 


for the upper state—were used to evaluate the transition 


integrals. The calculated distribution curve for the in- 
elastically scattered electrons is in fair agreement with the 
experimental curve in the region of small energy losses: the 
rapid rise from small values and the following two maxima 
exhibited by the experimental curve are all reproduced in 
their proper positions. The departures for higher losses 
are attributed to the fact that the free-electron-like wave 
functions used are not a proper approximation to the actual 


functions in this region. 


tribution of the low energy secondary electrons. 
In the present paper we consider the transition 
probabilities for electronic excitation on the basis 
of the general electron theory of solids. In experi- 
ments on inelastic collisions, the target is usually 
bombarded by a nearly parallel, homogeneous 
beam of electrons and the energy distribution is 
measured for those electrons which are scattered 
into a small solid angle about a definite fixed 
direction with respect to the target and the in- 
cident beam. Tests with polycrystalline targets 
have shown the energy distribution to be essen- 
tially independent of both this latter direction 
and of the angle of incidence. For this reason we 
the total probability of all 


have calculated 
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transitions with a given energy loss and arbitrary 
direction of the scattered electron, assuming an 
isotropic distribution in direction of the primary 
particles. It seems likely that the isotropic dis- 
tribution of primary and scattered electrons here 
involved is actually nearly established inside the 
target as a result of the frequent processes of 
elastic scattering by the massive ions in the 
lattice. 

While most experimental work in this field, for 
technical reasons, has been done on metallic 
conductors, the theory applies to solids in general. 
The present discussion divides into three main 
parts. In the first the wave functions and their 
properties are reviewed ; in the second the general 
expressions for the transition probabilities are 
developed. In the third part approximate cal- 
culations are carried out for the case of copper, 
for which the energy levels are now known and 
approximate wave functions available, and the 
results are compared with experimental data. 


- 


We wish to consider a collision in which a fast 
and a slow electron collide, the fast one losing a 
little energy and the slow one gaining it. The 
unperturbed problem is that in which each 
electron moves through a_ periodic potential 
field, representing the crystal, unaffected by the 
particular motions of other electrons. Each 
electron in this approximation is represented by 
a modulated periodic function of the type well 
known in the theory of metallic lattices : 


Yaro(r) =exp [iko-r |uar(r), (1) 


vat, having the periodicity of the potential field. 
We shall often use the equivalent Fourier reso- 
lution : 

Yoko) = >-a(ky+K,) exp [7(Ko+K,)-r]. (1a) 

- 

The stationary states are divided into a number 
of bands, and the index a represents a particular 
band. The vector Ky is different for the various 
states within a band. To satisfy Bloch’s condi- 
tion of periodicity,‘ the vectors K, must form a 
discrete set, constituting what is called the 
reciprocal lattice in the k-space. This reciprocal 
lattice is determined by the crystal structure, 


‘F. Bloch, Zeits. f. Physik 52, 555 (1928). 
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and is face-centered cubic for a body-centered 
cubic crystal, and vice versa. We can conveniently 
set up polyhedral cells in the k-space, one sur- 
rounding each point of the reciprocal lattice, and 
of such a shape that they fill the k-space com- 
pletely without overlapping. We may imagine 
a(Kyo+K,) =a(k) as a single-valued function of 
the continuous variable k in the k-space; then 
the Fourier coefficients of a particular state Ko 
of the band a will be the values of a at the ap- 
propriate points k of the various cells, points 
reached from Kp by the translations K,. It can 
be shown that the square of the magnitude of 
one of these a’s is proportional to the probability 
that the electron in question have a momentum 
(h, 27)k=(h/2r)(Ko+K,), when in the state 
(a, Ko). For this reason the a’s are sometimes 
called momentum eigenfunctions. It is evident 
that Yer is a periodic function of ko; for if Ko is 
increased by one of the vectors of the reciprocal 
lattice, say K,, the whole lattice of points ky) +K, 
will be shifted back into itself. The same Fourier 
coefficients must then be chosen, leading to the 
same function as before, since by arranging the 
states according to bands we have made y a 
single-valued function of kp within each band. 
At the same time the energy of this state will 
return to its initial value, so that the energy 
within a band may be regarded as a periodic 
function of the continuous variable Ky (or k). 
Often, on account of this periodicity, it is con- 
venient to restrict ky to the central polyhedron 
of k space, but this is not always desirable. 

If in particular the state a represents a free 
electron, its wave function is a plane wave, a 
constant times exp [7k,-r] (where k, need not 
be in the central polyhedron in k space), so that 
only one a is different from zero. The energy of 
the free electron is p*?/2m, where p is its mo- 
mentum. If we express energy in Rydberg units, 
distances in atomic units, this becomes simply 
2°, where k; is proportional to the momentum. 
For the sake of uniformity, it is often convenient 
to represent the state of a free electron by a 
function in k space, and to make this single- 
valued one must break up the free electron states 
arbitrarily into bands. We could do this by 
deciding to group together all those states for 
which k; was in a particular polyhedron in 
k space. For this particular band, a(ky+K,), 
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regarded as a function a(k) of the continuous 
variable k=k)+K,, would be zero outside this 
polyhedron, constant within it. The energy, 
within the polyhedron, would be given by &,’, 
and by the periodicity property this would be 
repeated in every other polyhedron. This par- 
ticular assignment, however, is not very reason- 
able when there is a small periodic perturbative 
potential. In this case, the wave function and 
energy will be slightly perturbed, and it can be 
shown that the periodic energy function which 
we have just discussed would be distorted by the 
introduction of certain planes of discontinuity, 
across which the energy would change discon- 
tinuously. Brillouin has shown’ that it is much 
better to take unperturbed functions by setting 
up certain zones in k space. Such a zone is formed 
from parts of many polyhedrons, all having 
about the same energy and |k,| value, but 
spread more or less uniformly about the surface 
of a sphere. These parts of polyhedrons are so 
chosen that they could be fitted together, by pure 
translations K,, to form a single polyhedron. We 
now assume a single energy band of the free 
electron to be one in which a(k) is constant 
within one Brillouin zone, zero outside it. The 
energy within the Brillouin zone is given by k,, 
and outside the zone it is a periodic function of 
k, coming back to the same value at correspond- 
ing points of each polyhedral cell. The Brillouin 
zone is chosen in such a way that this determines 
the energy at each point of k space uniquely. 
The energy, however, is a rather complicated 
function of k within an individual polyhedron. 
The polyhedron, as we have said, can be divided 
into parts which could be displaced by K, to 
form the zone, and within each of these parts the 
energy is represented by |k+K,/?, a parabolic 
function of k; but this function changes from 
one part of the polyhedron to another, since the 
proper translation vector K, is different for dif- 
ferent parts. Because of this, the slope or gradient 
of the energy as a function of k changes discon- 
tinuously from one part of the polyhedron to 
another, though the function itself is continuous. 
This is because the direction of the vector k, for 
the free electron, but not its magnitude or k,, 
changes discontinuously with k in going from 


5 L. Brillouin, Quantenstatistik (1931). 


AND 5. € 





SLATER 


one part to another. When small perturbations 
appear, however, the discontinuous changes of 
slope are rounded off, and the energy function is 
a smooth one. It is this property that leads one 
to say that the Brillouin zone is the correct zone 
to use for describing the free electrons. The 
resulting type of function is the sort actually 
found in the wave functions of fast electrons in 
real lattices. 

We shall have occasion to consider the ortho- 
gonality and normalization of our wave func- 
tions. The vector k as a matter of fact cannot 
take on arbitrary values in the k space, but only 
a discrete set lying on a sublattice much more 
closely spaced than the reciprocal lattice, such 
that the number of points within a single poly- 
hedron of the latter equals the number of atoms 
in the crystal. Each state is capable of holding 
two electrons, one of each spin, so that each 
band can accommodate two electrons for every 
atom of the crystal. This sublattice has the 
property that, if k and k’ are any two points of 
it, fexp [7(k—k’)-r jdr=0, if the integration is 
carried over the whole crystal, unless k=k’, in 
which case the integral of course is Vo, the volume 
of the crystal. Now consider two wave functions, 
Yaro and Yor. The normalization and orthogo- 


nality conditions demand that 


[ EXereo+K,)ode’ +K,) 
1, ifa=b, ky=k,’ 
Xexp [—7(Ko+K,—k,’—K,)-r]dr= 


0, otherwise. 


If ko#ky’, the property of the sublattice men- 
tioned above immediately shows that the integral 
is zero, so that the only remaining case is that in 
which ky=ky’. This yields the conditions 
0 if axb 

> a*(ko+K,,)b(ko+K,,) = wm (2) 

P 1 ] oil a=b. 
For a plane wave, when only one a is different 
from zero, we at once see that the wave function 


is Vy-' exp [7k,-r }. 


2. 


Returning now to our problem, we assume the 
fast electron both before and after collision to 
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be moving fast enough so that it can be repre- 
sented by a plane wave. If it has values K and K’ 
before and after collision, its wave functions will 
be Vo" exp [7K-r] and Vo" exp [7K’-r]}. It is 
immaterial whether these happen to belong to 
the same Brillouin band or not. We assume the 
slow electron to have wave functions yYaro and 
Ya'k’ before and after collision, described by the 
Fourier components a(Ko+K,), a’(Ko’+K,), re- 
spectively. Let the total energy of the two elec- 
trons before collision be W, afterwards W’, where 
W=W’ if conservation of energy holds. The 
unperturbed wave function of the two electrons 
together will be the product of the separate 
functions. Thus before collision, the function, 
including its time variation, is 


1 
en2riWey exp (iK-re lW«, ko(Y1), (3) 
V3 


where Tr; is the coordinate of the slow electron, 
re that of the fast one. We treat the collision by 
the perturbation method of variation of con- 
stants. At time t=0 we assume that the Coulomb 
interaction e?/|re—r;| commences to act. The 
solution of the time-dependent Schrédinger 
equation can then be expanded as a sum of terms 
> cla’ky’K’t)e2"'t! 

x... 


: exp [7K’-re ly ‘(r,). (4) 


If we assume that at ‘=0 the only one of the 
coefficients ¢ different from zero is the one mul- 
tiplying the function (3), and that this is unity, 
perturbation theory leads to the value 


c(a’k’ K’t) =w- "(1 —e2 tit") 7, (5) 


where w=W’—W’, the difference in energy 
between the final and initial states, and 


1 a e 
™ ah exp [i(K—K’)-r2] 
Vo ro—T 


x Wako( 11) W*a'ko’(¥1)d 71d 72. (0) 


The square of the modulus of c(a’ko’K’t) measures 
the probability of finding the system in the state 
(a’ko’K’) at time ¢; for large values of ¢ this will 
have appreciable values only near w=0, repre- 
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senting conservation of energy. (0/0/)|c|* meas- 
ures the transition probability. 

Let us consider first the interaction integral 
(6). Following Bethe,® we perform the integration 
first over the coordinates of the fast particle. 
Letting AK=K-—K’, r=r.-—r, and AK-r 
=AKr cos 6, we have 


exp [(AK-rp | 
f dr2=exp [iAK-r, 


fo—T, 
x | r—teiak 92rr*? sin 6d0dr 
=exp [/AK-r, |-- sin (AKr)dr. 
AKe 


The last integral does not converge. To avoid 
this difficulty, Bethe averages over a small range 
in AK, and after integration makes this range 
approach zero. Another method is to replace the 
exact Coulomb potential 1/r by e~*"/r, perform 
the integration, and then let a—0. Both methods 
give for the integral the value 1/AK. Hence 


4re? 
[= [voces exp [iAK-r, lW*arko'(¥i)d 7; 
(AK)? Vo 
4re? Z 
7 > a(ky+K,,)a’*(ko’ +K,) 
(AK)?Vo ”2 


x fexp Lido +K, -ky’—K,+AK)-1, ]d7:. (7) 


As mentioned above, the integral is zero unless 
k,)+K,—k,’ —K,+AK=0, or ko’+K,—(ko+K,) 
= AK, in which case it is Vo. If this condition is 


satisfied, we have 


4re? 
[= > a(ko’+K,—AK)a’*(ky’'+K,) (8) 
(AK)? ? 
or 
4re* 1 > 
I= Woro(¥}) exp [1AK-r; | 
(AK)? Vor vo 


4re? 1 ¢ 
X b*a'ke'(81)d71= erty | Waro(T) 
(AK )? Vor vo 


Xexp (7AK-r]y*eto(r)dr, (8a) 
where vo is the volume of one lattice cell of the 


6 H. Bethe, Ann. d. Physik 5, 325 (1930). 
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Fic. 1. Illustrating the possible transitions in the K space 
of the fast electron. 


crystal and the last integration is over such a cell. 
The last step depends on the space periodicity 
of the wave functions (1). In (8) the quantity /, 
whose square determines the transition prob- 
ability, is written as a sum of terms, of which 
each one the fast 
electron from momentum K to a value K— AK, 
and an equal and opposite transition of the slow 
electron to momentum k’=k)’+K, from a value 
k’—AK. This term has in the numerator the 
product of the amplitudes of the probabilities of 
finding the slow electron initially with momen- 
tum k’—AK, when in the state (a, ky’ — AK), 
finally with momentum k’, when in the state 
(a’, Ko’), and in the denominator the square of 


represents a transition of 


the momentum change AK. 

We are now ready to return to expression (5), 
giving the rate of transition from the initial to 
the final state. We really are not considering a 
single transition, however, but all possible transi- 
tions for which the fast electron originally has 
approximately a definite energy, but arbitrary 
direction, and the bound electron changes from 
state (a Kp) to state (a’ko’), gaining an amount of 
energy €; we wish to compute (0/dt)}.|c/|? for 
all these transitions. Let Fig. 1 represent a plane 
section through k space (the reciprocal lattice and 
the sublattice are not shown) parallel to the vector 
AK, the z axis lying along the direction of AK. 
The points K where the incident electrons have 


nonvanishing Fourier components are assumed 
to lie in a spherical shell of thickness 6K and 
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radius K. This shell will not in general be entirely 
in one Brillouin zone, but will include parts of 
several zones. It contains a stationary state for 
each point of the sublattice, so that it is popu- 
lated to a uniform density p. These incident 
electrons all have approximately the same energy 
K?. A second shell, inside the first, includes all 
vectors K’ for such states that the net energy 
increase of the system is between w and w+dw. 
Thus K”’+e=K?+w, and the thickness of the 
shell is dw/2K’. From the figure it will be seen 
that 
which can give rise to transitions in the range dw 
is dn=pdv, where dv is the volume generated by 
the small shaded area upon rotation around the 


the number of primary electron states 


z axis. Each of these states gives rise to one single 
transition. With the notations of Fig. 1, 


dv=27rK sin giK Kdg=2rKiKd(—2K cos ¢) 

and from the triangle 
K?+AK?*—2AKK cos ¢=K” 

so that 

d(—2K cos ¢)=(1/AK)d(K”) =(1/AK)dw 
and dn= prK6K(1/AK)dw. (9) 
The probability of all the specified transitions 
with a given AK is then 
0 11 —ertivt/h? or Kb Kdw 

In| 
at| ‘ w AK 


—_ ae 





0 | prKok mt} A 


= fF b = J*(Ko’, AK), (10) 
dtl AK h} AK 
where A =167°e!pK6K/h and 
1e¢ 
J(k,’, AK) =— | Yor(r) exp [tAK-r ] 
o’ Vo 
XW*arko'(r)dr (11) 
or 
J(Ko’, AK) = }-a(k,y’+K,—AK)a’*(k,’+K,). 
p (11a) 


This holds provided this type of transition is 
possible at all. From Fig. 1 it can be seen that 
transitions are possible only if 


2K —¢«/2K >AK >¢/2K =AK~,,. (12) 


In the experiments on inelastic scattering from 
solids, the lower limit AK,, in (12) is usually 
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quite small compared with the spacing of the 
reciprocal lattice. For small values of AK the 
exponential factor in (11) can be expanded and 
higher terms than the second neglected, giving 


1 . 


I(ko’, AK)=— | yu o(r){1+iAK-eldz 


if *,, + 
(roy A 


; 
= AKL COS (AK, RE (13) 
Vo 


the first term vanishing on account of orthogo- 
nality. Here 


Yaro(r) ry *arko'(r)dr, (14) 


the linear moment for this pair of states. 

Given the upper state (a’, Ko’) there will be a 
number of lower states a from which transitions 
to this state can take place with a change of 
AK and 


within the 


momentum of magnitude between 
AK+dAK AK) 
solid angle element dw. This number is propor- 
tional to AK*dAKdw; has a 
probability (10). The total probability of transi- 
ky’) will therefore be pro- 


and direction (of 


each transition 


tion to the state (a’, 
portional to 


J*(Ko’, AK) 
J d(AK )dw, 
“ AK? 


the limits of integration being determined by the 


(15) 


condition (12). 

It turns out that, when small values of «/2K 
are possible, the greater part of all transitions to 
the state (a’, ky’) comes from those terms in the 
probability sum (15) for which AK is very near 
its minimum value AK,,. In that region we can 
use the expression (13) for J. We should there- 
fore obtain a fairly good approximation to the 
actual distribution of small energy losses of the 
scattered primary electrons, by allowing each 
k,’ of an unoccupied state to receive transitions 
only from the occupied lowersstate with ko = ky’, 
assigning to each such!trawsition a probability 


proportional to p te 
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Fic. 2. Distribution in energy of the higher electronic levels 
in solid copper 


eit L? cos? (AK, L) 
d(AK)dw 
AK 


1L?-In (2K /e) (16 
or proportional to L? for very high K where « 
the energy difference between the two states. 
This approximation must be restricted to the 
region of small energy losses. When e is large, 
the approximation would introduce a consider- 
able error, except for very fast primary electrons, 
because AK,, would not be small. The effect in 
the region of small losses is essentially to make 
the peaks and valleys in the calculated energy 
distribution appear sharper than they should be, 
because of the sharp selection rule by which we 
have replaced the distribution of transition 
probabilities in the sum (15). For high AK and 
small energy losses, where the selection is quite 
strong, the transition probabilities for electronic 
excitation become practically proportional to 
those for optical absorption, both being deter- 
mined by the square of the moment L. 


3. 


We shall apply the method outlined in the 
preceding section to the calculation of the energy 
scattering from solid copper. The 
this case are from the 


losses on 
transitions involved in 
partly occupied band deriving from the 3d of 
the free atom to the unoccupied states above it, 
of which the nearest ones correspond to the 
atomic 4s. Fig. 2 shows the distribution in energy 
of the electronic states for solid copper; the 
area the part normally 


shaded represents 


7 From an extension of the work by H. Krutter, Phys. 
Rev. 48, 664 (1935); see also J. C. Slater, Phys. Rev. 49, 
537 (1936). 
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occupied by electrons. Since the energy levels 
are known for the different states of all the bands 
in question, the energy required for each par- 
ticular ‘‘vertical” transition, Ky=Ko’, can be 
found. To calculate Z from (14), however, and 
hence the weight L? to be attached to each such 
transition, we have to find approximate wave 
functions for the lower and upper states involved. 

The lower states a constitute a band of essen- 
tially bound electrons. If the solid is imagined as 
formed by bringing the free atoms together, the 
overlapping of neighboring atomic 3d functions 
will therefore be small for this band and Bloch’s 
approximation to the wave functions can be used. 
From the general properties of the Bloch function 
it can readily be shown that its introduction to 
represent the lower state in (11) and (14) results 
in replacing Yoto in these equations by (Vo/v0) Wa, 
where ¥, is the corresponding atomic wave 
function, in this case 3d. Thus for a given ky’, 
L is independent of AK, and L?=L?. The integral 
over the finite volume Vo of the entire crystal 
can then be replaced by an integral over all 
space, since the atomic function y,(r) has appre- 
ciable values only in a very small region of r. 

Actually the 3d band is 5-fold, as is the atomic 
d state, and accommodates 10 electrons per atom 
like the atomic state. This multiplicity arises 
from the five possible m values, 2, 1,0, —1, —2 
for a d state; in general, however, a single 
member of the band does not correspond to one 
particular m, but to a linear combination of all 
five. We shall only need the average L for all 
five components. To a good approximation the 
wave function for the mth component state of 
the 3d in copper may be written’ as the product 
of a radial function R(r) and the tesseral har- 
monic 7;"(¢, @) in the polar angles of r for /=2, 
each normalized to unity for the entire range of 
the variables : 


(-—|m|)!)3 
v¥.=R(r)- —+ emeP,|™ (cos @) (17) 
(1+ |m)\)! 
with 
R(r) =r°(75.8e—*-79" +-4.236e-2-85" 
+0.1185e—!:15"). (18) 


The unoccupied upper states behave much like 


~ 8 J.C. Slater, Phys. Rev. 42, 33 (1932). 
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states of free electrons, at least those immedi- 
ately above the 3d band, which may be said to 
derive from the atomic 4s. One might therefore 
try to represent such a state by a free electron 
wave function V)~! exp [7k,-r], where k; could 
be determined by the condition that k;? should 
equal the of the state in question, 
measured from a suitable zero, fixed by the loca- 


energy 


tion (in energy) of the lowest 4s level. This 
choice of function, however, does not satisfy the 
condition of orthogonality, and therefore makes 
(13) invalid; for the same reason it does not 
make J vanish for AK =0, as the proper function 
should according to (11). To represent the upper 
state (a’, Ko’) one would therefore have to set 
up a function yx, which is orthogonal to the 
atomic function y, and yet which resembles the 
plane wave V,~! exp [7k,-r]. Such a function can 
be found by using the expansion for a plane 
wave in a series of spherical harmonics, corre- 
sponding to s, p, d, f, --- components, and by 
assuming that the function y;, equals the plane 
wave, without the d spherical harmonic : 


1 


1 /2r\? 
pei(r) = _( ) > (+4)0'Jii(kir) 

Vo? \Rir I=0 

l2 
<Pif(cos@). (19) 
Here r, 6 are spherical polar (zonal) coordinates, 
the polar axis being the direction of k,. Since the 
other harmonics are orthogonal to the atomic d 
function (17), this removes the difficulty of 
nonorthogonality. As far as the calculation of L 
is concerned, it really makes no difference whether 
we use the simple plane wave or the expression 
(19), since the contribution of the omitted d term 
to the moment is zero. Of course (19) is not quite 
correctly normalized for a crystal of finite size, 
because of the omission of the d term; we have 
ignored this. The resulting function resembles a 
plane wave ineifhe. region where the atomic 
function ssizappfediable, but is distorted par- 
ticularly fordgieiysmall values of r where the d 
spherical haitmpniedsifelatively large. This is the 
type of behiavigs agtually shown by the wave 
function for ene@f thehigher states in the metal 
in the neighborhood: of rhe lattice points. In 
k space the simple plane wave has a single com- 
ponent for its Fourier resolution at k=k,. 
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Actually, with the d term missing, the Fourier 
resolution of (19) will still have a large com- 
ponent at the point k,, but in addition there 
will be fairly small components at all possible 
points k,+K, in k space. This in general is true 
of the correct wave functions of electrons in 
crystals, since they are never exactly plane 
waves. Of course our approximation is not very 
good. The correct function can be represented 
as a sum of s, p, d, components, which are 
not very well approximated by the Bessel’s 
functions, and the coefficients multiplying them 
are not well given by the coefficients in the plane 
wave expansion. Nevertheless, the general fea- 
tures are fairly well given by our assumption, in 
particular the way in which the wave-length of 
the oscillations of the function changes with 
energy, and this is the determining feature in 
the’ variation of L? for different upper states, 
which is the fundamental thing we shall get 
from our calculations. In the approximation 
adopted here the atomic d state is degenerate, 
the five m values all corresponding to the same 
energy. We are therefore free any 
orientation of the coordinate system in (17) ; we 
shall take a polar axis that coincides with the 
axis used in (19), the direction of k,. With this 
choice, the vector L is most conveniently deter- 
of rcos @, 


to choose 


mined from the matrix elements 
r sin 6e'?, 
components of the vector r, rather than x, y, z. 
The only terms in (19) which contribute any- 
thing to the moment are the p and f terms for 
l=1, 3, according to the same selection rules 
which apply for atoms. Furthermore, only the 
functions (17) for m=+1, 0 contribute any- 
thing, again in agreement with the selection rule 
for atomic transitions. For a case m=0, the 
only component of L different from zero is that 
of rcos 6; for m=1, that of rsin 6e**, and for 
m=-—1, rsin 6e~'*. The components have the 


following values : 


rsin 6e~‘*, representing the three 


m=0, 
L, cos 0 = Vo *((4 01/5) M+ (6201/5) My), 
m=-+1, 
Ly sin be*'? = Vo“*(2riy 6/5)(M,+M,y,), (20) 


where 


LE ( 
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Fic. 3. The average square moment L? as a function of 


the energy &,? for the upper state. 


i 


™ av 3 
M,= | PR()( ) Js;o(kir)dr 
e kur 


| 1—0.01516k,? 
=k} 0.0933 


+2.218 
(1+0.0253k,2)5 
1 —0).0856k;? 1 — 0.4545 2,7 
x +21.35 
(1+0.1425k,?) (1+0.757k,*)5 
and (21) 
” 2r \! 
My=f PR()( ) J; o(kyr)dr 
0 kyr 
0.00378 0.5050 


=k + 
(1+0.0253k;7)5 (1+0.1425k,)5 


25.85 
+ . -Fe 
(1+0.757k,2)5) 


To obtain L? we square the L value for each m, 
add and divide by 5, since all five m values occur 
with the same weight in the band of d functions. 
The result is that 


(22) 


L?=Const. (M,?+3M,M,+(21/16)M). 
In Fig. 3, L? is shown as a function of the energy 
k of the upper state. 

To construct the curve to be 


energy loss 


expected on this theory, we now proceed as 
follows. The energy range in question is divided 
into elements of the same size, Ae—in this case 
Ae=0.01 atomic unit. The transition probability 
L? is evaluated for each one of a series of upper 
states, picked with a uniform spacing as regards 
ky’. The energy involved in the corresponding 
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Fic, 4. Comparison of calculated and experimental dis- 
tribution curves for the energy losses in solid copper. 


vertical transition is determined from the 
tabulated energy levels for the different bands, 
and the L? value listed in the column for the cor- 
responding element Ae. The distribution curve 
for the energy losses is obtained by adding the 
L? contributions in each such column. The result 
is shown in Fig. 4 together with the experimental 
curve,’ which is given for comparison. The agree- 
ment in general is not unsatisfactory. Both curves 
start with low values for zero energy loss. In the 
theoretical this from the small 
number of possible transitions with small energy 
losses. The curves then rise rapidly, each one 
coming to a peak at about 7 volts, preceded by a 
subsidiary peak at about 4.3 volts. At higher 
voltages both curves fall off. The agreement here, 
however, is less satisfactory, since the experi- 
mental curve while the 
theoretical one shows marked further peaks. In 
regard to these differences one may make the 
following comments. Our calculation of the 
number of transitions for a given energy differ- 
ence has been based on actual calculations of 


case arises 


falls off smoothly, 


energy levels, and should be fairly reliable. It 
takes account of the distribution of density 
within the band of d electrons and of the devi- 


*E. Rudberg, Phys. Rev., preceding paper. 
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ation of the excited electrons from free electrons. 
It is this part of the calculations in which the 
peaks are introduced. Those at 4.7 and 7 volts 
arise in a rather fundamental way: the band of 
d electrons shows specially large concentration 
of states at its twe extreme energies, about 3 
volts apart, and it is roughly these two concen- 
trations which result in the two peaks. The two 
maxima at 14 and 18 volts also arise from these 
two concentrations, but here combined with 
peculiarities in the excited electron bands, about 
which there is much more uncertainty, on the 
basis of independent evidence. These peculiarities 
are connected with an especial concentration of 
excited ¢ components in the stationary states 
for these bands. But in calculating the transition 
probability L* we have represented the upper 
state by a free electron function without the d 
component. If, as it happens in this region, a 
large part of the correct wave function is actually 
d, our method must give too big values for L?’, 
since d—-d transitions should not occur. One thus 
infers that the proper transition probability to 
be used in this case should be small, which would 
tend to counteract the high concentration of 
transitions, thus reducing the violence of the 
fluctuations in this range. No such feature is 
present with the peaks at 4.7 and 7 volts. Of 
course our approximation to consider the transi- 
tions as vertical becomes increasingly poor for 
higher energy losses. Apart from this, one must 
bear in mind that there is no direct way of dis- 
tinguishing between the effects produced by 
single and multiple impacts in the experimental 
curves. A fast electron may have two or more 
inelastic impacts before escaping from the metal. 
Since it tends to lose only a few volts on one 
impact, it would then emerge with something 
like double the energy loss in a single collision, 
the number escaping with a 
loss. Since the experimental 


tending to raise 
moderate energy 
distribution curves show very few electrons for 
quite small energy losses, however, it is fairly 
certain that the first two peaks represent single 
inelastic impacts. 
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Neutral atom beans of noble gases of energies between 20 and 400 volts were produced by 


accelerating ionized atoms and then neutralizing the charges. The minimum energies which 


these neutral atom beams needed in order to ionize gas atoms by collision were determined. 


Neon, argon, krypton and xenon were found to ionize their own gases at energies about three 


times the actual ionization potentials. Attempts to observe ionization of one gas by neutrals 


of another have as yet given no results. 


NVESTIGATION of the minimum energies 

necessary for ionization of noble gases by 
collisions of positive alkali ions' showed three 
significant results. If we remember that a singly 
charged positive alkali ion has identically the 
same number of electrons and electron structure 
as the noble gas atom next beneath it in the 
periodic table, the results were as follows: First, 
the experimentally observed ionizing energies 
had to be multiplied by a factor m,/(m.+m,), 
this factor being the fraction of the incident ion 
energy available for outside work after deduction 
of the amount necessary for the conservation of 
momentum (m, and m; are the masses of the gas 
atom and alkali ion, respectively). The reduced 
energy thus obtained for an ionizing collision 
between an ion having A electrons and an atom 
with B electrons was almost identical with that 
obtained when the ion had the B electrons and 
the atom had the A electrons. Secondly, the 
reduced ionizing energy was less when an ion and 
an atom had the same number of electrons than 
when they had different numbers. Thirdly, if 
other gases than the noble gases were ionized at 
all by positive alkali ions, the ionization could 
not have greatly exceeded one percent of the 
efficiency of K* in A. 

The results of Beeck? and of Rostagni® on 
ionization efficiencies of noble gases by their own 
accelerated neutral atoms indicated that the 
determination of the actual inset energies for 
these processes might also yield valuable informa- 
tion, so the same method used for the alkali ion 
experiments was adapted to this problem. 

The apparatus is illustrated in Fig. 1. Gas 


'R. N. Varney, Phys. Rev. 47, 483 (1935). 
?Q. Beeck, Ann. d. Physik 19, 129 (1934). 
3A. Rostagni, N. Cimento 11, 34 (1934). 


flows from a reservoir over the heated tungsten 
filament through a deep, narrow slit, 3 mm X0.5 
mm (X4 mm depth), into a mercury pump. 
Electrons are accelerated from the filament 
toward the slit. If they produce any ions by 
collision in the field-free space inside the slit, 
these ions will be swept out at gas kinetic 
velocities by the flow through the slit. These 
ions are accelerated, as soon as they emerge from 
the slit, toward the negatively charged slit lying 
directly in line with slit one, and any electrons 
emerging from the first slit are stopped by this 
same potential.‘ Gas is allowed to stream out of 
slit number two into the pump orifice. The 
positive ions which were formed and accelerated 
before reaching slit two are now allowed to pass 
through the field-free space between slit two and 
aperture three. This region contains gas at about 
10-? mm Hg pressure. In these experiments, the 
same gas is used on both sides of the apparatus, 
so the positive ions are now traveling through 
their own gas. Kallmann and Rosen’ first 
observed the fact that when ions, particularly 
noble gas ions, travel thus through their own 
gas, they are extremely efficient in picking off 
electrons from the relatively stationary neutral 
atoms, proceeding without change of velocity 
(speed or direction) as neutral atoms, and leaving 
behing a very slow positive ion which may be 
removed by a weak cross field. Beeck® has 
reported observing, under similar conditions to 
those described above, a beam of neutral atoms 


‘ This method appears to be due to Goldmann, Ann. d. 
Physik 10, 460 (1931). It is obviously an excellent source 
of intense, monochromatic, positive gas ion beams. For 
details of the efficiency of this source see the reference 
above. 

> Kallmann and Rosen, Zeits. f. Physik 61, 61 (1930). 

§ QO. Beeck, Ann. d. Physik 19, 129 (1934). 
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Fic. 1. Diagram of tube. 


having 60 percent of the intensity of the original 
positive ion beam. No investigation was made in 
the present experiments to check the efficiency 
of neutralization, other than to observe that the 
of positive ions emerging from the 
and three 


intensity 
space between slits two 
dropped off considerably when the gas was 
admitted. A final aperture number four having a 
large positive potential assured the stopping of 
any unneutralized positive ions emerging from 
aperture three, but this potential obviously has 
no effect on the neutral atom beam. 

The neutral atoms, which are traveling with 
the energy they obtained as positive ions before 
neutralization, are now allowed to enter one of 


field-free 


the space charge cylinders described at length in 
reference 1. Between this cylinder as anode and a 
heated tungsten filament as cathode flows a 
space-charge-limited electron current. If the 
neutral atom beam ionizes any gas atoms inside 
.this cylinder, the positive ions thus produced are 
extremely effective in neutralizing the electron 
space charge, primarily because of their slower 
speed and more circuitous path. One positive ion 
has been observed to produce a space-charge 
effect equivalent to that of about 50,000 elec- 
trons.! The increased electron current reaching 
the anode as a result of this ionization by neutral 
atoms is measured by balancing out the ordinary 
that similar 


electron current against from a 


cylinder and filament arrangement outside of the 
path of the neutral atom beam. 

The neon, krypton, and xenon used in these 
experiments were purchased in Pyrex flasks as 
100 percent pure. The argon used was taken from 


N. 
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a cylinder. Spectroscopic investigations at 
Berkeley indicated that impurities in argon from 
this tank did not exceed one percent, so aside 
from a liquid-air trap to remove water and 
mercury vapors, no further purification was 
attempted. In these experiments where the effect 
observed was large, traces of impurities are 
quite unimportant since the effects of impurities 
fall off directly with the pressure. 

The ionization of each of these four gases by 
its own neutral atoms has been observed. The 
energies in electron volts at which ionization was 
observed to set in are given in column 1 of 
Table I. 

These results are probably correct to within 
+2 volts except for Ne for which +4 volts is 
Repeating the measurements 
under varying conditions, e.g., different 
sures, different potential between filament and 
the same results usually well 


more accurate. 


pres- 


slit one, gave 
within the limits given above. The errors are only 
estimated and are based chiefly on possible 
doubts of the exact speed of the beam ions at the 
time of their neutralization. 

Helium and nitrogen have been tried in this 
apparatus. Up to 400 volts no ionization was 
observed in either gas. Absolute measurements of 
ionization efficiency cannot be made with the 
space charge method, but it is estimated from 
the size of galvanometer readings that if the 
ionization of He in He had been five percent of 
that for A in A it could have been observed. 
In the case of Ne, it is possible that the experi- 
ments are confused by the fact that N* may 
have been formed rather than N.* in the ion 
source, that N* might not neutralize very well in 
going through Ne gas, and that barring the 
latter the N neutrals might not ionize Ne by 
collision as well as Nz neutrals might. A mass 
spectrograph is needed for further investigation 


of nitrogen ionization. 


TABLE I. Observed and reduced energies at which various 
gases are ionized by neutral atoms. 
OBSERVED REDUCED 
Ne in Ne 74 volts 37 volts 
AinA 48 24 
Kr in Kr 40 20 


Xe in Xe 35 17.5 
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The reduction factor mentioned above to 
correct the experimentally observed energies for 
energy losses in conserving momentum must also 
be applied in these cases. The factor is just one- 
half for all the cases in the table because the 
masses of the target and bullet atoms are equal. 
The resulting reduced values are given in column 
2 of Table I. 

These values are surprisingly low, being in 
each case only a few volts above the actual 
ionization potentials of the gases. The difference 
between these reduced potentials and the ioniza- 
tion potentials range from 16 volts for Ne to 6 
volts for Xe. The results are almost half those for 
the alkali ions in the noble gases. This is not 
surprising because the previous results indicated 
the importance of electron configurations in this 
process. A collision, for example, between K+ 
and A should be the same as between A and A 
except for the initially greater energy of the 
former combination with its positive charge 
which should make it more difficult to remove the 
electron. 

The result for A in A agrees quite well with 
that obtained by Rostagni’ at about 55 volts. 
Rostagni observed ionization of He by He 
neutrals, the efficiency being about ten percent 
that of A in A. Some doubt falls upon these 
results with helium, however, because ionization 
was indicated in the published curves at energies 
less than twice the electron ionization potential 
of He. Because of the conservation of momentum, 
which it seems certain must occur, only one-half 


7A. Rostagni, N. Cimento 11, 34 (1934). 
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of the energy of the incident atom at most is 
available in a collision between two such atoms 
of equal mass. Ionization by this mechanism 
should certainly cease therefore at a minimum 
value of twice the ionization potential. If the 
ordinates of Rostagni’s He curve be multiplied 
by ten, the He curve coincides almost exactly 
with his A curve. This implies that the helium 
may have been contaminated with argon. 
Further experiments have been conducted 
with a view to observing the ionization of one 
gas by the neutral atoms of another gas. Such 
experiments require one additional gas-filled 
chamber. In the present experiments the ions 
were neutralized in the same part of the tube in 
which they later ionized the gas. To observe the 
ionization of one gas by another, the ions must 
be neutralized in a chamber containing their own 
gas and then pass into one containing the 
different gas. Such an apparatus has been built 
and preliminary measurements miade. So far, no 
ionization has been observed by argon neutrals in 
Kr or Xe up to 250 volts. The experiments are 
not final, however, because the complications in 
forming ions in one gas, accelerating them in 
vacuum, neutralizing them in more of the gas, 
and finally passing them into another gas result 
in such a drop in intensity as to cause some 
doubts over the sensitivity of the method. Ain A 
was tested in the new arrangement and positive 
results were obtained. The efficiencies of A in 
Kr or Xe should thus not exceed 10 to 20 percent 
of those for A in A. Further experiments on 
ionization of one gas by accelerated neutral 
atoms of another may be reported at a later date. 
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The Total Ionization of Nitrogen by Electron Collisions 


GLapys A. ANSLOW AND MADELEINE DEBLOIs Watson, Smith College 
(Received February 14, 1936) 


Measurements of the total ionization produced in 1mm pressure. By dividing the total ionization produced by 
nitrogen by beams of electrons, whose homogeneity of electrons by their effective ranges as determined above the 
velocity had been tested, have been taken for electron average total ionization per cm of absorbing gas layer has 
energies up to 1500 volts. The ionization was found to _ been calculated: the results show a maximum above 100 
increase rapidly up to 75 volts, then less rapidly but volts in the region found for primary ionization per cm by 
linearly to 375 volts, for which there is a break in the Compton and Van Voorhis and by Tate and Smith. The 
ionization curve, indicating the onset of K electron ioniza- authors’ curves show a second maximum at 750 volts, 
tion; the ionization again increases rapidly near 750 volts _ twice the ionization potential of the K electrons of nitrogen, 
and after 1200 volts is found to increase less rapidly but these characteristics being similar to those discovered by 
linearly with increased electron energies. From determi- the senior author in an earlier study of the total ionization 
nations of the critical pressures for some of the electron produced in air. None of these curves agree with the results 
energies used, the mean effective ranges, the thicknesses of of Lehmann and Osgoode in the same field, the dis- 
the absorbing gas layer, have been measured and found  crepancies being ascribed to differences in the shapes of the 
nearly to obey the equation, R,=0.07+1.88 x 10-* V2 cm at ionization chambers used and in the fields within them. 


INCE the publication by the senior author' standing measurements being those of Hughes 
of determinations of the total ionization and Klein,? Compton and Van Voorhis,* Smyth,’ 
produced in air by electron collisions involving Bleakney,'® and Tate and Smith." Attention has 
energies up to 2250 volts, a limited number of been focused on this type of measurement 
experiments in the same field and involving the largely because in comparisons of experimental 
same range of electron energies have been re-_ results with the predictions of the collision 
ported, principally those of Lehmann and _ theories of classical theory’ and of wave me- 
Osgoode,? who used a method nearly equivalent chanics™ it has been difficult to evaluate the 
to that of Anslow' and accelerating voltages up amount of energy lost in secondary and further 
to 1000 volts, and those of Langmuir and Jones,’ collisions. 
who used an entirely different method and The necessity of examining the homogeneity 
voltages up to 100 volts. The best known of the velocity of the electron beam used in 
experiments with voltages in excess of 1000 ionization experiments was pointed out and 
volts are those of Eisl* and Buchmann,’ in which the homogeneity examined by Lehmann and 
the total ionization produced by swiftly moving Osgoode,“ and later by Lawrence. In_ the 
electrons was measured, and those carried out experiments, reported by the senior author'® at 
by Williams and Terroux,® who examined tracks the Washington meeting of the American 
of fast 8 particles in a Wilson cloud chamber. Physical Society in April, 1934 and described in 
In general, ionization measurements have dealt the present article, the homogeneity of the 


with that portion of the ionization produced in ~7 ffughes and Klein, Phys. Rev. 23, 111, 450 (1924). 
the first fraction of the path of the moving * Compton and Van Voorhis, Phys. Rev. 26, 436 (1925) ; 
electron, in which there is little chance of 27, 124 (1926). 

’ : . *Smyth, Phys. Rev. 35, 472, 478 (1930). 


1° Bleakney, Phys. Rev. 35, 139 (1930) ; 36, 1303 (1930). 
1 Smith, Phys. Rev. 36, 1293 (1930); Tate and Smith, 
Phys. Rev. 39, 270 (1932); 46, 773 (1934). 


secondary collisions, some of the most out- 


1 Anslow, Phys. Rev. 25, 484 (1925). 2]. J. Thomson, Phil. Mag. 23, 449 (1912); Bohr, 
* Lehmann and Osgoode, Proc. Roy. Soc. A115, 609 Phil. Mag. 25, 1 (1913); 30, 581 (1915); Henderson, Phil. 
(1927); Lehmann, Proc. Roy. Soc. A115, 634 (1927). Mag. 44, 690 (1922); Fowler, Proc. Camb. Phil. Soc. 21, 
3’ Langmuir and Jones, Phys. Rev. 31, 357 (1928). 521, 531 (1923); Phil. Mag. 47, 257 (1924). 
4 Eisl, Ann. d. Physik 3, 277 (1929). ‘8 Bethe, Ann. d. Physik 5, 325 (1930); Massey and 
> Buchmann, Ann. d. Physik 87, 509 (1928). Mohr, Proc. Roy. Soc. A132, 605 (1931) ; A139, 187 (1933); 
§ Williams, Mem. Man. Phil. Soc. 71, 23 (1927); Proc. A140, 613 (1933). 
Camb. Phil. Soc. 24, 315 (1928); Proc. Roy. Soc. A125, 4 Lehmann and Osgoode, Proc. Camb. Phil. Soc. 22, 
420 (1929); A128, 459 (1930); A130, 310, 328 (1931); 731 (1925). 
A135, 108 (1932); Williams and Terroux, Proc. Roy. Soc. 15 Lawrence, Proc. Nat. Acad. Sci. 12, 29 (1926). 
A126, 289 (1930). 16 Anslow, Phys. Rev. 45, 750 (1934). 
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Fic. 1. The filament tube and ionization chamber. 


electron beam was tested at various voltages and 
measurements taken of the total ionization in 
nitrogen produced by electron collisions with 
energies up to 1500 electron volts. 


EXPERIMENTAL PROCEDURE 


The pure nitrogen used in the experiment was 
produced according to the method of Waran,"’ 
by the action of dilute bromine water on dilute 
ammonium hydroxide, and was stored over 
phosphorus pentoxide, being allowed to leak 
through a capillary into the ionization chamber 
at a rate which maintained a nearly constant 
pressure in the chamber. 

The filament tube and ionization chamber are 
shown in Fig. 1, being the same as described in 
the earlier paper' with one alteration, necessary 
for the measurements of the homogeneity of the 
beam; the opening at the center of the hemi- 
spherical ionization chamber was decreased to 4 
mm by attaching to its lid an aluminum plate, s, 
with this circular opening, taking care to center 


'? Waran, Phil. Mag. 42, 246 (1921). 
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the aperture opposite the end of the capillary of 
the anode through which the accelerated elec- 
trons the 


ionization chamber. The capillary was 1.1 cm 


passed just before they entered 
long and 0.2 mm in diameter and the distance 
between the lid of the chamber and the anode 
was adjusted to 2 mm. 

The ionization measurements were also taken 
according to the compensation method described 
in the first experiment, the accelerating potential 
applied to the electrons being taken in all cases 
as that between the midpoint of the filament and 


the anode. 


The homogeneity of the beam 

To test the homogeneity of the electron beam 
the fraction of the electrons entering the ioniza- 
tion chamber with a velocity corresponding to 
the full drop of potential between the filament 
and the anode was determined by applying a 
retarding potential between the anode and the 
plate, s, of the ionization chamber, the pressure 
in the chamber being maintained by diffusion 
pumps at less than 10-* mm before starting a 
run. In spite of the previous outgassing of the 
metal parts, the pressure usually rose during a 
run 2X10-* mm 
production of ions, which had to be considered in 


to about with a consequent 
the interpretation of the measurements. The 
fraction of the electrons emitted by the hot 
filament and accelerated to the anode, which 
enter the chamber through the aperture, s, with 
velocities equivalent to the full drop of potential, 
was determined in the following manner, the 
Fig. 2. 


Connecting the rods and chamber together and 


electrical connections being shown in 


o— 











Fic. 2. Circuit used in homogeneity and 
ionization measurements. 
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Fic. 3. The homogeneity ratio for 108, 165 and 400-volt 
electron beams. 


to one pair of the quadrants of the electrometer 
by using the switch // in the position 1 shown in 
the diagram, and accelerating 
potential between the filament and the anode 
while the shielding cylinder and anode were 
connected to earth by switch K, used in position 
1, the number of electrons reaching the rods and 
chamber in a given time interval was measured 
by the rate of charging of the electrometer. 
Then, with switches K and // in positions 2 and 


applying an 


1, respectively, a retarding potential was applied 
between the anode and the chamber by keeping 


AND 
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the anode at a positive potential and connecting 
the chamber and the rods to earth through the 
electrometer; the latter collected both the total 
number of electrons which entered the chamber 
with velocities in excess of that corresponding 
to the retarding potential and also the positive 
ions produced by collisions with molecules of the 
gas in the chamber; the negative ions and any 
secondary electrons emitted by the walls of the 
chamber were attracted to the anode. Finally, 
with switches K and // in positions 1 and 2, 
respectively, an equal retarding potential was 
applied between the anode and the chamber 
with the anode connected to earth and a negative 
the had 
disconnected from the rods. The rods, connected 


potential on chamber which been 
to earth through the electrometer, collected the 
the the 


negative ions formed, and the secondary electrons 


electrons which entered chamber, 
produced, while the positive ions flowed to the 
chamber. 

The Townsend null method was used in taking 
the measurements, by balancing in turn the 
currents flowing to the electrometer by opposing 
currents from a potentiometer system, which 
charged a standard air condenser shunted to the 
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Fic. 4. Typical ionization-pressure curves. The ionization is proportional to the pressure up to a region 
of critical pressures, for which the thickness of the absorbing gas layer is nearly equal to the radius of 
the chamber. The pressure at C in a curve is used to calculate the mean effective range of the electrons. 
Total ionization values for a given energy are taken as the mean of the values on the horizontal part of 
the curve. 
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ras_e I. Test of homogeneity of electron beams 


ACCELERATING 
POTENTIAL 400 volts 165 volts 


RETARDING 


POTENTIAL 310 volts 152 volts 
Pressure (10-* mm 1.2 0.8 1.5 1.8 1.1 2.0 
r) 0.74 0.97 0.60 0.59 0.61 0.67 


ro 1.24 1.00 1.44 1.17 1.16 1.06 
Homogeneity ratio 0.99 0.98 1.02 O.88 O.88 0.86 
mean 1.00 mean 0.87 


electrometer. Since the pressure on the chamber 
side of the anode rose slightly and the electron 
beam varied somewhat in intensity during a set 
of readings, the three measurements to be made 
were taken in order and then the first and the 
second repeated in reverse order. If 7; represents 
the ratio of the current due to the retarded 
electrons and the positive ions within the 
chamber to the unretarded electron current, and 
r. the ratio of the current due to the retarded 
electrons, the negative ions and the secondary 
electrons to the same unretarded electron current, 
the homogeneity ratio may be taken as the mean 
of r; and re. Two sample sets of results are given 
in Table I, and the values of the homogeneity 
ratios for different retarding potentials are shown 
graphically in Fig. 3 for 108, 165 and 400-volt 
electrons. It is evident that these electron beams 
possessed a high degree of homogeneity. 


Ionization measurements 


The ionization produced by electrons of a 
definite energy was measured exactly as in the 
previous experiment. For a number of voltages 
the ionization was determined for a series of gas 
pressures: typical ionization-pressure curves are 
shown in Fig. 4, and indicate that for small 
pressures a linear increase in ionization occurs 
with increasing pressure up to a region of 
critical pressures, for which the radius of the 
chamber is comparable with the mean thickness 
of the absorbing gas layer traversed; for still 
higher pressures the ionization is constant, this 
constant value giving the total ionization 
measurement. Determinations of the total ioniza- 
tion produced were made at frequent voltage 
intervals; the values obtained are given in the 
first two columns of Table II, and the results 


are shown in Fig. 5. 
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Fic. 5. Total ionization curve for electrons with energies 
ranging from 0—1500 volts 


It is evident that the total ionization produced 
is not a linear function of the energy of the 
colliding electron for energies less than 1000 


TABLE II. Measurements of total ionization produced by 
electrons of vartous energies in nitrogen. 


Ions Per 
ELECTRON 
PER cm OF 
EFFrectivt 


RANGE Al 
ENERGY Ions PER EFFECTIVE RANG! 1 mm 
volts ELECTRON AT 1 mm PREssuR! PRESSURI 
22 0.11 *0.071 1.55 
42 0.73 0.073 9.95 
60 1.20 0.077 14.2 
108 1.72 *0.092 18.7 
165 2.25 *0.123 18.3 
225 2.78 0.165 16.9 
306 3.59 *0.241 14.6 
321 3.78 0.264 14.3 
350 3.85 0.316 12.3 
376 4.06 0.335 12.1 
386 4.18 0.350 11.9 
401 4.17 *0.372 11.2 
412 4.20 0.398 10.5 
487 4.46 *0.513 8.6 
514 4.39 0.574 7.65 
560 4.88 *0.640 7.40 
680 6.20 0.939 6.60 
690 6.52 0.960 6.79 
700 7.05 *0.975 7.23 
732 10.4 1.097 9.66 
754 11.5 1.110 10.40 
804 13.2 1.255 10.52 
870 15.1 *1.470 10.17 
935 17.1 1.693 10.10 
990 18.3 *1.880 9.65 
1030 19.2 2.053 9.36 
1100 19.5 2.275 8.54 
1120 20.0 2.348 8.52 
1240 21.2 "2.400 7.71 
1460 23.1 *3.618 6.36 


Note. The values of the effective range marked * were obtained from 
the ionization-pressure curves, the others by interpolation from the 
range curve in Fig. 8. 
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volts, contrary to the results of Lehmann and 
Osgoode ;? it is believed that this lack of agree- 
ment can be attributed to two important 
differences in the experimental conditions. Firstly, 
Lehmann and Osgoode used a cylindrical rather 
than a hemispherical ionization chamber. In the 
former type of chamber the electrons deflected in 
primary collisions may hit the sides of the vessel 
before their energy is exhausted, since their 
ranges at the pressures used are greater than the 
radius, though than the length of the 
cylindrical vessel. In the hemispherical type of 
chamber these deflected electrons spend all their 
energy in collisions before reaching the sides of 
the chamber. Secondly, in the experiments of the 
authors the anode was kept at earth potential, 
thus maintaining a field free region within the 
ionization chamber; the other 
experiments the filament was at earth potential 
and the anode at a positive potential, thus 


less 


whereas, in 


causing a retarding field between the anode and 
the ionization chamber, which decreased the 
energy of the electrons passing through it. It is 
believed that these variations in the shape of the 
ionization chamber and in the field within it are 
mainly responsible for the discrepancies between 
the total ionization and the range curves here 
reported and those of Lehmann and Osgoode. 
Fig. 5 indicates that the ionization of nitrogen 
sets in at about 17 volts, the ionization potential 
of the nitrogen molecule,'* increases rapidly to 
about 75 volts, then less rapidly but linearly up to 
375 volts, the ionization potential of the K 
electrons!’ of nitrogen, at which there is a 
sudden decrease in the number of ions produced 
because of the large amount of energy expended 
in releasing one of the K electrons, leaving less 
energy available for the secondary collisions. 
The fact that a bend in the ionization curve is to 
be expected has been emphasized by Darrow,” 
and probably is more evident in total ionization 
measurements, the effects of 
secondary collisions, than in primary ionization 


which include 
experiments, in which these are eliminated. The 
break has never been recognized in any of the 
latter experiments. 


18 Smyth, Proc. Roy. Soc. A104, 121 (1923). 


1° Mohler and Foote, Sci. Papers Bur. Standards 17, 471 
(1922). 
*0 Darrow, Electrical Phenomena in Gases (1932), p. 43. 
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Fic. 6. Total ionization curve for nitrogen and air. 
Both show breaks at 375 volts and similar rapid increases 
in ionization near 750 volts. 


As the energy of the accelerated electron 
approaches that equal to twice the ionization 
energy of the K electrons, the probability of this 
type of ionization increases rapidly, and finally 
a second linearity is approached at and beyond 
1200 volts. 

It is interesting to note that these features were 
some of the characteristics of the total ionization 
curve for air given in the previous paper, and 
which is shown with the nitrogen curve in Fig. 6; 
the breaks at 250, 375, and 495 volts in the air 
curve were attributed to argon, nitrogen, and 
oxygen, respectively, since these are their 
ionization potentials.*': '*- 2? The rapid increase 
in the ionization of air 1000 volts was 
attributed in the earlier paper to the ionization 


near 


of oxygen. 


Range calculations 


Since the path of an electron which has suffered 
multiple collisions is tortuous and in three 
dimensions, the mean depth of the absorbing gas 


layer, known as the effective range, is consider- 


*1 Int. Crit. Tab. 6, 71-72 (1929). 
*? Kurth, Phys. Rev. 18, 461 (1921). 
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Fic. 7. lonization-pressure curves for small energy electrons 
in air. 


ably smaller than the actual length of path 
traversed by the electron. In a separate paper 
the senior author* has given the derivation of a 
formula relating these quantities, and reported 
numerical calculations of the actual ranges, made 
by suitably choosing constants which appear in 
Bethe’s" equations for the scattering of electrons 
in both elastic and inelastic collisions. 

In this experiment a mean value of the 
effective range of the electrons may be obtained 
from an examination of the region of critical 
pressures in their ionization-pressure curve. 
Those electrons, for which the fraction of energy 


Anslow, Phys. Rev. 49, 480A (1936). 
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lost in collisions is consistently large, are scat- 
tered most; when the pressure in the ionization 
chamber is such as that at the point A in the 
curve for the 700-volt electrons in Fig. 4, the 
ends of their paths must be on the surface of the 
ionization chamber and the radius of the chamber 
will equal the effective range of these electrons. 
Similarly, the point B on the same curve repre- 
sents a pressure for which the radius of the vessel 
is very nearly equal to the effective range of those 
electrons and for which scattering is a minimum 
and straggling a maximum. For some intervening 
pressure the radius of the chamber must equal 
the mean value of the effective range of the 
electrons in the beam, such that half the electrons 
have effective ranges in excess of this value and 
the others less. Since the region of critical 
pressures is very narrow for the lower energy 
electrons and amounts to less than ten percent of 
the mean pressure within the critical region for 
the highest energy electrons used, a mean value 
of the effective range of the electrons in a beam 
may be obtained by calculating the pressure, C, 
at which the linear portions of its ionization- 
pressure curve intersect: multiplying this pres- 
sure in mm by the radius of the chamber gives 
the effective range in cm at 1 mm pressure. The 
error involved in such a determination cannot 
be large for energies somewhat in excess of those 


used in this experiment, for an inspection of the 
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Fic. 8. Effective range vs. voltage relation for energies up to 1500 volts in N» and air. 
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Fic. 9. Effective range vs. voltage relation for energies up 
to 175 volts in Ne and air. 


similar curves obtained by Buchmann’ for elec- 
trons with energies from 4—13 kv and by Eisl* for 
energies from 9.8—3.6 kv indicate that the region 
of critical pressures remains narrow up to at 
least 20,000 volts. The region narrows down to 
practically a single pressure for beams with 
energies less than 150 volts, as can be seen in 
Fig. 7, where are plotted three of a series of 
curves obtained in a detailed study of the range 
of such electrons in air, measured since the 
publication of the previous paper. Since electrons 
with these energies have ranges scarcely twice 
the mean free path of electrons in the gas, they 
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participate in only one or two collisions and those 
in a single beam must suffer practically equal 
scattering. Hence for these low energy electrons 
a single critical pressure should be anticipated. 
The effective ranges thus obtained were plotted 
against the square of their energy in volts for 
energies up to 1500 volts, the results being shown 
in Fig. 8, together with the similar result obtained 
in the previous experiment for air. In Fig. 9 the 
same is shown for energies less than 175 volts. 
It is evident that there is a fair linear relation 
between the effective range of the colliding 
electron and the square of the energy expended 
for ranges less than 2 cm at 1 mm pressure, with 
a slight falling off in the range values for higher 
energies, and that there is a deviation from this 
linearity for small energy collisions in air. The 
empirical equation, derived from the lower 
portion of the curve in Fig. 8, between the mean 
effective range, R,., measured in cm at 1 mm 
pressure, and the electron energies, measured in 
electron volts, V, is for nitrogen 
R.=0.07+1.88X10-° V2. (1) 
The slope of this curve, 1.88X10~°, decreases 
with the higher energies, and probably approaches 
the value obtained by Williams* for slow £- 


24 Williams, Phil. Mag. 2, 1109 (1926). 
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Fic. 10. Average total ionization per cm of absorbing gas layer for Nz and air, calculated by divid- 
ing total ionization measurements by the mean effective range of the ionizing electron. 
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particles in nitrogen, which reduced to 1 mm 
pressure, is 1.54X10-®. The effective ranges of 
the various electron beams, estimated from this 
curve, are given in the third column of Table IT. 


Ionization pet cm of absorbing gas layer 

The average ionization produced per cm of the 
absorbing gas layer was obtained by dividing 
the experimentally determined value of the total 
ionization by the effective range of the colliding 
electron. The values obtained are given in the 
last column of Table II and plotted against the 
energy of the accelerated electrons in volts in 
Fig. 10, together with the same for air as given 
in the first paper, the results for small voltages 
having been calculated from the later determi- 
nations of the range in air. It is seen that in the 
curve for nitrogen a maximum value of the 
average ionization per cm of the effective range 
occurs between 100 and 150 volts, a break at 375 
volts, and a second maximum at 750 volts, this 
latter maximum being at twice the ionization 
potential of the K electrons of nitrogen. In the 
air curve this maximum came at twice the K 
ionization potential of oxygen. 

In Fig. 11 the characteristics of the nitrogen 
curve are compared with the well-known curves 
for primary ionization per cm of Compton and 
Van Voorhis* and of Tate and Smith.”° At the 
low pressure used in these experiments not more 
than one collision is probable in the electron 
path, and the thickness of the absorbing gas 
layer is identical with the actual distance 
traversed by the electron. In the figure there are 
also given graphs of the results obtained in 
Kossel’s** early experiment and values calculated 
by the present authors from the total ionization 
and range curves published by Lehmann and 
Osgoode. It is evident that both total and 
ptimary ionization measurements indicate a 


** Tate and Smith, Phys. Rev. 39, 270 (1932). 
26 Kossel, Ann. d. Physik 37, 393 (1912). 
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Fic. 11. Average total ionization per cm of effective 
range compared with primary ionization per cm results for 
actual range of Compton and Van Voorhis (C-V), Tate 
and Smith (T-S), Kossel (K), and average total ionization 
per cm of effective range calculated from Lehmann and 
Osgoode’s determinations (L-—O). 


maximum probability of ionization for energies 
slightly in excess of 100 volts, but that the second 
maximum at 750 volts does not appear in any of 
the curves except that of the authors. In the 
experiments of von Hippel®? with mercury and 
Funk®* with sodium and potassium similar 
maxima have been located at approximately 
twice the ionization potentials involved. The 
form of the Lehmann and Osgoode curve is at 
almost complete variance with the general type 
of curve obtained in these experiments, and 
strengthens the belief of the authors that the 
Lehmann and Osgoode results were seriously 
affected by the shape and the field within their 
ionization chamber. 

A reduction of the authors’ results to the 
ionization per cm of the actual range, together 
with a comparison with the values predicted by 
Bethe’s theory may be presented in a future 
communication by the senior author. 


7 von Hippel, Ann. d. Physik 87, 1035 (1928). 
*’ Funk, Ann. d. Physik 4, 149 (1930). 
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The noise level of a system composed of an ionization 
chamber and a linear amplifier is calculated as a function of 
the input circuit parameters and the characteristics of the 
amplifier. The frequency spectrum of a pulse caused when 
an ionizing particle passes through the ionization chamber 
is also calculated as a function of the circuit parameters 
and the collection time. From this analysis the shape of the 
output pulse is calculated as a function of the frequency 


OST of the information available concern- 

ing the design of ionization chambers and 
their associated circuit and pulse amplifier have 
been of a strictly empirical nature. It is the 
purpose of this analysis to indicate more defi- 
nitely the important factors which limit the 
sensitivity of this type of detector in the in- 
vestigation of the products of nuclear dis- 
integrations. The problem is essentially one of 
determining the effect of circuit parameters and 
the characteristics of the amplifier used upon the 
signal to noise ratio. 

In the design of an amplifier for the measure- 
ment of voltages whose amplitude is of the same 
order of magnitude as the noise of the first tube 
and its associated input circuit, the ratio of 
signal voltage to noise voltage depends upon 
their relative frequency spectrum and the trans- 
fer characteristics of the amplifier. In order to 
design an amplifier for the measurement of pulses 
from an ionization chamber, therefore, it is 
necessary to determine the frequency spectrum 
of a pulse, and the frequency spectrum of the 
noise from the input circuit and first tube. 
Before considering the limitations imposed in 
practice by tube characteristics, these two fre- 
quency spectra will be determined. 

There are three principal sources of noise! 
from a circuit of the type shown in Fig. 1 
which is the equivalent circuit of an ionization 
chamber and the first stage of the amplifier. 
The noise is due to thermal noise from the 
resistance, shot noise from the grid current in 
the resistance-capacity network, and the zero 
noise of the tube. In Fig. 1, C is the combined 


1G. L. Pearson, Physics 5, 233 (1934). 


characteristic of the amplifier. It is conéluded that the 
cut-off frequencies of the amplifier for optimum pulse shape 
and amplitude depend upon the collection time. The signal 
to noise ratio is calculated for the case of optimum pulse 
shape. By proper choice of tubes it is shown that it should 
be possible to count 100 ion pairs and that for a 1-cm deep 
chamber and fast protons the signal to noise ratio should 
be 38 : 1. 


capacity of the ionization chamber, tube, and 
wiring, and R is the resistance in the grid circuit 
due to the internal dynamic grid resistance and 
any external grid resistance. 

In the analysis that follows an ideal amplifier 
is assumed which has constant gain between the 
frequencies f; and fz and zero gain outside of 
this band, and no phase distortion in the band. 
Such an amplifier is difficult to obtain in practice, 
but the conclusions reached will be only slightly 
affected by the approximation. 


CALCULATION OF NoIsE LEVEL 


The thermal noise from the input circuit due 
to thermal agitation of charge in the resistance 
is given by 


E*=4kTR { df (1+47*R°C*f*), (1) 
“si 
where K is Boltzmann’s constant, 7 is the 
absolute temperature. By integrating 


E?=(4kT/27rC) tan“ 
[2eRC(fe—fi)/(1+42?R°C*Ffof;) J. (2) 


The shot noise in the input circuit due to 
statistical fluctuations in the grid current J, of 
the first tube is given by 
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Fic. 1. Input stage of pulse amplifier. 
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. 


E=2e1,R° | df/(1+4R°Cf2), (3) 


vf 
E,*=(2el,R/2xC) tan 
[2rRC(fo—f:)/(1+4e2R2C%fof:)], (4) 


where e is the electronic charge and J,=|J+ 
+\|I—| where J+) and |J—| are the absolute 
magnitude of positive and negative components 
of grid current. From (2) and (4) we obtain the 
relation for T= 300°K, 
E,*=(19.41,R)E,*. (5) 
The third source of noise is due to the zero 
noise level in the tube itself and is given by: 
Ee=4kT | Raf (6) 
vst 
=4kTR:(fe—fi), (7) 
where R, is the resistance which, if placed in the 
input of the first tube, would give the same 
reading in the output meter of the amplifier if 
there were no noise from the tube itself. The 
total mean square voltage is then given by 
E,Z=EZ+E?+E-?, 
E,*=4kT}(1/2rC) tan 
-2eRC(fe—fi)/(1 +42? R?2C*f fe) | 
X(1+19.47,R)]+Rilfe—fi)}. (8) 


The case of interest to us is when 22f2RC 
>2nf:RC>1 and for this condition Eq. (8) 


reduces to 
E,2=4kT(fo—f;) 
XLRi+(1+19.47,R)/RC*An*f fo). (9) 


If we plot £,? vs. R from Eq. (8) we obtain a 
curve of the type shown in Fig. 2 where the 
circuit constants assumed are: f;=500 c.p.s., 


fe=5000 c.p.s., C=10-" farad, J,=10-'° ampere 


and R,= 20,000 ohms. 

We see that when R is very small the zero 
noise from the tube itself predominates and 
determines the noise level. As R is increased the 
thermal noise increases and is the determining 
factor and would increase indefinitely except for 
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AND PULSE 


ut ™ 


Fic. 2. Input noise as a function of input resistance. 


the shunting effect of the capacity C. Because of 
the shunting effect of C which becomes more 
effective as RwC is large compared to unity, 
the thermal noise decreases and the total noise 
would again be reduced to the noise level of the 
tube alone were it not for the fact that the shot 
noise due to the grid current becomes of im- 
portance. The shot noise decreases with an 
increase in R to an asymptotic value and the 
minimum noise is given by 


E,2=4kT(fe—fi)[Ri+19.41, 42°Cfi fe] (10) 


and is determined only by the zero noise of the 
tube, by its grid current and by the parameters 
C, fi, and fz and does not depend upon R if R is 
sufficiently large. 

From (10) we conclude that the noise level is 
made a minimum by choosing a tube with as 
low values of grid current and of zero noise as 
possible and making R as large as possible. 
The parameters f;, fe, and C also determine the 
noise level but we shall see that they are fixed 
by the requirements of pulse definition. 


FOURIER ANALYSIS OF PULSE 


The voltage on the grid of the first tube due 
to the discharge of the ionization chamber is 
given by? 


—x <t<Q, 


E,=0, 


E,=(Rq A\)(l—e VF°), O0<t<d, (11) 
E, =(Rq/d)(1—eN RO) eIRE, NES 4, 

2G. Ortner and G. Stetter, Zeits. f. Physik 54, 449 
(1929). 
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where gq is the charge collected on the plate of 
the ionization chamber connected to the first 
stage of the amplifier, due to the passage of a 
single ionizing particle through the chamber, 
and X is the time required to sweep the ions to 


the plate. 


fe+co 


E,,(t) is continuous and the integral | E,(t)dt 


2 x 
is absolutely integrable. Thus the conditions for 
a Fourier integral are satisfied. The spectrum of 


(11) is then given by 


d 
G(w) = (1, (27))(Rq »| [ e “dt 


aX Ace 
sit | e— (U/RC+jw) ttt ve | e (URC +ja) ty 
e e 


r 


-{ e(URC+ie at.| (12) 
r 


Integrating and substituting the limits we have 


G(w) = A(1—cos wA+/ sin wd) /jw(1+jwRC) 


=A(1—e-/*) /jw(1+jwRC), (13) 
where A = Rq/(27)!n. 

We note that the frequency spectrum is of the 
oscillating type with an envelope that approaches 
Rq/(27)} asymptotically at the origin. The 
average energy per cycle consequently falls off 
at the higher frequencies and it should be noted 
that there is no maximum in the frequency 
spectrum corresponding to a frequency 1/4, as 
is commonly supposed. The effect of varying X is 
merely to change the number of oscillations 
inside of the envelope of the energy spectrum 
within a given frequency range. The effect of 
changing RC is to shift the envelope to higher 
or lower frequencies as RC is decreased or in- 
creased. 

The best means of evaluating the effect of 
the circuit parameters is to determine the time 
response in the output of the amplifier. The time 
response in the output of the amplifier referred 
to the input, i.e., divided by the gain of the 
amplifier, is given by 


EW =A f (1 —€7?* )e?*"dw/jw(1+jwRC). (14) 
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It can easily be shown that the value of E(t) 
will be a maximum when weRC>w,RC>1 and 
(14) becomes * 


y2 


E(t) =A(t—2) RC| cindy /y? 


—AU/RC{ &*dx/x*, (15) , 
where x= w/t and y=w(t—A). By integrating and 
substituting the limits, 4 
E(t) ={A(t— d) RC} { [cos w:(t—X) ] wi (t—d) 

—_ [cos we(t— d) ] weo(t—d) + Siw;(t— ) 
— Siwe(t—r)}—{At/ RC} { [cos wit | ‘wit 
—[cos wet ] wet+ Siw t— Siwet} 
{7A (t—d) /RC} { [sin wi(t—d) ] /wi(t—2) , 
—[sin wo(t—X) ]/we(t—d) — Ciw;(t—d) 
+ Ciws(t—r)} — | jAt/RC} {[sin wit) wit ; 
—_ [sin wot | wot — Ciw t+ Ciwst! ° (16) 
in any practical case the amplitude of the com- 
plex part is small and can be neglected, and (16) 
can be written: ’ 
E(t) =(q/Cr(2r)! Lt X — Y)+aV] 
=[¢/C(2m)! JE(A, 2), (17) 
where 
' 
E(x, ()=(4X-—Y)+AVY]/A, X=X1-Xze, 
X n=(cOs wrt |/wnt+Siwnt,*¥ Y=Y,— Fz, 
Y,,=[cos wn(t—d) ]/wn(t—A) + Siw,(t—d). 
Then the ratio of signal to the root-mean-square 
value of noise is given by 
E(t)/E,=qE(a, t)3.14 XX 10°/((fe—fi) (18) 
X[R.C?+(1+19.47,R)/42°Rfife])!. 
One of the additional requirements in counting 
’ 


pulses is that the amplitude of all except the first 
pulse in the train of output pulses due to the 
causal input pulse be small. This requirement is a 
fundamental one, since otherwise a single pulse 


* The best tables of this functien are found in Vol. I of 
the Mathematical Tables published by the British Associa- 
tion for the Advancement of Science. 











IONIZATION CHAMBERS 
might be counted more than once in the output, 
and the sharpness of the first output pulse is a 
deciding factor in the resolution of the counting 
system. 
From a 
determine the relative noise level when all of 


consideration of (18) we can now 
the parameters are given. Before discussing this 
equation it is necessary to consider the function 
E(x, t). 


E(x, t) depends only upon f;, fe, A, and ¢. 





"2 
‘ 


AND PULSE AMPLIFIERS 1 
In order to determine its dependence upon these 
quantities E(\, t) has been plotted for a large 
number of different combinations of the param- 
eters. A typical family of curves of E(A, ¢) is 
shown in Fig. 3 for various values of f. with 
f,;=240 c.p.s. and \=2X10~ second. It is to be 
noted that only the first pulse is essentially 
affected by the upper band limit, the secondary 
and negative pulses are primarily determined by 
the lower band limit. This does not hold exactly 
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Fic. 5. Pulse shape as a function of collection time 
(fi = 240 c.p.s.; fe= 2400 c.p.s.). 

Fic. 6. Pulse amplitude as a function of collection time 
(fo= 10f1; To9=2/fe). 


when f2 is very near to f,; but does for the cases 
in which we are interested, since from Fig. 3 it 
is evident that f must be large compared to /f;. 
However, it is seen that no appreciable change 
occurs if f, is increased beyond 2400 c.p.s. in 
this particular case. Similar families of curves 
are obtained for different values of f; and \ and 
in each case the critical value of fe occurs at 
fo=1/2X. 

Fig. 4 shows a similar family of curves of 
E(x, t) vs. ¢ for various values of f; and with 
fo=2400 c.p.s. and }\=2X10~ second. From 
these curves and similar families we conclude 
that the optimum value of f; is obtained when 
fi equals 1/20. 

The optimum pulse shape is thus seen to be 
obtained when f2 equals 10 f; and fe=1/ 2X. Fig. 5 
is a plot of E(A, ¢) vs. t for various values of \ and 
with f; = 240 c.p.s. and fz= 2400 c.p.s., and shows 
the effect of variation of collection time for a 
fixed band width. The same data are plotted in a 
different way in Fig. 6 in which the maximum 
amplitude of E(A, ¢) is plotted against \/A» where 
Ao=1,/ 2f2. This curve is perfectly general and it 
is seen that no further change in pulse amplitude 
occurs if \ is less than the critical collection 
time Ao. From this discussion we see that if the 
longest collection time of the ionization chamber 
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is known the characteristics of the amplifier are 
determined. 

In any practical case the collection time X is 
determined by the geometry and the gas of the 
ionization chamber as the collection 
voltage. It should be noted that the curve of 
Fig. 6 provides a means of determining the 


well as 


collection time of a chamber, for since \ can be 
varied by varying the collection voltage across 
the chamber, the voltage at which the critical 
collection time is obtained can be determined, 
and the collection time at any other voltage 
found from the curve. 


If we assume the amplifier fe=10f,; and 


fe=1/2d and write for E(A, t) its maximum value 


we obtain as a measure of the signal to noise 


ratio 
E,, (t)/ Ex =q7.3 X10°(A, [RiC2+ 


x (1+19.47,R)/0.99R})!. (19) 


If we differentiate (19) with respect to \ and 
set equal to zero in order to find the optimum 
value of \ for a given set of tube and circuit 
parameters, we find that (19) is a maximum 
when X has the value : 


A= C(0.99R.R (1+19.47,R))}. (20) 
The maximum signal-to-noise ratio is then 
given by 


[En (t)/Enjmax=N11.3X10-" 


X(R/R.C7(1+19.47,R))', (21) 


where JN is the equivalent number of singly 
charged ion pairs collected by the chamber. 
If we regard R, and J, as parameters fixed by 
the tube from (21) the signal-to-noise ratio is 
made a maximum when C is a minimum and R 
is large. If the tube is operated with a floating 
grid the factor (R/R,C?(1+19.47,R))! is a meas- 
ure of the merit of the tube. In this particular 
case R is the dynamic grid resistance at the 
floating potential, and C is the sum of the 
ionization chamber capacity and the dynamic 
input capacity of the tube, and is given nearly 
enough by 

C=C.+Cos+(1+5) CoC. 


(22) 


where: 
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C; =ionization chamber capacity, 
C,;=grid-to-filament capacity, 
Cp, = grid-to-plate capacity, 

C,=stray capacity to ground, 

ue =effective voltage amplification of first stage. 


In triodes R,; may be as low as 4000 ohms' 
while in pentodes the best values of R, usually 
fall between 10,000 to 20,000 ohms.': * However, 
in triodes the value of (1+ 4.)C,, is usually 
prohibitively large so that pentodes are to be 
preferred in spite of their higher values of R;. 

For example, if we assume an ionization 
chamber plus stray capacity of 5 mmf, the 
relative merit of the RCA 238, the WE 259-B 
and WE 262-A has been calculated from pub- 
lished data and is given in Table I. 

The 238 and 259-B are pentodes and the 262-A 
is a triode. There are a large number of other 
types of tubes, both pentodes and triodes with 
factors of merit in the same range of values: 
and the pentodes in general have the higher 
values. 

Since the collection time \ determines the 
amplifier which in turn determines the resolution, 
for fast counting rates it may be desirable to 
operate at a different value of \ than the one 
that gives the best noise level. For this case from 
(19) we see that the best ratio of signal to noise 
is still obtained when R, and C are a minimum 
and R is large. When counting fast protons it 
may not be necessary to keep the secondary 
extremely low compared to the first pulse, and 
in that case the signal to noise may be improved 
slightly by the use of a narrower band width. 

By a careful choice of tubes it may be expected 
that the factor of merit of the input tube may be 
increased beyond the values given in Table I. 
If the circuit of Fig. 1 is used and a tube with the 
constants R,=10,000 ohms, C=10-" farad, 
I,=10-" ampere, and an input resistance of 
10" ohms is used, the factor of merit will be 


TABLE I. Relative merit of several input tubes for 
pulse amplifiers. 


( R Ig (R/RiC%1 +19.41,R))4 
TUBE Ry x10'? K 107" «10! x10~° 
RCA 238 14000 14 § 0.2 8.9 
WE 259-B 19800 11 1.2 0.8 6.5 
WE 262-A 7700 28 0.8 1.8 4.3 


’E. A. Johnson and C. Neitzert, Rev. Sci. Inst. 5, 196 
(1934). 





75 


AMPLIFIERS | 








Fic. 7. Grid voltage vs. grid current curve of RCA type 
38 tube. 


24X10°. From Eq. (21) if we assume that pulses 
whose amplitude is twice the mean noise level 
can just be detected, the minimum number of 
ion pairs that can be counted will be approxi- 
mately 100. If we assume that fast protons 
produce 2000 ion pairs per centimeter of path, 
the signal-to-noise ratio for a one-cm deep 
ionization chamber will be 38 : 1 from Eq. (21) 
for fast protons and will be higher for slow 
protons and alpha-particles. 

The circuit proposed in Fig. 1 has several ad- 
vantages. If we consider the curve of E, vs. I, of 
the 238, from Fig. 7 we see that by operating at 
the point E, = —2 volts instead of at the floating 
potential the grid current is reduced to about 
half, since the electronic component is reduced 
to practically zero. The resistance of the circuit 
is practically the external resistance at E,= —2 
volts since the dynamic resistance is very large. 
The noise level can consequently be improved 
slightly by the use of this circuit. It has one 
further advantage as is shown by the following 
analysis. 

If RC is large and the number of impulses 
large, the grid will assume a voltage different 
from the voltage with no signal impressed. We 
assume first that the resistance in the grid 
circuit remains constant. It must be pointed out 
that this may be far from the actual circum- 
stance, and may seriously affect the operating 
conditions. 

We assume that each pulse decays as 


Ee“V®° where E=Ra(t—e V/FC)eRC/Q, (23) 


If the pulses occur at an average interval of 1/n, 





176 Ge i 


then the voltage of the grid at any time ¢ will be 
'=Efe t/RC 4 @—(t—1/n) RC 4 ¢ (t—2/ n) RC4...), 


Assuming a large number of pulses, the grid will 
be driven to a voltage: V= E/(e!/"®©-—1) =nRCE 
since RC will usually be large compared to 
unity. 

V may be of the order of 0.1 volt if the counting 
rate is high, when counting alpha-particles. With 
a floating grid, if we consider the curve of E, vs. I, 
of Fig. 7 it is seen that if the sense of V is such 
as to make the grid more negative J, is decreased 
and R, is increased. This results in an increase 
in V which carries the operating point still 
further negative, and the process continues until 
the grid voltage is driven so far negative that 
the tube becomes inoperative. For this reason 
the ionization chamber collection voltage is 
applied so that V drives the grid positive. 
However, this increases J, and decreases R, and 
the grid voltage is established at a new equi- 
librium at a more positive voltage, with J, 
greater and R, less than with no signal impressed. 
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Both of these factors result in a decrease in the 
signal to noise ratio which may be quite 
appreciable. 

The circuit of Fig. 1 can be used to prevent 
this condition either by manually correcting for 
the change in grid bias at high counting rates, 
or by automatically biasing the grid. 

In conclusion it should again be pointed out 
that this analysis holds only when an amplifier 
which has no phase distortion is used. If unequal 
phase changes occur between the frequency 
components during amplification, they will not 
add in the output to give a maximum signal, 
but interference will take place which will result 
in a decreased signal. It is therefore essential 
that the overall characteristic of the amplifier 
shall contain no phase distortion. 

The authors are indebted to Dr. John A. 
Fleming, the Director of the Department of 
Terrestrial Magnetism, for the opportunity to 
carry on this investigation, and to Dr. L. C. 
Van Atta and Dr. L. R. Hafstad for criticism of 
the manuscript. 
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Magnetization of Nickel Under Compressive Stresses and the Production of 
Magnetic Discontinuities 


C. W. HEAPS 


, Rice Institute 


(Received May 26, 1936) 


When nickel is compressed its hysteresis loop approaches a rectangular shape. Several curves 
for stresses increasing up to the elastic limit have been obtained. At the maximum stress the 
steep part of the loop is still not vertical and the Barkhausen effect is not appreciably changed 
by the pressure. Much smaller stresses in an elastically bent wire produce large magnetic 
discontinuities, hence the existence of contiguous compressed and stretched regions appears 
to be necessary for the production of the large discontinuities. 


HEN a straight piece of nickel wire is bent 

into the arc of a circle large discontinuities 
of magnetization are found to appear in its 
hysteresis loop.' It has been supposed that a 
discontinuity of this kind is caused by the sudden 
and complete reversal of the saturated mag- 
netization of the compressed region of the wire. 
The compressed nickel is believed? to have an 
approximately rectangular hysteresis loop, while 


'R. Forrer, J. de physique 7, 109 (1926); 10, 247 (1929). 
2 M. Kersten, Zeits. f. Physik 71, 553 (1931). 


the region of the wire which is under tension has 
an oblique loop with hardly any knee visible in 
the respective branches. The resultant loop for 
the entire wire would be obtained by super- 
position of the two loops. 

The compressed half of the material of the 
wire maintains an almost constant intensity of 
magnetization as the field is decreased from a 
saturation value to zero; during this same field 
change the stretched part of the wire decreases 
its magnetization from saturation to a relatively 
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Fic. 1. Hysteresis loop of a nickel wire subjected to an 
elastic bending strain. 


small value. The resultant effect, therefore, is 
that the residual magnetism for the entire wire 
should be about half the saturation value. Fig. 1 
shows the experimental results obtained by 
M. Kersten for a nickel wire under bending 
stress. Not only does the residual magnetization 
have approximately the half-saturation value, 
as expected, but according to Kersten the part 
AB of the curve is quite accurately given by the 
theory. 

In order to understand more perfectly the 
nature of the processes involved in the discon- 
tinuity BC it is desirable to know more about the 
behavior of nickel under high, purely compressive 
stresses. The rectangular hysteresis loop with 
magnetic discontinuities has not been observed 
experimentally. Experiments have been per- 
formed on compressed nickel bars by Ewing.’ 
His data show that the loop tends towards a 
the compressive 
are no 


rectangular form as stress 
increases, but there 
However, Ewing used the method of reversals 
and would therefore not discover a discontinuity 
even if it existed. This consideration made it 
seem desirable to repeat Ewing's experiment. 
Commercial nickel was obtained in the form 
of a bar, specified by the makers as “‘pure, grade 
A, annealed.”’ It was turned down into a rod 
10.0 cm long and 0.795 cm in diameter. This rod 


discontinuities. 


3J.A. Ewing, Phil. Trans. 179, 333 (1888). 
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was fitted into a yoke of soft iron, the cross- 
sectional area of the yoke on each side of the 
nickel bar being 56.2 cm?®. Each end of the bar 
extended 2.5 cm into the yoke so that a length 
of 5 cm was available for the coils. An inner 
induction coil of 478 turns was used. Over this 
was placed a magnetizing coil of 286 turns. 

A hard steel plunger was put through a hole 
in the yoke so as to rest on top of the vertical 
nickel bar, and by means of a lever resting on 
the steel plunger various pressures could be 
applied by the use of weights. The hysteresis 
curve was obtained by using a ballistic gal- 
vanometer and the step-by-step method. In the 
region where discontinuities might be expected a 
liquid resistance could be used to vary the 
current slowly and continuously. 

Fig. 2 shows one branch of each of the various 
hysteresis loops obtained. A common saturation 
value of J=450 was obtained for all the curves 
in a field of about 120 gauss. No large discon- 
tinuities were detected for any of the curves. 
Furthermore, the use of a telephone receiver and 
amplifier showed that the general character of 
the Barkhausen effect was about the same on the 
steep parts of all the curves. This result is in 
decided contradistinction the effect 
duced on the Barkhausen phenomenon by slight 
elastic bending of a wire. In the latter case the 


with pro- 


Barkhausen jumps are made larger and fewer, 
whereas pure compression appears to have no 
effect of this kind. 

When a stress of 12.1 kg/mm? was applied to 
the nickel the curve obtained was not so steep as 
the one for a stress of 10.2 kg/mm. It appears, 
therefore, that plastic flow must have occurred 
during the adjustment of the weights for the 
higher stress, and the nickel had become 
hardened. 

The presence of internal strains in the nickel 
initially will, of course, make the actual operative 
somewhat different from the 
listed in Fig. 2. The nickel, however, had been 


stresses values 
annealed and appeared to be quite soft, so 
probably it did not have large internal strains. 
If internal compressive strains existed they would 
have the effect of making the values listed in 
Fig. 2 smaller than the actual maximum of stress 
in the material. 

The experiments lead to the conclusion that a 
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Fic. 2. Descending branches of hysteresis loops for a nickel 
bar subjected to various compressive stresses. 


rectangular hysteresis loop in nickel cannot be 
obtained by compression alone, although an 
approximation to that shape is secured by using 
sufficiently high pressures. It is now of some 
interest to determine whether compressive 
stresses in bent nickel wires need to be higher 
than those used in this experiment with the bar 
before the large discontinuity of Fig. 1 appears. 

In the case of a wire elastically bent into the 
arc of a circle the compressive stress may be 
calculated from the equation ¢,=Ex/p. Here E 
is Young’s modulus, p is the radius of curvature 
of the arc into which the wire is bent, and ¢, 
is the compressive stress at a point distant x 
from the neutral surface where the strain is zero. 
The maximum compressive stress, o, is obtained 
when x=r, the radius of the wire. 

A wire of radius 0.006 cm was drawn down 
from a piece of the nickel bar and heated to 
redness to remove internal strains. A straight 
piece of this wire was then thrust into a capillary 
tube curved to make p= 23 cm. Taking E=2X 10+ 
kg/mm?’ we get o,, = 5.2 kg/mm?. The Barkhausen 
effect in this wire was investigated using a small 
coil, amplifier, and telephone receiver. Outside 
the tube when unstrained, the wire showed the 
usual large number of very small Barkhausen 
jumps. Inside the tube no jumps were observed 
until the demagnetizing field reached a suffi- 
ciently large value, when a single large discon- 
tinuity occurred. This large discontinuity was 
followed by a number of small ones. It appears, 


HEAPS 


therefore, that large stresses are not necessary 
in a bent wire for the production of the large 
discontinuity. These discontinuities will appear 
even though the maximum compressive stress in 
the bent wire is considerably smaller than some 
of the compressive stresses recorded in Fig. 2, 
where no discontinuities are to be observed. 

In explaining the jump BC of Fig. 1, therefore, 
it seems necessary to suppose that the stretched 
portion of the wire plays an important role in the 
production of instability. A somewhat simpler 
distribution of strains than is found in elastically 
bent wires may be obtained by the following 
device. A nickel bar is fitted into a nickel tube 
which has plugs brazed in at each end. The bar 
is inserted through a slot in the side of the tube 
and is strongly compressed by means of a heavy 
bolt threaded through one of the plugs. In this 
way the tube is strongly stretched and the bar 
strongly compressed. Such an assembly shows 
a Barkhausen effect different in no noticeable 
way from that observed when the strains are 
removed by unscrewing the bolt. A continuous 
lateral, mechanical connection between the 
stretched and compressed portions of the 
specimen appears necessary if the discontinuity 
is to be produced. 

The main features involved in the production 
of a discontinuity are probably as follows. As the 
field is decreased from A towards B in Fig. 1 the 
compressed portion of the wire decreases its 
intensity of magnetization by a very small 
amount, the stretched portion by a much larger 
amount. At B, however, a certain small region 
of the compressed portion where strain is a 
maximum, will have reached and passed the knee 
of the loop. The magnetization of this small 
region will be zero although the rest of the wire 
still retains residual magnetism. A slight increase 
in the demagnetizing field will now cause this 
small region to acquire an increment of mag- 
netization in the direction of the field. Whereas, 
previously, during the demagnetization process 
it had been expanding, it now begins to contract. 
This contraction will reduce the strain of the 
contiguous stretched region and allow a decrease 
of residual magnetism here. The net result is 
that the small region acts as a nucleus of 
instability and a wave of reversal of magnetism 
travels along the wire in both directions out 
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from the nucleus. This wave of magnetic 
reversal will be accompanied by a wave of 


magnetostrictive strain, of eddy currents, and of 
temperature rise, so that the interaction of the 
various factors is decidedly complex. 

In many cases the waves run completely to 
the ends of the wire; in other cases a magnetic 
wave may be blocked before it reaches the end. 
This blocking, or freezing, of the wave is probably 
due to unfavorable local strains at the point of 


‘Raymond E. Reinhart, Phys. Rev. 46, 483 (1934). 
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blocking. When it occurs a magnetometer as 
ordinarily of an 
accurate picture of the magnetic state of the 


used cannot, course, give 
wire. Some device must be used to detect the 
existence of consequent poles,—usually moving 
the wire parallel to its length past the mag- 
netometer needle serves to indicate the presence 
of such anomalies. 

The writer is indebted to Mr. H. B. Slayback 
of The Metal Goods Corporation and to the 
Driver Harris Company for supplying the nickel 
used in these experiments. 
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Measurements were made in the temperature range 
from —185°C to +100°C on bismuth crystals containing 
Pb, Sn, Sb and Te. The effects were found to be compli- 
cated but the following generalizations may be made. 
Decreasing temperature results in greater effectiveness of 
the impurity present except when the impurity concen- 
tration is such that a separation of phases occurs at low 
temperatures. Extremely small amounts (less than 0.03 
percent) of all impurities cause a sharp increase in the 


resistance. The effect of larger amounts depends on the 


INTRODUCTION 
HIS study of the effect of impurities upon 
the electrical resistivity of single crystals 
of bismuth, was originally undertaken at the 
California Institute of Technology as a compan- 
ion problem to the magnetic study reported by 
Goetz and Focke.! The data obtained at that 
time have not been used in the present paper, 
although the results obtained here are in very 
good agreement with the previous ones. It was 
considered much more satisfactory to complete 
the entire study using but one supply of bismuth. 
In this way the undetermined quantity of 
residual impurities would be the same for all 
crystals, so that direct comparisons, could be 
made. 
Previous measurements of the electrical re- 
sistivity of bismuth single crystals by Kapitza,? 


! A. Goetz and A. B. Focke, Phys. Rev. 45, 170 (1934). 
Kapitza, Proc. Roy. Soc. A119, 358 (1928). 
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nature of the impurity. Pb and Sn continue to increase the 
resistance, additional Te forces the resistance to fall to 
lower values than that of pure Bi, additional Sb causes the 
initial rise of resistance to fall off and then to be slowly 
reestablished. More than 0.03 percent Sn and more than 
0.3 percent Pb result in a negative temperature coefficient 
for the resistance parallel to the principal axis. In other 
cases the coefficient remains positive. The results may be 
explained on a combined theoretical and empirical basis. 


Bridgman*® and Goetz and Focke,* have clearly 
demonstrated the variations which may occur 
in different samples of ‘chemically pure’’ bis- 
muth. These variations suggested the desirability 
of a quantitative investigation of the electrical 
resistivity of bismuth as a function of the 
amount and kind of impurity present. 

In the magnetic problem mentioned above, the 
impurities which were found to cause the prin- 
cipal changes were lead, tin and tellurium. These 
metals were therefore the ones which have been 
used in the present work. Antimony has also 
been used as an impurity because it is an iso- 
morph of bismuth and might be expected to act 
as a control. The case of tin is especially im- 
portant because of the work of Ufford® on the 
resistance of polycrystalline alloys of the tin- 


3 P. W. Bridgman, Proc. Am. Acad. Sci. 63, 351 (1928). 
* A. Goetz and A. B. Focke, Phys. Rev. 37, 1044 (1931). 
>C. W. Ufford, Proc. Am. Acad. Sci. 63, 309 (1928). 
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bismuth series. Jones® has used the results 
obtained by Ufford as a basis for a theoretical 
discussion of both the electric and magnetic prop- 
erties of bismuth. Thomas and Evans’ made 
measurements on the temperature coefficient of 
resistance of polycrystalline lead-bismuth alloys 
within the range from —20°C to +30°C. 


EXPERIMENTAL PROCEDURE 

The bismuth used was specified as ‘Merk, 
Pure 99.5 percent to 100 percent.’’ No chemical 
tests were made to determine its absolute com- 
position because preliminary tests indicated that 
the probable maximum total impurity was not 
greater than 0.01 percent to 0.02 percent. The 
results obtained were in agreement with this 
allowance. The metals used as impurities were 
‘chemically pure.’’ Any minute amounts of con- 
tamination in them would be entirely negligible 
because of the small quantities used. 

The masses of metals added to the bismuth 
were such as to result in crystals containing 
atomic concentrations of 0.01, 0.03, 0.09, 0.27 
and 0.81 percent of each impurity studied. Two 
crystals were grown from each concentration. 
One of these had the principal axis [111] parallel 
to the length of the rod, while the orientation of 
the other was 90° from this. In the notation 
introduced by Goetz,* these orientations are 
designated as P; and P;, respectively. 

Ten grams of bismuth was heated in a Pyrex 
crucible to a temperature of 500°C under an 
atmosphere of hydrogen, so that the oxide 
could be reduced. The bismuth was then cooled 
to about 350°C and the impurity was added. 
The mixture was thoroughly stirred and poured 
into a small steel mould which was heated to 
about 270°C and was also in a hydrogen atmos- 
phere. Thus bars of all alloys were formed which 
had dimensions of 2X2 X150 millimeters. The 
bars were cut in half, each half later being formed 
into a single crystal. Magnetic tests showed no 
evidence of iron in the resulting samples, which 
is sufficient proof that no contamination of the 
alloys was caused by contact with the mould. 

The single crystals were grown by progressive 
crystallization from the molten metal, the orien- 





®°H. Jones, Proc. Roy. Soc. A147, 396 (1934). 

7W. R. Thomas and E. J. Evans, Phil. Mag. 16, 329 
(1933). 

8 A. Goetz, Phys. Rev. 35, 193 (1930). 


tation being predetermined by a seed crystal. 
The procedure followed closely that described by 
Goetz and Focke! with the one difference that 
here the furnace was moved while the crystal 
remained at rest. Much of the trouble due to 
vibrations during crystallization was thus elimi- 
nated, and single crystals were obtained with 
practically no failures. 

The electrical measurements were made by the 
potentiometer method of comparing the poten- 
tial difference observed across a standard resist- 
ance, with that observed across the sample 
under test when the same current was passed 
through both. 

Measurements were made at the temperatures 
of boiling liquid air (—185°C), carbon dioxide 
snow and acetone (—77°C), ice and water 
(0°C), water at room temperature (22°C) and 
boiling water (100°C). In all cases the crystals 
were placed directly into the temperature baths. 

The results obtained in this manner were 
satisfactory with the exception of two crystals 
containing tin as the impurity. For these a 
nearly continuous change was desired between 
about — 150°C and room temperature. This was 
obtained by inserting the crystal being studied, 
into a hole drilled in a massive cylinder of brass. 
The brass was cooled with liquid air and covered 
with a thermos flask so that its temperature 
would increase very slowly. A_ thermocouple 
placed beside the crystal was used to determine 
its temperature. By this means it was possible to 
obtain measurements at 5° intervals from 
—150°C to +20°C, the entire change requiring 
about two hours. 

RESULTS 

A summary of the results is given in Table I. 

The group of isothermal curves showing the 
dependence of the resistivity parallel to the 
principal axis (111), upon the amount of impurity 
present for crystals containing tin and lead is 
given in Fig. 1. It is evident from these curves 
that the effects produced by tin and lead are 
qualitatively the same, the principal difference 
being in the relative effectiveness of the different 
atoms. One atom of tin produces the effect of 
from three to five atoms of lead. Since this 
similarity occurs for both P; and P; orientations, 
further curves and discussion will treat only the 
case of tin. 
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TABLE I. Electrical resistivity of bismuth single crystals. 
All resistivities multiplied by 10+. 
Resistivity parallel to the principal axis 
Pure BIsMUTH 
— 185 —77 0 22 100° 
0.36 0.88 1.20 1.28 1.82 
CRYSTALS CONTAINING ANTIMONY 
0.016 0.57 0.86 1.19 1.28 1.67 
0.03% 0.68 0.83 1.17 1.28 1.64 
0.09°, 0.57 0.83 1.26 1.36 1.77 
0.27% 0.45 0.86 1.22 1.35 1.78 
0.81% 0.47 0.86 1.23 1.35 1.76 
CRYSTALS CONTAINING LEAD 
0.01% 0.77 0.98 1.32 1.41 1.88 
0.03% 1.28 1.17 1.41 1.64 1.96 
0.09%; 4.05 2.15 1.93 1.97 2.18 
0.27% 3.37 3.24 2.57 2.48 2.52 
0.81% 1.79 5.25 4.10 3.62 3.14 
CRYSTALS CONTAINING TIN 
0.01% 0.88 0.91 1.28 1.50 1.82 
0.03% 3.14 2.10 1.96 1.94 2.10 
0.09°, 1.98 4.27 3.07 2.88 2.62 
0.27% 1.38 5.40 5.40 4.84 3.64 
0.81 1.21 5.70 7.25 6.70 5.42 
CRYSTALS CONTAINING TELLURIUM 
0.01; 0.42 0.74 1.16 1.26 1.68 
0.03% 0.27 0.55 0.98 1.05 1.52 
0.09°, 0.25 0.53 0.82 0.92 1.40 
0.27% 0.29 0.55 0.82 0.91 1.43 
0.81% 0.53 0.94 1.22 1.36 1.66 
Resistivity perpendicular to the principal axts 
PuRE BisMUTH 
—185 —77 0° 22 100° 
0.33 0.74 0.98 1.06 1.36 
CRYSTALS CONTAINING ANTIMONY 
0.01% 0.53 0.78 1.06 1.16 1.47 
0.03% 0.44 0.70 0.98 1.06 1.41 
0.09%, 0.41 0.66 0.96 1.05 1.40 
0.27% 0.42 0.69 0.97 1.05 1.37 
0.81% 0.48 0.76 1.03 1.10 1.44 
CRYSTALS CONTAINING LEAD 
0.01% 0.56 0.76 1.00 1.08 1.40 
0.036; 0.72 0.86 1.06 1.14 1.44 
0.09°; 0.94 1.05 1.16 1.22 1.44 
0.27% 0.91 1.27 1.37 1.40 1.69 
0.815; 0.82 1.77 1.97 2.00 2.07 
CRYSTALS CONTAINING TIN 
0.01% 0.86 0.98 1.08 1.14 1.40 
0.03% 0.97 1.10 1.23 1.27 1.58 
0.09% 0.76 1.44 1.57 1.60 1.75 
0.27% 0.56 1.68 1.99 2.05 2.11 
0.81 0.85 1.80 2.36 2.43 2.51 
CRYSTALS CONTAINING TELLURIUM 
0.01% 0.48 0.72 0.96 1.04 1.37 
0.03; 0.28 0.62 0.91 0.98 1.32 
0.09°; 0.24 0.51 0.78 0.88 1.28 
0.27% 0.25 0.52 0.76 0.85 1.22 
0.81% 0.38 0.69 0.95 1.04 1.41 


Three important characteristics can be seen in 
Fig. 1: first, the tremendous effect produced by 
extremely small amounts of the impurities; 
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second, the marked difference between the meas- 
urements at liquid air temperature and for those 
at higher temperatures; and third, the existence 
of a negative temperature coefficient of resistance 
within certain temperature ranges in the cases of 
some of the crystals. 

These abnormal temperature coefficients are 
more clearly shown in Fig. 2, in which the 
resistivity is plotted as a function of temperature 
for P; crystals containing tin. In order to deter- 
mine the true shape of these curves, measure- 
ments were taken at intervals of 5°C in the range 
—150°C to +20°C in the cases of the 
crystals containing 0.09 percent and 0.27 percent 


from 


tin. 

The isothermal curves for P3 crystals con- 
taining tellurium and antimony are given in 
Fig. 3. It may be seen from these that the effects 
produced by tellurium and antimony are very 
much smaller than those produced by tin and 
lead. Furthermore, while tin and lead increased 
the resistivity, tellurium tends in general to 
decrease it. Antimony seems to be a hybrid, 
which decreases the resistivity at high tempera- 
tures and increases it at low temperatures. For 
amounts greater than 0.09 percent antimony the 
changes produced are negligibly small. 

Fig. 4 shows the variations produced in the 
resistivity perpendicular to the principal axis. 
In the case of tin the difference in behavior at 
—185°C is again apparent; otherwise the only 
irregularity is the initial increase in the resistivity 
observed both in the case of tellurium and of 
antimony. 

The results given in Figs. 3 and 4 indicate that 
there is no great advantage to be gained from 
plotting these resistivities as a function of tem- 
perature. The temperature coefficients will be 
positive at all times although they would ob- 
viously be greatly increased between —185°C 
and 0°C for those crystals containing more than 
0.03 percent tin. 


DISCUSSION 


A nearly complete explanation of the results 
may be obtained through a consideration of the 
theoretical work of Jones® together with the 
phase equilibrium diagram of the bismuth-tin 
system. 
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Fic. 1. Resistivity of bismuth crystals, parallel to the principal axis, as a_function of the amount of tin and 
lead present. 


Jones has attacked the problem of the electric 
and magnetic properties of bismuth by making 
use of the Bloch theory and has determined that 
the number of electrons within the inner Bril- 
louin zone in the bismuth lattice is very nearly 
equal to five. Since the number of valence elec- 
trons for bismuth is five, this leaves, on the 
average, only a very small fraction of one electron 
per atom in the outer zone which is therefore 
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Fic. 2. Resistivity of bismuth crystals, parallel to the 
principal axis, as a function of temperature for bismuth 
crystals containing tin. 


available for participation in normal metallic 
conduction. 

If some of the bismuth atoms are replaced by 
foreign atoms having a valence of four, the aver- 
age number of these electrons per atom, for the 
whole lattice, will be reduced below five. If this 
is carried far enough there will be no electrons 
remaining in the outer zone, and ordinary metallic 
conduction will be replaced by one which is 
characterized by a negative temperature coef- 
ficient. In this way Jones has explained the 
results of Thomas and Evans’ and Ufford.® 

This same reasoning would explain the initial 
effects observed for all temperatures with either 
tin or lead as the foreign atoms. It would also 
explain the results obtained at temperatures 
above —10°C for all crystals containing these 
impurities. 

A possible explanation of the decided changes 
that occur at lower temperatures may be ob- 
tained from the phase equilibrium diagram for 
the bismuth-tin system. This is similar to the 
bismuth-lead system and is given in Fig. 5. The 
important thing to notice here is that the 
boundary between the solid solution 
phase A, and the mixture of solid solutions 
A+B, is not a vertical line but approaches the 
pure bismuth axis as the temperature decreases. 
This indicates that a single crystal of solid solu- 


single 
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tion A is not stable at all temperatures, but that 
if the temperature is sufficiently reduced, a 
recrystallization must take place so that some 
of the tin, together with some of the bismuth, 
must leave the bismuth lattice. 

If this occurs, the average number of valence 
electrons per atom remaining in the bismuth 
lattice must be increased, thus forcing the resist- 
ance back toward that for the pure metal. Since 
this is just the type of behavior that is observed, 
this is probably a valid explanation of the results. 

It is interesting to note here that this is also 
a probable explanation of the similar peaks 
observed in the magnetic study.! 

The curves for antimony given in Figs. 3 and 
4 cannot be explained on the basis of Jones’ 
theory, because antimony possesses five valence 
electrons. Thus there can be no change made in 
the average number of electrons per atom by the 
addition of antimony. It may be that the presence 
of the antimony atoms increases or decreases the 
existing mechanical strains within the crystal, 
the type of change depending upon the tem- 
perature and the orientation in the crystal. 

This type of distortion would probably be 
present regardless of the kind of impurity added 
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Fic. 3. Resistivity of bismuth crystals, parallel to the 
principal axis, as a function of the amount of tellurium 
and antimony present. 
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Fic. 4. Resistivity of bismuth crystals, perpendicular to 
the principal axis, as a function of the amount of impurity 


present. 
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Fic. 5. The phase equilibrium diagram of the bismuth-tin 


system. (After Endo.) 
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Fic. 6. The resistivity of bismuth crystals containing 
antimony subtracted from the resistivity of crystals 
containing tin and tellurium in an attempt to eliminate 
mechanical strain effects. 
to the bismuth, so that in general the observed 
effect would be a combination of the mechanical 
and electronic disturbances in the lattice. 

If the mechanical distortions due to atoms of 
nearly the same atomic weight can be assumed 
to be approximately the same, the pure electronic 
effect would be obtained by subtracting the ob- 
served antimony result from that of the corre- 
sponding quantity of tin or tellurium. In the 
case of tin such a procedure would not cause any 
particular change in the character of the curves 
obtained, but in the case of tellurium the com- 
plicated curves given in Figs. 3 and 4 would be 
replaced by those given in Fig. 6. These curves 
show the same general characteristics as the 
curves for tin at —185°C except that here the 


resistivity is decreased instead of being increased. 
This may be explained by the same reasoning as 


AND J. 
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that given for the case of tin when the fact that 
tellurium has a valence of six instead of four is 
considered. The tellurium would 
increase the average number of electrons per 
atom thus decreasing the resistivity of the 
crystal containing it. The sharp breaks in the 


presence of 


curves would again be explained by recrystal- 
lization. 
At first 
taking place with such small amounts of impurity 
seems rather questionable. The work of one of 
the authors,® however, has definitely shown that 
in the case of polonium, the foreign atoms are 


the possibility of recrystallization 


segregated into groups and are not distributed at 
random throughout the bismuth crystals. The 
recrystallization would thus take place in the 
regions of high impurity and there would be no 
question of the need of migration of foreign 
atoms through the solid, to permit the formation 
of the new lattice. 

The results may be summarized as follows: 
The addition of tin or lead to bismuth causes a 
marked increase in the electrical resistivity. This 
increase is greater, the larger the amounts of 
impurity added, for all temperatures, so long as 
no recrystallization takes place: The addition of 
tellurium usually causes a decrease in the re- 
sistivity so long as no recrystallization occurs: 
The addition of antimony causes an increase in 
the resistivity except in the instance of the re- 
sistivity parallel to the principal axis [111] at 
temperatures of the order of 100°C. 

Approximate limits of solubility of tin, lead 
and tellurium in solid bismuth may be estimated 
as follows: tin, 0.03 percent at —185°C, 0.09 
percent at —95°C, 0.27 percent at —42°C and 


0.81 percent at —20°C; lead, 0.09 percent at 
—185°C, 0.81 percent between —185°C and 
—77°C; tellurium, 0.03 percent at —185°C, 


0.12 percent at +100°C. 

The results may be explained quite well by 
the theoretical work of Jones® together with the 
assumption that the effects due to antimony 
give an indication of the changes produced by 
mechanical factors alone. 

We wish to express our thanks to Professor 
Goetz of the California Institute of Technology 
for his suggestion of the problem and his help 
and interest in the early part of the work. 


, A. B. Focke, Phys. Rev. 46, 623 (1934). 
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Use of the Reflection Echelon for Interferometric Wave- 
Length Comparisons Extending into the 
Schumann Region 


Interference patterns have been photographed which are 
] gral 


suitable for the comparison of wave-lengths of lines in the 


spectrum of Cull from 1488A to 1621A with the wave- 
2100A) of 


Burns and Walters.' A reflection echelon 


lengths (greater than lines measured inter 
ferometrically by 
was used with a concave grating as shown in the plan in 
Fig. 1. 
chamber. Eight wave-lengths shorter than 2100A computed 
Cu Il 


compared with values deduced from the patterns 


This optical system was enclosed in a vacuum 


been 
The 


greatest discrepancy was 0.004A, and five of the values 


by Shenstone? from term values for have 


differed by less than 0.001A from the computed values. 
hese wave-length values were based on the assumption 
t of a zero average error in Burns and Walters’ values. The 
precision of measurement was low because of the difficulties 
to be discussed below. Consequently no wave-length values 


based on these patterns will be published. 


The only source found to radiate with sufficient intensity 
and homogeneity was the hollow cathode copper discharge 
in helium described by Shenstone.? The radiation entered 
the chamber through the slit, S. After collimation by the 
platinized concave mirror, C, the beam fell on the twenty 
five plate, platinized, Hilger fused quartz reflection echelon, 
i, placed with edges of the steps vertical. The reflected 
beam fell on the two meter, platinized, concave, glass 
grating, G, having 30,000 lines per inch and placed with the 
rulings vertical. The resulting parallel dispersions were 
found to be advantageous in permitting long fringes 
without increase of overlap of adjacent patterns. The 
grating was rotated about a vertical axis to place various 
wave-length regions on the plate, P, whose distance from 
beam from the grating 


the grating and inclination to the 


could be adjusted for focus. Ilford Q plates five inches long 


were used, recording simultaneously all of the patterns 


within any 1000A wave-length range. The beam incident 
on the echelon was exactly perpendicular to the edges of 
the echelon steps, a necessary condition for reliable wave- 
rhe 
determinations in 


Williams. 


dispersions, in which case the edges of the echelon steps 


length comparisons. fiducial mirrors were used fo 


wave-length essentially the manner 


described by some work was done with crossed 


were horizontal, and the center ol ¥ was be low the line 


rhe 


necessary for crossed dispersions can be secured only when 


joining the centers of & and G, stigmatic image 


P is on the normal to the grating. The dimensions of the 
vacuum chamber prevented the fulfillment of this condition 
for wave-lengths longer than 2400A, 

Ihe most serious disadvantage of the reflection echelon 
is the concentration of intensity near the single order fringe 
position. As a consequence there are never more than two 
measurable Iringes, and errors of measurement cannot be 
eliminated by averaging a number of fringes. This concen- 
tration of intensity also causes displacements of all but 
extremely sharp fringes toward the single order position by 
amounts proportional to their distance from the single 
order position and to their widths. The widths of the lines 
of the uncooled hollow cathode discharge were observed to 
be 0.01A, which caused displacements amounting to about 
twenty percent of the distances of the fringes from the 
single order position. The extremely high order of inter 
100,000 at 1460A) in 


makes the measurement of the patterns and the identifi- 


ference the extreme ultraviolet 
cation of the order of interference increasingly difficult for 
shorter wave-lengths. 

It can be concluded that the reflection echelon is suitable 
for wave-length comparisons in the Schumann region when 
the optical system described is used. This outline is of work 
fully described in a thesis submitted to the faculty of the 
Massachusetts Institute of Technology and performed 
under the direction of Professor G. R. Harrison to whom 


grateful acknowledgment is now made. 


Davip L. Mac Apa 

Massachusetts Institute of Technology, 

June 18, 1936. 

Burns and Walters, Publ. Alleg. Obs. 8, 27 (1930 

\. E. Shenstone The First Spark Spectr ( er P 
rr Ns Is Lan n 235, 145 (1936 

\\ ! W s ‘ literter Il. ¢ 
iz 1 ul Use t Retles I elonus tor e \ ‘ l 
Regions,’ Proc. Phys. Soc. 45, 699 (1933 
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FIG. 1. Mass sx 


The Isotopic Constitution of Strontium and Tellurium 


In the Physical Review for May 15, J. P. Blewett and M. 
3. Sampson report the observation of a new isotope at mass 
84 in strontium. I have recently confirmed this observation, 
the mass spectrum of the ions from a spark between a pal- 
ladium and a strontium electrode showing the mass at 8&4 
in addition to the much stronger lines at 86, 87, and 8&8. 

Mass spectra of the charged atoms from a spark between 
tellurium and palladium electrodes showed a new faint 
isotope at mass 120. A reproduction is shown in Fig. 1. 
Only the four heaviest isotopes at 130, 128, 126, 125 were 
observed by Aston. Later Bainbridge! found three addi- 
tional isotopes at 124, 123, and 122, and indications of an 
extremely faint one at 127. The present photographs show 
no trace of this line at 127, and it may have been due to 
iodine present as an impurity. The new isotope at 120 can- 
not be ascribed to tin as the other strong tin isotopes are 
not present. In a chart of the elements it extends the limit 
of stability of the tellurium nucleus to make it more in line 
with the lightest isotopes of tin, xenon, and barium. 

A. J. DEMPSTER 


University of Chicago, 
June 23, 1936. 


1K. T. Bainbridge, Phys. Rev. 39, 1021 (1932 


Evidence from Efficiency Curves for the Nature of the 
Disintegration Process for Boron 


The existence of a continuous distribution and a homo- 
geneous group of a-particles resulting from bombarding 
boron with protons has been definitely established,': * and 
confirmed by integrated numbers-vs.-range curves in this 
laboratory. 

This letter reports an investigation of the efficiency of 
production of a-particles from this reaction as a function of 
the voltage accelerating the protons. A transformer-kene- 
tron-condenser source of high potential up to 275,000 volts 
with voltage measured by a high resistance voltmeter was 
used. The energy of the magnetically resolved proton beam 
could be controlled to within three kilovolts and except for 
collision energy losses after resolution was highly monochro- 
matic. The a@-particles emitted from the fused B,O, thick 
target were observed by an ionization chamber and amplifier 
of the Dunning* type and a thyratron counter adjusted to 
count all @-particles reaching the ionization chamber. The 
apparatus gave absolute yields from a thick target of LiCl 
to within 10 percent of the results of the Wisconsin 
laboratory.4 

The observed efficiency curves for production of all 


e-particles having a range greatel than a given range are 


122 124 126 128 
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shown in Fig. 1. It is seen that the main group of range 
greater than 18 mm follows the ordinary efficiency curve 
with no evidence for any peculiarities of excitation. This is 
also true for the Li efficiency curve shown by the dotted 
line. However, when the majority of a-particles recorded 
are of the long range tvpe, >38 mm, it is seen that the curve 
is not exponential but shows a decided break with the 
maximum of the differentiated curve (i.e., thin target 
180 kv. The feet of these curves are of a quite different 
character from those of the shorter 1 inge curves and the 
intermediate sections show a more rapid rise in vield with 
increasing voltage. 

rhese curves suggest that the disintegration which has 
been proposed 

i1H'+;B'"'—,Be’+.Het 


or any further reaction postulated to account for this long 
range group, is a case of resonance disintegration. 
Reasons for the non-constancy of the vield from thick 
targets at higher voltages are: (a) the onset of further modes 
of disintegration; (b) the change of shape of the numbers 
vs.-range curve so that the continuous distribution extends 
to greater range and is observed with increasing abundance 
at higher voltages. From this point of view the short range 
continuous distribution efficiency curve is not “regular 
rhere is evidence for the latter suggestion from the differ 


ential curve, Fig. 5 of Cockcroft and Lewis.” 
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An attempt to overcome the latter effect by working near 
the end of the homogeneous group range, absorber more 
than 45 mm, met with inherent experimental difficulties. 

These observations are substantiated by the observations 
of Cockcroft and Lewis? that: “‘the ratios of the numbers 
counted at the peak of the main curve at 2 cm to the num- 
bers at the peak of the end group at 4 cm changes from 
22: 1 at 185 kv to 90: 1 at 360 kv.” 

An explanation of the observed efficiency curve for the 
long range homogeneous group is dependent on a satisfac- 
tory explanation of the proton-boron disintegration prod- 
ucts, for which there have been several mechanisms 
devised.* 

Qualitative results for the yield of a-particles per proton 
calculated for a pure boron target by multiplying the 
observations by 34/10 at 150, 175, 200 and 225 kv are: 
for range greater than 18 mm, 1.68X10~-°, 7.06X10~°, 
18.4X10-°, and 32.3 x 10~°; for range greater than 40 mm, 
1.79X10-", 17.410", 49.7 10-" and 59.7 X 107". 

We wish to express our sincere thanks to Professor E. L. 
Hill for many stimulating discussions and to Messrs. W. G. 
Shepherd and R. O. Haxby for invaluable assistance with 
the recording apparatus. This research was supported by a 
grant from the Graduate School to Professor John T. Tate 
and the senior author. 


University of Minnesota, 
Minneapolis, Minnesota, 
June 19, 1936 


Joun H. WILLIAMS 
Wittiam H. WELLS 


1933); Kirchner and Neuert, Physik 


1 Kirchner, Naturwiss. 21, 473 
Roy 


Zeits. 34, 897 (1933); Oliphant, Kempton and Rutherford, Proc. 
Soc. A150, 241 (1935); Neuert, Physik Zeits. 36, 629 (1935) 
? Cockcroft and Lewis, Proc. Roy. Soc. A154, 246 (1936 


3 Dunning, Rev. Sci. Inst. 5, 387 (1934). 
‘ Herb, Parkinson and Kerst, Phys. Rev. 48, 118 (1935); Hevdenberg, 
Zahn and King, Phys. Rev. 49, 100 (1936) 


5 Dee and Gilbert, Proc. Roy. Soc. A154, 279 (1936). 


Photon Theory and Compton Effect 


Recently Shankland! reported coincidence experiments 
on the Compton effect with Ra C y-rays, which led him to 
the conclusion that the photon theory fails to explain such 
experiments. On the other hand Bothe? and Geiger had 
shown some years ago, that the Compton effect is in full 
agreement with this theory for x-rays of about 70 kv. 
This result has been carefully established and is to be 
considered as correct, although meanwhile there has been 
some progress in the experimental technique. The result of 
Shankland has been discussed theoretically by several 
authors.® 

We have reexamined Shankland’s result, to find out 
whether really the behavior of radiation changes so funda- 
mentally if the frequency is altered by one order of magni- 
tude. Our experiments do not confirm Shankland’s results. 
We found the photon theory able to describe the scattering 
of hard y-rays as well as of the x-rays. 

The arrangement was similar to the one used by Shank- 
land except some simplifications. As primary radiation we 
used the filtered y-rays of Th C” (20 mg RaTh), which 
have the advantage of being more homogeneous than the 
y-rays of Ra C. The scatterer was a Cellophane foil 0.028 
mm thick. The recoil electrons were counted by one thin- 
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TABLE |. Summar f result 
| Coincidences in 14.5 Hours 
= . » , 1 
Exper Scattering Recoil 
No. and Recoil | Electrons 
Angle per min Without Wit! Surplus 
Foil Foil ' 
1* 30 40 5 25 20 
: 30 52 34 95 61 
3 30 40 30 39 9 
21 40 416 7 
4 30 38 39 95 s6 
45 28 50 11 
* Only one tube counter for scattered y-rays 
walled Geiger-Miiller counter 1.2 cm in diameter, the 


scattered photons in most cases by two connected lead 
counters of the same size. The 8 and y counters were 
mounted symmetrically on opposite sides of the primary 
beam at a distance of about 2.1 cm from its axis; only for 
the experiment number 3 the azimuth was not 180°, but 
90° (designated by j). The foil could be moved in the 
direction of the primary beam. According to the photon 
theory coincidences were to be expected if the foil was 
placed in such a manner that the direction from the foil 
to both the recoil electron counter and the photon counters 
formed an angle of 30° with the primary beam. 

The results of experiments are given in Table I. The 
second column gives the angle between the primary beam 
and the line from the foil to each one of the counters. The 
third column gives the increase of 8 counts caused by the 
foil. The fourth column shows the zero effect; to obtain it 
the foil was removed and the counting rates of the 8 
counter were, by means of a weak source of 8-rays, adjusted 
to the same amount as found when the foil was present. 
The y counting rates were not changed noticeably by the 
foil. In the fifth the 
counted with the foil are shown, in column 6 the difference 
caused by the foil. It is seen that in the position predicted 
by theory the number of coincidences is far greater than 
the zero effect, while in every adjacent position there is 


column numbers of coincidences 


only a small fraction of this effect. 

These experiments show, firstly that scattered photons 
and recoil electrons are emitted simultaneously, and 
secondly that their directions are 
predicted by Compton's theory. By a 
discussion, which will be published elsewhere, it is seen that 


with 


interconnected as is 
more detailed 
the observed number of coincidences is consistent 
the assumption of this connection holding quantitatively. 

To get a Shankland's 
negative result it would be necessary to consider some more 


satisfactory explanation of 
experimental data which were not contained in his paper. 
We believe that the inhomogeneity of the Ra C y-rays 
played an important role in his experiment. 

BOTHE 
MaAteER-LEIBNITZ 


Institut fiir Physik W. 
am Kaiser Wilhelm—lInstitut 

fiir medizinische Forschung, H. 
Heidelberg, June 6, 1936. 


' R. S. Shankland, Phys. Rev. 49, 8 (1936). 

?W. Bothe and H. Geiger, Zeits. f. Physik 32, 639 (1925) 

3 P. A. M. Dirac, Nature 137, 298 (1936); E. J. Williams, Nature 137, 
614 (1936). R. Peieris, Nature 137, 904 (1936) 

‘The experiments 1, 2, 3 have been published in Géttinger Nach- 
richten, No. 2, 127 (1936). 
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Band Spectrum of OH* 
Loomis and Brandt! have found that the constant A of 


the multiplet splitting varies from —66 to —97 cm~. This 
according to the authors is accounted for by the uncer- 
tainty of the energy formulae. Indeed Loomis and Brandt 
used the energy formulae of the molecule in the state which 
corresponds to the pure case a. Meanwhile the molecule 
OH* is in the state which is intermediate between a and } 
cases even when J is small. 

Hill and Van Vleck? gave an exact equation of the third 
which holds 
We have 
solved this equation using the Cardan formula.’ For 
K =10 and 11 and when a=4/B, has values —4, —5, —6, 
—7, —8, we get the solutions summarized in Table I. In 
Table II the IIo — II, 
‘II, —*Il, are given. These differences are plotted against a 


degree (for the computation of triplet terms) 
true for the limiting and intermediate cases. 


corresponding differences and 
in Fig. 1. On the other hand these differences for K = 10 and 
11 can be determined by using in the following formulae the 


experimental data given by Loomis and Brandt. 


371,(10) —11, (10) 
=[?Qi2(11) — P2(11)+2R12(10) —Q2(10) ]/2 
= — 22.095 cm™ (for the (1, 0) band). 


Similarly, 


s11.(11) —11, (11) 
=[?Q12(12) — P2(12) +¢@Rio(11) —Q2(11) J/2 
= — 21.475 cm™ (for the (1, 0) band). 


According to Loomis and Brandt B,’=12.381 cm™. 
TABLE I. 
K 
10 11 
a 
¥ +1 4 0 ) —1 +1 . 0 , -1 
J=11 J =10 J=9 J=12 J=11 J=10 
‘Ile Ty Il Ile Thy TI 
B B B B B B 
—4 108.3323 | 109.1537 | 109.6343 130.3983 | 131.1295 | 131.5830 
—5 108.0363 | 109.2036 | 109.9233 130.1166 | 131.1722 | 131.8570 
—6 107.7003 | 109.2577 | 110.2730 | 129.8066 | 131.2191 | 132.1693 
—7 107.3250 | 109.3155 | 110.0523 129.4600 | 131.2693 | 132.5190 
—8& 106.9043 | 109.3978 | 111.0566 129.0710 | 131.3414 | 132.8940 
z==+1 >=0 z —1 ==+1 =() z —1 
TABLE II. 
K 
10 11 
a 
Il; —*Ile Ilo —*Ily Il; —*Ile ‘Ilo —*Ihy 
B, B B B, 
} 0.8214 0.4806 0.7312 0.4535 
5 1.1673 0.7197 1.0556 0.6848 
—6 1.5574 1.0153 1.4125 0.9502 
7 1.9905 1.3368 1.8093 1.2497 
-8 2.4935 1.6588 2.2704 1.5526 
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Fic. 1. Term differences plotted against a 


Therefore: 


11.(10) —*11,(10) = —1.8330 B,’ cm 
To(11) —*11,(11) = —1.7345 B,’ cm™. 


and 


Plotting these values on our diagram and dropping the 
perpendiculars on the a-axis we get the corresponding 
—6.82, A 
—82.2 cm™ and —84.1 cm”, respectively. 


values of a, namely a= —6.64 and a= being 
Let us now look at the differences “II, —“IIy. Averaging 
the a and + separations and using combination relations 
we get: 
M1, (10) —*119(10) =[05(10) — R3(10) +2052(10) — Q2(10 
+[P2(11)+R3(11) —*#Q32(11) —P3(11) ]/2 
1, 0) band). 


tN 


= — 17.135 cm™ (for the ( 


—7.18 and 


These values of a and A lie indeed more 


Returning to the diagram we find that a= 
A = —88.9 cm“. 
closely to the true ones than the values found by Loomis and 
Brandt. 

The similar computation of *II,(11)—*I)(11) is impos- 
sible because of lack of experimental data for certain 
frequencies when K = 12. The computation for lower values 
of K would enable us to get more precise values of A and a. 

The Hill and Van Vleck equation does not take into 
account the perturbations. As long as the formulae used in 
the calculation involve the perturbed branches Qs, P3 and 
R;, the discrepancy which still exists in the values should be 
accounted for by perturbations. 
to thank Nikolsky (Moscow 


University) for his consultation and help. 


I wish Professor K. V. 


M. 


KOVNER 


Woronezh University, U. S. R. S., 
May 16, 1936. 


Loomis and Brandt, Phys. Rev. 49, 55 (1936) 
? Hill and Van Vleck, Phys. Rev. 32, 259 (1928). 
Proc. State Optic. Inst. U. S. S. R. (in press). 
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Some Experiments with Photo- Neutrons 


We have performed experiments with the neutrons 
obtained from the photo-disintegration' of Be and D by 
the y-rays from radon. The y-ray source was about 700 mc 
of radon in a 2 mm diameter spherical bulb of Gundlach 
glass. The neutrons were detected with a boron lined 
ionization chamber previously described.? 

The source of photo-neutrons from Be was a rectangular 
parallelepiped of Be metal with the Rn bulb at its center 
(59 g of Be about 2.5X2.5X5.5 cm‘). This source was 
placed at the center of cylinders of water of various sizes 
and in each case the number of neutrons detected was 
counted with and without a } mm sheet of Cd in front of 
the ionization chamber. The growth in number of slow 
neutrons thus observed is shown in Fig. 1. In order to 
compare the slowing of these photo-neutrons with that of 
neutrons from a Kn—Be source of the sort used in earlier 
neutrons® the 
was obtained 


experiments on the production of slow 
corresponding growth curve, also in Fig. 1, 
under identical geometrical conditions. 

From the shape of the curves two conclusions may be 
drawn. The increasing portion of the curve indicates that 
the neutrons from the photo-disintegration of Be have less 
energy than those from the alpha-disintegration since they 
are slowed down to detectable speeds in a thinner layer of 
water. The relative mean free paths in water of neutrons 
from the two sources may be estimated from the decreasing 
parts of the curves since beyond the peaks of the curves 
there should be an approximate equilibrium between the 
fast neutrons and the slow neutrons which are strongly 
absorbed by Cd. The rate of decrease observed for the slow 
neutrons should then be the same as that of the fast 
neutrons. A comparison of the curves shows that the MFP 
of the Rn, — Be photo-neutrons is about 5 cm as compared 
with 9 cm for the Rng—Be neutrons. 

Chadwick and Goldhaber' using a radio-thorium source 
have measured the maximum energy of photo-neutrons 
thus obtained from Be. From their data they obtained 
1.6 MEV as the minimum y-ray energy required to release 


1.76 MEV 


y-lines' from Ra C should therefore give neutrons of 


the photo-neutron from Be. The 2.2 and 
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Fic. 1. Number of neutrons os. radius (R =ht.) of water cylinders 
The Rng — Be curves are marked @ and the Rn, — Be are marked y. In 
each case the curves 1, 2 and 3 are, respectively, all neutrons, slow and 
fast neutrons 
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8/9X0.16 and 8/9X0.6 MEV neutrons. As the relative 
efficiency of these two y-lines in producing photo-neutrons 
is not known we can only say that the above MFP of 5 cm 
is that of the mixture of these two groups of neutrons 

A rough estimate of the relative numbers of neutrons 
per mc of Rn is obtained from the maxima of the two 
curves. This indicates that } as many photo-neutrons as 
a-particle neutrons may be expected under these experi- 
mental conditions. 

Similar experiments were made with the Rn bulb in the 
center of a 200 cc sphere of pure D,O. In this case the 
vield of photo-neutrons is so low that the neutrons produced 
by the a-particle bombardment of the light elements in the 
be taken These 


into account 


glass of the bulb must 
“glass’’ neutrons were in our experiments } as numerous 
as the photo-neutrons from the deuterium. In order to 
separate the two effects in each case a count was taken 
with and without the D,O around the bulb of Rn. The 
count with the D.O minus that without the D,O (all other 
conditions being the same) gives the count from D,.O only 
The maximum of the growth curve shows that the number 


of photo-neutrons from the D.O was 2 of the 


percent 
number of neutrons obtainable from an equal amount of 
Rn mixed with Be powder. This in turn indicates that 
1 percent of the neutrons from the usual Rn—Be source 
come from the glass of the bulb. The peak of the growth 
curve for Rn—D photo-neutrons was at 7 cm (—3} for 
the space occupied by the D,O) in two experiments and 
8 cm (—3}) in a third. A growth curve for cme usual 
Rn,— Be source in the same D,.O and cylinders of water 
had its maximum at 10} cm (—3}). 

Using the net thickness of H,O through which the 
neutrons diffused when the number of detected neutrons 
was a maximum, we obtain: 

53 cm for Rn,—Be and 8} for Rng—Be under the 
same conditions, 
3} cm for Rn,—D,O and 7 cm for Rng— Be. 

This is at least a good indication that the photo-neutrons 
obtained from Be. 
If the binding energy of D is 2.14 MEV® and the hardest 
Ra C is 2.198 MEV* these photo- 


neutrons are of very low energy. That the effect with D is 


D are much slower than those from 


strong y-ray from 
observable suggests that there may be some y-radiation 
from Ra C that is harder than the above. 

We wish to acknowledge our debt to Dr. J. J. Duffy of 
Memorial Hospital for the loan of the Kn and to Professor 
H. C. Urey of Columbia University for the loan of the 


D.O used in these experiments. 
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Viscosity of Air and the Electronic Charge 


As one of the possible explanations for the discrepancy 
” and the ‘‘x-ray value”’ of the 


between the ‘‘oil drop value 
electronic charge it has been suggested! that Millikan has 
used too low a value for the coefficient of viscosity of air n, 
which enters as a factor of the power 3/2 in calculating e by 
the oil drop method. In order to examine this hypothesis, 
the author has undertaken a new determination of , using 
the rotating cylinder method. As a preliminary result was 


published? the value 


>= (18200+30)-10-8 corresponding to 


(18348 +30) -10-%. 


N23 

This is 0.67 percent higher than the value used by 
Millikan in 1916 (18226-10-%), and would give a value of 
e 1.0 percent higher than Millikan’s. 

Since this result was published, the apparatus has been 
improved so as to insure greater constancy of the tempera- 
ture and of the speed of rotation of the outer cylinder. The 
mean of 54 determinations, carried out with the improved 
apparatus, is 


nyo? = (18201+27)-10-8 or o3°= (18349427) -10-8 


in good agreement with the previous value, and with 
somewhat greater accuracy. 





ERRATUM 


B. TAYLOR 
If we use Millikan’s data with this value of » we would 


find for e 


e= 


18349\! 

( 4.770-10-" = (4.818+0.011)-10~ e.s.u. 
18226 

the uncertainty stated being due only to the viscosity. 

It should be emphasized, that the error in » has been 
calculated so as to give a reasonable /imit of error rather 
than a probable error, as some systematic errors have also 
been considered. (The probable accidental error, computed 
by the usual least square method from the individual 
deviations would be only 1.5-10-* or 0.008 percent.) 

From this result the conclusion can be drawn, that the 
lower value of e obtained by the oil drop method of 
Millikan can very well be explained by the error in the 
n-value used. New measurements by Backlin and Flemberg® 
seem to indicate that other errors may also exist. 

A detailed report has been sent for publication elsewhere. 


GUNNAR KELLSTROM 
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University of Uppsala, Sweden, 
June 20, 1936. 
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Erratum: Heat Conductivity of Tungsten and the Cooling Effects of Leads 
upon Filaments at Low Temperatures 


IRVING LANGMUIR AND JOHN BRADSHAW TAYLOR, General Electric Company 
(Phys. Rev. 50, 68 (1936)) 


In the last sentence of the abstract the equation was misprinted. The correct equation is: 


log A\=0.9518—0.30 log T. 








